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XXVII. BAND ERSTES HEFT 1959 


Beitrage zur Berechnung von Kaplanturbinen 
Von J. Raabe 


1. Einleitung. Die zunachst empirische Entwicklung der Kaplanturbine ergab infolge der 
nahezu gleichzeitigen Ausbildung der Theorie des gebundenen Wirbelfadens, welche von L. Prandtl 
ausging, einen der Priifsteine fiir diese Methode. Durch die Arbeiten von W. Birnbaum! aus dem 
Jahre 1923 angeregt, gelang es erstmals M. Schilhansl? 1927, diese Singularitatenmethode auf das 
einfache Modell einer koaxial zylindrischen Laufradstrémung anzuwenden. Die experimentellen 
Untersuchungen von K. Hahn*im Jahre 1933 an Kaplanlaufradern verschiedener Fliigelzahl sicherten 
den Anwendungsbereich dieses Modells. P. Ruden* verfeinerte 1937 die WVonseeiiing von der Lauf- 
radstrémung im Zusammenhang mit Versuchen an Verdichterlaufradern von verschiedener Form und 
Fliigelzahl, die an Hand einer 1931 von A. Betz® entwickelten Methode unter zusiitzlicher Beriick- 
sichtigung der Verdrangungswirkung der Fliigel und unter Beachtung der freien Wirbel und der 
Eulerschen Gleichungen gedeutet wurden. 1951 gelang es M. Strscheletzky®, die fiir die Kaplanturbine 
typische Strémung in Leitrad und Ubergangsraum in Ubereinstimmung mit Versuchen zu erklaren. 

Die vorliegende Arbeit’, die sich auf Versuche des Verfassers an einem nach der Singularitaten- 
methode entwickelten Laufrad stiitzt, beriicksichtigt vor allem die Vertraglichkeit der Strémung 
in den Querschnitten aufeinanderfolgender Bauteile (Leitrad, Laufrad, Saugrohr) beim Entwurf. 
AuBerdem bringt sie methodische und grundsatzliche Erganzungen, die sich auf die angeniherte 
zwei- und dreidimensionale Behandlung der Strémung sowie die Berechnung der Kennlinien einer 
Teilturbine beziehen. 


2. Die ebene inkompressible Strémung durch das gerade Gitter. Nimmt man innerhalb des 
Laufrades einer Kaplanturbine rotationssymmetrische, aquidistante Stromflachen an, dann kann 
man bekanntlich dessen Durchstrémung auf das Modell der ebenen, wirbelfreien Strémung durch 
ein gerades Gitter aus Profilen zuriickfiihren. Auf er der Teilung t und der Profillange L ergeben 
sich aus der ebenen Gittertheorie® als wesentliche Parameter des Strémungsvorganges der Win- 
kel £9, den die. auftriebslose Anstrémrichtung mit der Gitterachse einschlieBt, und der Gitter- 
verstarkungsfaktor K, welcher definiert ist durch 


rota (SL: (loom: ‘ 
6,—> 0 9 /Gitler ‘0 / Hinzelftiigel 
Dabei bezeichnet 6, den Nullanstellwinkel, den die um f,, gegen die Gitterachse (gleich Umfangs- 
richtung) geneigte Anstrémgeschwindigkeit w.. mit der auftriebslosen Anstrémrichtung einschliebt 
und c, den Auftriebsbeiwert. 

Es werde im folgenden eine Strémung von geringer Reibung angenommen, wobei ein kleiner 
Gleitwinkel ¢ vorliegt und wobei die obigen Gitterkennwerte mit guter Naherung mittels bekannter 
Methoden® aus der ebenen Potentialstré6mung entnommen werden kénnen. Bei den hier betrach- 


1 W. Birnbaum, Z. angew. Math. Mech. 3 (1923), S. 290. 

2 M. Schilhansl, Jahrb. wiss. Ges. fiir Luftfahrt, Miinchen (1927), 5. 151. 

3 K. Hahn, Untersuchung der Stromung durch eine Fliigelradturbine bei verschiedenen Schaufelzahlen, 
Dissertation Karlsruhe 1935. 

4 P, Ruden, Luftfahrtforschung 14 (1937), S. 327. 

5 A. Betz, Ing.-Arch. 2 (1931), S. 359. 

6 M. Strscheletzky, Ing.-Arch. 19 (1951), S. 309. 

7 Die Arbeit ist eine Kurzfassung der von der Fakultat fiir Maschinenwesen und Elektrotechnik der TH 
Miinchen genehmigten und am 13. 7. 1955 eingereichten Habilitationsschrift des Verfassers. Der Verfasser 
dankt an dieser Stelle Herrn Prof. Dr.-Ing. K. Hahn fiir seine wohlwollende, stetige Unterstiitzung, die cinen 
AbschluB dieser Arbeit erst ermoglicht hat. 

8 Siehe auBer den Arbeiten von Birnbaum, Betz und Schilhansl die folgenden Arbeiten: E. Konig, Z. angew. 
Math, Mech. 2 (1922), S. 422; F. Weinig, Strémung in Schaufeln von Turbomaschinen, Leipzig 1935; E. Pisto- 
lesi, L’Aerotecnica, 17 (1937) S. J. Ackeret, Schweizer Bauztg. 120 (1942) S. 120; H. Schlichting, Berech- 
nung der reibungslosen, inkompressiblen Strémung fiir ein vorgegebenes Schaufelgitter, VDI Forschungs- 
heft 447 (1955). 
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teten diinnen Profilen kann der Dickeneinflu8 entweder vernachlassigt oder durch die in Abschnitt 3 
dargestellte Naherungsmethode bestimmt werden. Mit dem oben definierten Nullanstellwinkel 
Boo —Bpo» einem die Zirkulationsminderung in wirklicher Stroémung beriicksichtigenden Beiwert 
m <2, welcher annahernd aus der Kennlinie des Einzelprofils in ebener Str6mung entnommen 
werden kann, und dem Gitterverstarkungsfaktor K ergibt sich bekanntlich der Auftriebsbeiwert 
im Gitterverband zu 

Cig KGsin (Pes 6)e (2) 
Fiir die hier betrachteten kleinen Anstellwinkel darf man! annehmen, da K nur von dem Winkel der 
auftriebslosen Anstrémrichtung gegen die Gitterachse /,,) und dem Teilungsverhdltnis t/L abhangt: 


K = K (6,9, 4/L) - (3) 


Gleichung (2) kann man in der folgenden Form 
schreiben: 


c, =m K (sin Bao 608 Bg — 08 Boo SiN Byo)- (4) 


Mittels der bekannten Gitterbemessungsgleichung 
fir Turbinen? 


12 ij Acu 
ae eT, Wo (l—tg € ctg Bo) (5) 


Cus und den beiden aus den Geschwindigkeitsdrei- 
Abb. 1. Geschwindigkeitsdreiecke. ecken (Abb. 1) zu entnehmenden Beziehungen 
Wag S10 Boo) == Cnr (6) 
1 
We COS Pog = U— Cm Ctg a + ai Ac, (7) 


kann man ctg Bx auf cp, u, Ac, und x, zuriickfiihren. Auf diese Weise erhalt man eine Gleichung 
zweiten Grades fiir die Ablenkungsgeschwindigkeit Ac,, deren Wurzeln mit der Abkiirzung 


At 
ae eae = (8) 


die folgende Form haben: 
Ac, = [ctg 8,9 + ctge (1 + %) + 2 ctg oa] o, —2u 


(+4) /c?, {[etg Boot ctg e (1+ )|?-+ 4 ctg « (x ctg x, — ctg Bn o)} — 4% Cte € Cn tt « (9) 


Da Ac, wesentlich positiv ist, gilt nur das obere Vorzeichen vor der Wurzel. 

Formel (9) stellt die Ablenkungsgeschwindigkeit Ac, als Funktion der DurchfluBgeschwindig- 
keit cm, der Umfangsgeschwindigkeit des Gitters u, der Gitterparameter /,,) und t/L und des Gleit- 
winkels ¢ dar. Da die bei den iiblichen kleinen Gleitwinkeln sehr groBen Glieder mit ctg é sich 
nahezu aufheben, kann der Einflu8 des Gleitwinkels nur bei sehr genauer numerischer Auswertung 
von (9) erfaSt werden. Man entwickelt daher zweckmafig die Wurzel von (9) nach Ausklammerung 
der Glieder mit ctg ¢ in eine binomische Reihe, deren Summanden wieder in geometrische Reihen 
zu zerlegen sind. Auf diese Weise bekommt man mit den Abkirzungen 


g=1+x, (10a) 


u 


x (ets on Se ce (10 b) 


die folgende Abhangigkeit der Ablenkungsgeschwindigkeit von den Potenzen von tg ¢: 


f=4 


coon 1/2 
=, o — zl Le ey = n+ tol n—m ot —n— m 
Ac, = 2 ctg 0&4 Gm —2u + cm i = (—1)"4 ileal tg” € ctg Boog Lf OTe da) 
Da in (11a) eine alternierende Reihe der Potenzen von tg vorliegt, ist der Fehler, den man bei 
Vernachlassigung des Gleitwinkeleinflusses begeht, kleiner als das im allgemeinen sehr kleine, in 
tg e lineare Glied der Reihe. Weil man den Einflu8 des Gleitwinkels auch durch geeignete Wahl 
des Beiwertes m (2) beriicksichtigen kann, wird im folgenden der Zusammenhang (lla) stets auf 


* K. Panteli, Wasserkraft u. Wasserwirtschaft 28 (1933) S. 241. 
> W. Bauersfeld, Z. VDI 66 (1922) S. 461. 
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den eines verschwindenden Gleitwinkels ¢ spezialisiert. Man erhalt hierfiir aus (lla) die Be- 
ziehung 


2 
AG a a [(ctg a, + ctg Boe) Cm — U] (11b) 


welche fiir den Sonderfall des geraden Streckengitters in ebener Potentialstrémung als strenge 
Lésung bekannt ist.1 : - 

Mit Beachtung von (3) und (8) stellt (11) bei vorgegebenen Geschwindigkeitsdreiecken eine 
Beziehung zwischen der durch f,9 gekennzeichneten auftriebslosen Anstrémrichtung des Gitters 
und dem Teilungsverhaltnis t/I, dar. 

Wie spater gezeigt werden wird, ist fiir den Zylinderschnitt des Fligels entweder t/L oder 
Bpo vorgeschrieben, Damit kann bei gegebenen Entwurfsdaten, wie dies in Abschnitt 7b angegeben 
ist, Gleichung (11) zur Bestimmung von fo bzw. t/L verwendet werden. In den Gleenien, Ab- 
schnitten soll gezeigt werden, wie man aus t/L, 6,9 und den Geschwindigkeitsdreiecken ein még- 
lichst kavitationssicheres Profil und dessen Staffelungswinkel 8, welchen die Skelettsehne des 
Profils mit der Gitterachse einschlieBt, erhalt. 


3. Die genaherte Bestimmung der auftriebslosen Anstrémrichtung im Gitterverband bei diinnen, 
schwach gewolbten Profilen. Bei sehr diinnen Profilen ist der Winkel, welchen die auftriebslose 
Anstrémrichtung und die Skelettsehne miteinander bilden, nur von den Wélbungsparametern des 
Skeletts und den Gitterparametern abhangig.? 

Bei maQig diinnen Profilen kann man in erster Naherung den Dickeneinflu8 dadurch beriick- 
sichtigen, dab man den Winkel Byo, den die auftriebslose Anstrémrichtung des Skeletts im Gitter- 
verband mit der Gitterachse einschlieSt, um einen Betrag 


iL 
Ad = — K (y, By + % By) (12) 


vermehrt. Hierin hangen B, und B,, wie im folgenden gezeigt werden wird, gemaB Abb. 2 und 
Abb. 3 von den Gitterparametern f,9 und t/L ab. Die Winkel y, und x, werden gemaB Abb. 4 


| 


G05) 
Gitferachse | 


Og 10 15 20 25 30 


Abb. 2. Verlauf der EinfuBfunktion B, (t/L, Bo) und Erliiuterungsskizze zu den GitterkenngréBen. 


von den gegeniiberliegenden Seiten eines Polygons gebildet, welches dem Profiltropfen einbe- 
schrieben ist. Durch (12) wird auch der Einflu8 des endlichen Hinterkantenwinkels 7, (Abb. 4) 
erfaBt, da dieser infolge der SchlieBbedingung fiir den Polygonzug wie folgt mit den beiden in 
(12) enthaltenen Winkeln y, und y. zusammenhangt: 


Ye i) Oy ON (13) 


1 Siehe die Arbeit von Weinig in Fufnote 8 von Seite 1. 
2 Sjehe die Arbeit von Pantell in FuBnote 1 von Seite 2. 
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Ersetzt man den Profiltropfen diinner Profile mit rektifizierter Mittellinie analog Abb.4 durch 
einen Polygonzug, dessen gegeniiberliegende Seitenpaare die Winkel y; einschlieBen, dann be- 
kommt man fiir die Neigungsanderung der auftriebslosen Richtung infolge endlicher Profildicke s; 
die zu (12) analoge, allgemeine Bezichung 

i 
apes aay 14 
IS) NSS (14) 
worin K und B; naherungsweise nur von den Parametern eines aus den Skeletten der Profile ge- 


bildeten Gitters abhangen. Zu (14) gelangt man, wenn man zur Erzeugung des Profils jeden der 
n durch die Polygonabschnitte begrenzten 
E;-Bereiche mit Einheitsquellen von jeweils 


konstanter Starke q; belegt, welche bekannt- 
lich mit y;, der Polygondicke s;, ~ s; und 
@ der ungestérten Anstrémgeschwindigkeit w., 
naherungsweise durch die Beziehungen 
p= heey 7 GIO) 15 
— q . is (15) 
- 
ss te 
Kia dé; (16) 
02 verbunden sind. Der aus (16) folgende Ge- 
radenzug schlieBt sich dann zum Polygon, 
wenn man die Neigungswinkel y; gegeniiber- 
liegender Seiten der Bedingung 
015 
n—l1 
Dy Asi 0 (17) 
i= 
Ot unterwirft, die sich fiir die gewahlten drei 
Abschnitte zur obigen Gleichung (13) spe- 
zialisiert. 
Betrachtet man ein gerades Profilgitter 
ua vom Staffelungswinkel 6 ~ 6,0 und der 
° o5 


Abb. 3. Verlauf der EinfluBfunktion B, (t/L, Bp o Abb. 4. Ersatz des Profiltropfens durch ein Polygon. 


Teilung t, und belegt man die i-Abschnitte aller Skelette mit Einheitsquellen q;, dann erzeugen 
diese am Punkt € =, von der Sehne des Aufskeletts normal zu dieser bekanntlich! die Ge- 
schwindigkeit 


3 l cos 6 Gin 2° sin] — sin sin 2a conf 
fodi — 9 {== are : ES d& (18) 
Coj a oe in p — cos 2 7 29 >! cos B 
Siu ’ : 
Mit q; aus (15) kann man Vs, ,,, geschlossen darstellen als 
Tq | (& — Ein) si 2 ae | 
4 5 Woo Xe g i 0 iu) sin B Tq (&5 €io) sin B 
ber (ag ee AEC CON eee ——— - | — are tg ———— «nt L9) 


ae ; 
tg F (§) — &iu) cos B tg 


| — (Go — Ei) cos 5 


wobei &€;, —&;, = AE; die Erstreckung des i-ten Polygonabschnittes in €-Richtung ist. 


* NN. Scholz, Strémungsuntersuchung von Schaufelgittern, VD1-Forschungsheft 442 (1954) S, 29 


XXVII. Band 1959 | J. Raabe: Beitrage zur Berechnung von Kaplanturbinen BY 


Diese mit einer tangentialen Umstrémung des Skeletts (bzw. Profils) unvertragliche Geschwin- 
digkeit erfordert eine zusatzliche Belegung der Skelette mit Einheitswirbeln der, Starke y,, deren 


Induktion normal zur Skelettsehne im Punkt & die Geschwindigkeit by Veg; gerade aufhebt. 


Es muB also gelten 


er anes a= 0), (20) 


Wahlt man fiir y, die ersten drei Glieder der Birnbaumschen Reihe, dann kann der Bedingung 
(20) an drei Aufpunkten geniigt werden. Es seien hierfiir die von Birnbaum benutzten Punkte 
gewahlt'. Dabei wird nach dem Vorgang von Pantell? naiherungsweise 8, durch Ppo ersetat. Fir 
die Freiwerte a, b, c der ersten drei Glieder der Reihe von yq erhalt man so drei lineare Gleichungen. 
Damit ergibt sich bei Beachtung der SchlieBSbedingung (17) und mit den bekannten Beziehungen 


=a(2atb+4), (21a) 
AS Ca 
A Sea (21b) 


der in (14) angegebene Zusammenhang, welcher sich fir den hier behandelten aus drei Seiten- 
paaren bestehenden Polygonzug zu (12) spezialisiert®. 


4. Die genaherte Bestimmung der fiir die Kavitationssicherheit optimalen Schaufelform beim 
geraden Gitter in ebener Potentialstromung. Nach der Bemerkung am SchluB des Abschnittes 2 
ist es wichtig, bei einem Gitter von gegebenen Parametern t/L und f,o in einer Strémung von 
gegebenen Geschwindigkeitsdreiecken die kavitationssicherste Fliigelform zu finden. Sie folgt aus 
der Forderung nach konstantem Druckverlauf auf der Saugseite des Fliigels. Die Aufgabe werde 
gelést unter der Annahme, dab es 


sich um diinne, schwach gewélbte : a6 
Profile von kleinem Anstellwinkel +g : ee 
handelt. 


Wenn die Wirbelbelegung y 
und die Quellenbelegung q vom 
Restgitter* auf der Sehne des Auf- 
profils die Langsgeschwindigkeiten 
Ve,,r und VzgRr erzeugen, dann 
fiihrt die gewiinschte konstante 
Geschwindigkeit langs der Saug- 
seite w, auf die bei kleinem Anstell- 
winkel giiltige Beziehung® 


y 
Ww, = Wao + on +- Vey,R 


+ Veq,r = konst. (22) 


Abb. 5. Modell eines geraden Profilgitters zur angendherten Bestimmung 
der vom Restgitter an einem Aufpunkt erzeugten Geschwindigkeiten. 


Zur gendherten Bestimmung 
der Langsgeschwindigkeiten durch 
das Restgitter Vz,,,p und V:4,n diene ein in Abb. 5 dargestelltes Gittermodell. Bei diesem wird 


die Verdrangungswirkung eines Profils von der maximalen Dicke s,,,, in einer Strémung von der 
ungestérten Anstrémgeschwindigkeit w.. durch eine Quelle und Senke von der Starke 


1 Siehe die Arbeit von Birnbaum in FuBnote | auf Seite 1. 

2 Siehe die Arbeit von Pantell in FuBnote | auf Seite 2. 

3 Die formelmaBige Abhangigkeit der Gré®en B, und B, von den Gitterparametern, welche aus der vor- 
hergehenden Darstellung ohne Schwierigkeit folgt, ist hier wegen ihrer Weitlaufigkeit unterdriickt. 

4 Unter Restgitter werden alle Fliigel eines Gitters verstanden mit Ausnahme des sogenannten Auffliigels, 
welcher den Aufpunkt enthalt. 

5 Siehe die Arbeit von Schilhansl in FuBnote 2 auf Seite 1. 
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erzeugt, welche im Abstand + 1/2 und — uw L/2* vom Ursprung & = 0 (in der Mitte der Skelett- 
sehne) liegen. Die Zirkulation um ein Profil, welche bekanntlich mit der Ablenkungsgeschwindig- 
keit Ac,, der Teilung t und der Wirbelbelegung y in dem Zusammenhang 

+ L/2 

Pe | vdea ews (24) 

Sip 
steht, soll durch einen Stabwirbel verursacht werden, welcher im Druckmittelpunkt des Profils 
an der Stelle £ =€, angeordnet ist. Dabei falle die -Achse gemaB Abb. 5 in die Richtung der 
um den Winkel 8 gegen die Gitterachse geneigten Skelettsehne. 

Entwickelt man die konjugierte Geschwindigkeit? der Singularitatenanordnung von Abb. 5 in 
eine Potenzreihe von é, dann ergibt bekanntlich das erste zu € umgekehrt proportionale Glied die 
Wirkung der einzelnen Singularitat des Aufprofils wieder. Das folgende in & lineare Glied liefert 
dann in erster Naherung die Wirkung der iibrigen Singularitaten. Auf diese Weise gewinnt man 
die Langsinduktion des Restgitters zu 

m (§ 


ee 
Venn + Vepn=— =e Ac = w cos 2p. (25) 


Um die Gleichung (22) zu erfiillen, hebt man zweckmafig das in € langs der Fliigeloberflache 
veranderliche Glied der letzten Gleichung (25) durch den Ansatz 


g 
Y=Vo bles (26) 
auf. 
Uber die Beziehung 
; + L/2 
i= wr | yé dé (27) 
Sap 


bekommt man mit (24) und (26) die folgenden Gleichungen fiir y), y, und €,: 


t 
Yo a Ac,, ns ? (28) 
ae Lb 7 
a Ac, sin 2B, (29) 
ae ff IGP ; 
£=— 3% (7) Lome. (30) 
Damit wird die saugseitige Geschwindigkeit gemaB (22) 
= mak Srey IC t ae ff ILXS 
w, tao (1 bag cos 2) + Ae, [55 4 ato (+) sint2 6). (31) 
: 2 2 5 o 
Mit 205 praae UO) 3 ergibt sich auch der in der folgenden Bestimmungsgleichung fiir den 
Thomaschen Kavitationsbeiwert o bei drallfreiem Austritt 
a we 
+) _2 (32) 


Oo = sm 
iu 2gH 2g H 


1 Legt man Quelle und Senke in den Schwer < 
; gt punkt der Quellen- und g i 1 
wird bei einer den GroBausfiihrungen naherungsweise par er merynre cee ee Eee 


3 
iE re Se 
'S == AV Sma 7 L 


gekennzeichneten Profiltropfenform der Wert = te é 
3 


2 Siehe et : : 
ee ee etwa W. Kaufmann, Technische Hydro- und Aeromechanik, S. 298 ff, Berlin/Géttingen/Heidel- 
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enthaltene Beiwert 2 an dem betrachteten Profil zu 


mallee, i|'Acy : 
(HY sine p | | Greet Sinan Fe cos B + 1 


= ue max a= ll PN: eo 2L Woo 9 t t 
a Gia sane amen mo 
aa ae — ri cos fp, +1 


in Abhangigkeit von den Profil- und Gitterparametern fb, t/L, Smax wnd den GréBen Ac,, w,, und 
Boo der Geschwindigkeitsdreiecke. Bei gegebener Ortshdhe h, Saughohe h,, Fallhéhe H cael Dae 
druckhéhe hy erfordert der kavitationsfreie Betrieb bekanntlich einen Kavitationsbeiwert von 
héchstens Om. = H-1 [10,33 — h, — 0,0012 h — h,] (wobei alle Héhen in [m] vorausgesetzt 
werden). Dadurch ist mit den Geschwindigkeitsdreiecken und dem Saugrohrwirkungsgrad Ps 
der Wert A gegeben. Liegen s,,,,:t und f fiir einen Zylinderschnitt fest, dann folgt das Uberdeckungs- 
verhaltnis L:t bei einem Fliigel von konstantem saugseitigen Druck mit Riicksicht auf Kavi- 
tationssicherheit aus Gleichung (33a), welche niherungsweise in der einfachen Form 


t Ac, 2 
Se E tat 1 1 (33b) 
verwendet werden kann. 


Die mit der Wirbelbelegung y nach (26) verbundene Skelettform 7(&) wird bei y <w,, und 
unter den Voraussetzungen von Birnbaum! 


Dime Acu ate | 1 oy ME We mw /L\2 /2 &\2 Ge? ff du SPAS 
fa L ea ! 13 (7) 8im28 Fa] — ay (7) (T} 09526 a (7) 1 in| 


TEES mace ones 
eae ; op in 2B. (34) 


Hierin sind F, und F, die in Tabelle 1 tabellierten folgenden Funktionen: 


ee b+ | in(1 ! =} y ! L =| In —y\a +e, (35) 
2é 
; (epi 
2 é\2 ip 2€ 
F, = L— (>) In mae +25. (36) 
siald by 


Tabelle 1. Verlauf der Funktionen F, und F, tiber 2 €/L 


2 | | 

oe WWE 45 Ws Wi 20-4 | ee eee | 22) 1 cael 
Be 40,088 |) '=-90,786 493865 | +0403. | £0,041 +0 
Fr | — 23605 | —2392 | —2088 | —1's12 | —o7789 +0 


F, ist eine gerade, F’, eine ungerade Funktion in &, 

Die in (34) angegebene Lésung fiir das Skelett mit konstantem Unterdruck auf der Saugseite 
gilt nicht an den Skelettenden. Man kann diesen Ubelstand durch Extrapolation an den Enden 
beseitigen. ErfahrungsgemaB ist der Fehler dabei nicht gro. Am Profilkopf versagt ja ohnedies 
infolge der stets auftretenden starken Kriimmung der Profiloberflache die einfache Bedingungs- 
gleichung (22) fiir konstante saugseitige Geschwindigkeit?. 

Die nach der obigen Methode ausgelegten Profile weisen auf Grund von vergleichenden Be- 
rechnungen mit einer genaueren Methode® eine angenahert konstante Geschwindigkeitsverteilung 
auf der Saugseite auf. 

Als Beispiel ergeben sich fiir ein Gitter mit Bpo = 20,5°, L/t =0,86 und ein durch ¢:w. = 
0,287 sowie Bx = 27° gekennzeichnetes Ablenkungsdreieck bei = 2/3 die in der Abb. 6 dar- 
26 


L° Man erkennt, daB das Einzelskelett von konstantem 


gestellten Verlaufe des Skeletts tiber 


1 Siehe die Arbeit von Birnbaum in FuBnote 1 auf Seite 1. 
2 Um Geschwindigkeitsspitzen an der Eintrittskante, welche zu Kavitation fiihren konnen, zu unterdricken, 
wird man diese entsprechend abrunden und so ausbilden, daf} der eintrittsseitige Staupunkt in die Mitte des 


Bogens von der Profilnase fallt. 
3 Siehe die Arbeit von Schlichting in FuBnote 8 auf Seite 1. 
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saugseitigen Unterdruck symmetrisch zur 7-Achse verlauft, daB das Skelett im geraden Gitter- 
verband von unendlich diinnen Profilen einen S-Schlag hat und da® das Skelett im Gitterverband 
von diinnen Profilen einen dagegen etwas geringeren S-Schlag aufweist. 


7] 
Te 


Abb. 6. Skelettformen bei angenihert konstantem Druck auf der Saugseite. 


5. Die Abhingigkeit der Meridiangeschwindigkeit von der Drehzahl bei einer Teilturbine von 
gegebener Fliigelform, Leitradstellung und gegebenen Verlustbeiwerten bei zweidimensionaler Be- 
trachtung. Unter der Annahme rotationssymmetrischer Stromflachen innerhalb der aus Leit- 
paparat, Laufrad und Saugrohr gedachten Turbine und unter Voraussetzung koaxial zylindrischer 
Stromflachen (und damit ebener stationarer Potentialstrémung im Laufrad) soll untersucht werden, 


wie sich die auf die AusfluBgeschwindigkeit bei der Fallhéhe H bezogene Meridiangeschwindig- 


keit ¢m,, = Cm! / 2g H von einer zwischen zwei benachbarten Stromflachen liegenden Teilturbine 


mit der bezogenen Umfangsgeschwindigkeit u;g = u:|/ 2 gH des zugehérigen Teilfliigels andert, 
wenn fiir diese Teilturbine der vor allem von der Leitradstellung abhangige absolute Zustrém- 
winkel zum Laufrad «, (Abb. 1), der auftriebslose Anstrémwinkel 6,9 vom Laufradgitter, dessen 
Teilungsverhaltnis t/L und Gleitwinkel ¢ und die Saugrohrverlustbeiwerte ¢,,, und 4m (welche den 
Druckverlust im Saugrohr zu h = Cm ¢m/2 g + Cu Cu2/ 2g bestimmen) gegeben sind. 
Beriicksichtigt man den Einflu8 der Reibung durch geeignete Wahl des Beiwertes m in (2), 


dann ist der Zusammenhang zwischen der bezogenen Umlenkungsgeschwindigkeit Ac,,, = Ac, : / 2gH 
und ¢,,, gemaB (11b) durch 
2 
AG == ant [(ctg a, + ctg Bpo) ¢m,, — via] (37) 
gegeben. 
Die Turbinenhauptgleichung lautet in bezogenen Geschwindigkeiten 
2 u,, Ae,.4 = Mr - (38) 


Hierin ist bekanntlich’ der hydraulische Wirkungsgrad 7, in erster Naherung in nachstehender 
Weise von den bezogenen Geschwindigkeiten ujq, Aeu,,> €m,q> der Zustrémrichtung x, , dem Gleit- 
winkel e und den Saugrohrverlustheiwerten ¢,,, und C,, abhangig: 

Veg ee gee ae perenne oly 

Th Cin sm “mi su Cmiq cls &, — Cu; 4) . (39) 


Damit ergibt sich aus den letzten beiden Gleichungen 


Cs u ar Cid : oe oe osu i af } 
' C. Keller und H. Bleuler, Escher-Wyss-Mitteilungen 10 (1937) S. 85, 


Psu 


Usa a l a ea, ; 2 etg ot, ¢ seats ae Ce. ; ue 
Aeu;, on eer : -L ctg Oy Cae fis a xt - mid id sm “My y “id (40) 
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Durch Elimination von Ac,,, aus (37) und (40) bekommt man die gesuchte Abhangigkeit in der 
Form 


3 2 2 
A Cmzq + 2 Bere, Wid — C Om, , Wid — Cm, —- € Wig = 0. (41) 


Dabei sind die Beiwerte A, B und C in nachstehender Weise durch die Parameter der Teilturbine 
bestimmt: 


A= Con alr es [ctg == N (ctg Oy -+-+ ctg Bro) |? : . (42) 
B = N [(ctg Oy = ctg Bpo) (1 oan N) pe Gast ctg | 9 (43) 
C2 Nee (44) 
mit 
2 
== ease (45) 


wobei x aus (8) folgt. Fiir ¢ = 0 erhalt man die bezogene Meridiangeschwindigkeit zu 


B Bie Ce as 
Cmia,0 a A’ U4 as yl ra ia A Utd ae 2 (46) 
wahrend fiir kleine Werte von ¢ aus (46) die angenaherte Beziehung 


C =p € ae 
Ait eA 2 Cmid, 0 (A Cmid,o B ug) 


(47) 


folgt. 
Da der Gleitwinkel im allgemeinen klein ist und dessen Beiwert in (47) etwa die GréBe der Einheit 
hat, ist der KinfluB des Gleitwinkels auf den Verlauf der Meridiangeschwindigkeit tiber der Um- 


a,b, A By?-10' 


Hyperbel 


Cmid,o 


—— 
—— 
——— 


Abb. 7. Verlauf von y, A, B, b, a iiber a, bei L/t = 0,9 Cn = 0,2, oper = 1,0 und schematischer Verlauf von ‘mig, (usq) bzw. Smid, (44): 
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fangsgeschwindigkeit gering. An Stelle der ¢,,,, (u,,)-Linie werde daher die den GleitwinkeleinfluB 


nicht enthaltende ¢n;,,. (u,,)-Linie diskutiert. Letztere ist mit (46) als eine Hyperbel' gegeben, 


deren Halbachsen 


a = (A cos?y + Bsin2 y — Csin® y) 1? (48) 
b = (Asin? y + Bsin2 y + Ccos*? y) 1? (49) 

um den Winkel 
y => are tg qi (50) 


im Gegenzeigersinn zur u,,- bzw. Cm;4,97Achse gedreht sind, wie dies Abb. 7 erlautert. 


Tabelle 2. Die Gréfen A, B,y, a und b in Abhdngigkeit von bei Csm = 0,2 und Csy = 1,0 


fiir Bpo = 20° 50° | 1,237| 1,278 | 24,00 | 0,782 | 0,744) 0,7 
20° | 0,9113| 3,182 | 36,72 | 0,546 | 0,553 | 0,7 a ee | ne Pe nae Ree He 
20° | 3,036 | 3,360 | 29,14 | 0,579 | 0,452 | 0,9 | See Bs Asearat ye 797 0, : 

: 9 ; 2 | | 
RS eee eae cae eed ate 60° | 1,3480) 1,038 | 20,35 | 0,831 | 0,760 | 0,7 
30° | 1,496 | 2,323 | 30,96 | 0,649 | 0,578 | 0,7 60° | 1,964 | 1,001 16,36 | 0,833 | 0,666 | 0,9 
30° | 3,530 | 2,398 | 23,01 | 0,696 | 0,470 | 0,9 60° | 2,291 | 0,966 | 14,53 | 0,847 0,630} 11 | 
30° | 4,930 | 2,358 | 18,85 | 0,712 | 0,418 | 1,1 ' 
i = 40° | 
40° | 1,886 | 1,864 | 26,22 | 0,726 | 0,598 Sir B po | 


S 
un 
uw 
nS 
= 
. 

fo 
on 
— 


0,7 é 
40° | 3,872 | 1,884 | 18,53 | 0,765 | 0,472 | 0,9 oh ened ae 38,32 
40° |5,068 | 1,826 | 15,17 | 0,775 | 0,424) 1,1 30°. 0-615. 2.83517 33 kOe aeuo eee 


50° | 2,165 | 1,565 | 22,41 | 0,784 | 0,599| 0 
50° | 4,008 | 1,551 | 15,63 | 0,803 | 0,475 | 0 
1 


: 30° | 0,2182) 1,783 | 34,86 | 0,654 | 0,828 
50° |5,170 | 1,479 | 12,50 | 0,818 | 0,427 1 


0,7 
30° | 0,3108) 1,833 34,09 | 0,646 | 0,803 | 0,9 
30° | 0,4048) 1,845 | 33,34 | 0,646 | 0,786 | 1,1 
60° | 2,390 | 1,344 | 19,62 | 0,830 | 0,593] 0 
60° | 4,152 | 1,303 | 13,15 | 0,839 | 0,474| 0 
60° | 5,230 | 1,223] 10,45 | 0,840 | 0,428 | 1 


w 
— 
_ 
Oo 
S 
ae 
Ww 
bo 
S 
oO 
lo) 
Ne) 


7 
9 40° | 0,2504) 1,279 0,7 
1 40° | 0,3412| 1,306 | 30,08 | 0,723 | 0,954] 0,9 

40° | 0,4187| 1,310 | 29,26 | 0,722 | 0,933 1,1 


far’ Bip = 30° 50° | 0,2802/ 0,949 | 27,08 | 0,790 | 1,131 | 0,7 
20° | 0,5875| 3,015 | 37,34 | 0,545 | 0,588 | 0,7 30° | 0,3636) 0,962 | 25,84 | 0,785 | 1,099 | 0,9 
20° cc 3,105 | 33,63 | 0,550 | 0,550) 0,9 50° | 0,4278) 0,959 | 24,96 | 0,776 | 1,079 | 1,1 
20° | 2,019 | 3,115 | 31,40 | 0,569 | 0,505 | 1,1 


60° | 0,3063} 0,704 | 22,53 | 0,848 | 1,294 | 0,7 
30° | 0,8920| 2,086 | 32,44 | 0,652 | 0,658 | 0,7 60° | 0,3803} 0,706 | 21,27 | 0,832 | 1,235 |) 0,9 
30° |1,682 | 2,120 | 28,15 | 0,662 | 0,596 | 0,9 60° | 0,4348) 0,699 | 20,39 | 0,827 | 1,200} 1,1 
30° | 2,147 | 2,106 | 25,83 | 0,674 | 0,562] 1,1 


40° |1,094 | 1,599 | 27,99 | 0,715 | 0,709 | 0,7 
40° | 1,816 | 1,598 | 23,59 | 0,736 | 0,631 | 0,9 
40° | 2,215 | 1,573 | 21,38 | 0,745 | 0,595 | 1,1 


Tabelle 2 gibt die Abhangigkeit der Gré8en A, B, a,b und y von den Parametern einer Teil- 
turbine bei m = 2.7, €,, = 1,0 und €,m = 0,2 wieder. Abb. 7 veranschaulicht die Abhangigkeit 
dieser GréBen von x, bei L:t = 0,9. 

Aus dem Verlauf von ¢,,, (u,,), wie ihn Abb. 7 prinzipiell und Abb. 23 fiir ein bestimmtes 
Beispiel zeigen, ist zu erkennen, daB im allgemeinen Cm;, mit wachsendem u,, zunachst abnimmt 
und nach Erreichen eines Minimums wieder zunimmt. Wegen der bekannten Proportionalitat 
von ug mit der Einheitsdrehzahl n,, =n D H-'? (worin n die Drehzahl, D den Laufraddurchmesser 
und H die Fallhéhe bedeuten) und wegen der ebenfalls bekannten Proportionalitat von Cm;q mit 
dem Einheitsgesamtstrom Q,, =Q D-? H-1/? (worin Q der Gesamtstrom ist) ergibt sich aus dem 
besprochenen Cmn;, (Wia)-Verlauf die Erscheinung, daB die Q,,(n,,)-Linie mit wachsender Drehzahl 
steiler ansteigt. 


' Ein derartiger Zusammenhang ergibt sich auch bei Zugrundelegung der Stromfadentheorie in Verbindung 


mit den tblichen Ansiatzen fiir den StoBverlust, Laufradverlust und Diffusorverlust. Hieriiber siche C. Pflei- 
derer, Stromungsmaschinen, S, 186, Berlin, 1957. 
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6. Die Abhingigkeit der Meridiangeschwindigkeit von der Form des Leitapparats und dem Ge- 
samtstrom bei freien Wirbeln geringer Starke. Die achsensymmetrische Meridianpotentialstrémung 
durch einen Rotationshohlraum, dessen Randstromflachen denen des Leitapparates entsprechen, 
sei bekannt und durch ihre Meridianstromlinien sowie die vom gegebenen Gesamtstrom abhangige 
Geschwindigkeit ¢n,po: gekennzeichnet, die wirbelfrei angenommene Zustrémung zum Leitrad 
ebenfalls. Damit ist bei gegebener Leitradstellung und bestimmbaren Meridianpotentialstrom- 
linien in erster Naherung der Winkel x), den die Stromlinien an der Leitradaustrittskante mit 
der Umfangsrichtung einschlieBen, berechenbar.! Hieraus folgt im allgemeinen bei konstant iiber 
die Leitschaufeleintrittsbreite angenommenem Zustrémdrall ein langs der Leitradaustrittskante 
veranderlicher Drall, welcher mit der Meridiangeschwindigkeit der Potentialstrémung am Leitrad- 
austritt Cm p,9 anndhernd 2u ry ctg X Cm Pot,o Angesetzt werden kann. Infolge einer demzufolge 
langs der Schaufelbreite veranderlichen Zirkulation um eine Leitschaufel treten von deren Aus- 
trittskante freie Wirbel in die Strémung mit der Geschwindigkeit c, von welcher angenommen 
werde, da sie sich aus der durch die freien Wirbel induzierten Geschwindigkeit ) und einer Ge- 
schwindigkeit ¢p,,; zusammensetze, entsprechend 


C= Cre +. (51) 


Das durch ¢p,; gekennzeichnete Strémungsfeld verlaufe auf den Stromflachen der Meridianpotential- 
str6mung und habe axiale und radiale Geschwindigkeitskomponenten c, Por? Tr Pop? Cie aus der 
Zerlegung von Cm pot folgen und 
eine Umfangskomponente, die 
sich unter Voraussetzung kleiner 
Wirbelstarken tiber den Drallsatz 
aus dem Leitradaustrittsdrall zu 


Cu pop = 2 Cm pot, CY %q ergibt und 
damit auch durch die Meridian- 
potentialstrémung gegeben ist. Die 
zur Berechnung der Induktion p 
erforderliche Zirkulation d2/’ um 
das Element eines freien Wirbel- 
fadens kann aus den bekannten 


GréBen Cm pot, und Cr pop.g 2D der 


Leitradaustrittskante dann  be- 
rechnet werden, wenn man an- 
nahernd annimmt, daB die Wirbel- 
faden den durch ¢p,; beschriebenen 
Stromlinien folgen und sich, von 
einem sehr kleinen Abstand von 
der Leitradaustrittskante begin- 
nend, in Umfangsrichtung gleich- 
maBig iiber eine Potentialstrom- 
flache verteilen. 

Schneidet man gemafs Abb. 8 
aus der Stromschicht zwischen 


7 ° Abb. 8. Zur Induktion cines von der Lvitschaufelaustrittskante ausgehenden 
zwei benachbarten Potentialstrom- Wirbelfadenelements W an einem Aufpunkt P. 


flachen eine Stromréhre heraus, ; . 
welche zugleich Wirbelréhre ist, die an dem Punkt W(0,r, x) in Umfangsrichtung die Erstreckung 
+ d9 und normal zu den Stromflachen die Lange dn hat, dann ergibt sich die gesuchte Zirkulation 


um diese Stromréhre zu 


il d (Cm pot,g Ct& %o) 


"Crips? dn di. (52) 


Cr Pot,0 dx 


Dabei beziehen sich die mit dem Zeiger 0 versehenen GréBen gemals Abb. 8 auf den Schnittpunkt 0’ 
der durch den Punkt W gehenden Potentialstromlinie mit der Leitradaustrittskante. 


1 Siehe die Arbeit auf Strscheletzky in FuBnote 6 von Seite 1. 
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Sind 1°, t° und x° die Einheitsvektoren am Punkt W in tangentialer, radialer und axialer Rich- 
tung und ¢y 5,5 Cr po,> © pos die analogen Komponenten der Geschwindigkeit Cpo, dann wird der | 


Vektor eines zu Cp, parallelen Wirbelelementes im Punkte W 


i= je len cos } -+- sin 0) + r° be ee sin @ + cos ») ae a dr. (53) 


T Pot T Pot Cr Pot 


Der Abstandsvektor von einem Aufpunkt P (9 = 0, r = rp, x = xp) zum Punkt W(0, r, x) eines 
Wirbelelements ist 


a=u°rsin & + 1° (r cos # — rp) + £° (« — xp). (54) 


Die von einem mit freien Wirbeln erfiillten Gebiet G in P induzierte Geschwindigkeit folgt mit dem 
Biot Savartschen Gesetz zu 


oh 1 S GRIP Gail) tl f " dl [a, dl] : 1 PP aR oe aih 
t=a){{ Sar ae] || yp eee ee ets ea 
G 


lal 0 


Da das letzte Glied von (55) verschwindet, ergeben sich fiir die Komponenten der induzierten Ge- 
schwindigkeit in radialer, axialer und tangentialer Richtung nach Integration tiber den Bereich 
0<0< 2a die Ausdriicke 


ee 4 5 1 d(cm Pot,0 CS Oy) To Cm pot,9 CLE % (x — xp) J, dn ds 56 
ey | Cr Pot,0 dx pe Z eo) 
G 
: Laarae oe ad (cm por,o ctg %) To Cm pot, Ct & (Tp Jy —1r Jz) dn ds P 
oO Oh ee Cr Pot,0 dx jae are et’ Ce 
G 
FA 1 i r (Cm por,9 CtE %) {[r Cx pop — (% — Xp) Cr py] Jy — Tp Cx po, Ja} dn ds 
re Dies Frets dx p® Co 


mit s als Koordinate in Richtung der Meridianstromlinien und n als Koordinate normal zu diesen. 
Dabei bedeuten 


De aes ae Oe (59) 
fe a ey 
1—qVl+q ° (ce) 
2rrp 


a Ce (61) 


1/2 


= dp eae ae 
B | Weaeanoe © (vollstandiges elliptisches Integral erster Gattung), (62) 
0 


m/2 


f= i 1 — k*sin® y dg , (vollstandiges elliptisches Integral zweiter Gattung), (63) 


ale 
5 \/ Ta ; (64) 


PPT eke ne | (65) 


i eerie src ee) ae (58) erstrecken sich iiber das wirbelerfiillte Gebiet G in der 
i nter dem Leitrad und sind im allgemeinen nur dann eraphisch 1é i 
auf die Leitradaustrittskante folgende Begrenzung des Gebietes G Hee Lares 
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Laufradentwurf, wie iiblich die d i 
ur 5 as r ay . 
ea ee _ eee ch das Laufrad nicht gestérte Zustrémung vom Leitapparat 
S - ; R rstreckt sic von der Leitradaustrittskante stromwiarts theoretisch bis ins Un- 
en me 3 da den Gleichungen (56) bis (58) rotationssymmetrische Stromlinien zugrunde gelegt 
sind, ist G in der Praxis unter Umstanden durch den Eintrittsquerschnitt in ein nicht rotations- 
symmetrisch ausgebildetes Saugrohr zu begrenzen. 

7. Der Laufradentwurf bei zweidimensionaler Strémung. a) Die hinsichtlich Kavitations- 
sicherheit und Wirkungsgrad optimalen KenngréBen des auBeren Fligelschnittes 
Es werde die spezifische Drehzahl n, wegen der im allgemeinen vorgegebenen Anlagedaten, wie 

9 
der Fallhéhe H, dem Wasserstrom Q und der Drehzahl n als ein durch die Beziehung n, = n oa 
2 15 
-3/4 
cs /4 gegebener Parameter angesehen. Dasselbe gelte vom hydraulischen Wirkungsgrad 1, dort, wo 
os * . L 
seine Anderung vernachlassigbar ist. AuBerdem werde drallfreier Austritt angenommen. Der in 
(32) definierte Thomasche Kavitationsbeiwert o, welcher mit dem nahezu konstanten, im Wesent- 
lichen nur vom Nabenverhialtnis y abhangigen Beiwert : 


_ ag (yd — v2 
A 30 75a g/t "9/4 (66) 


wie folgt von der DurchfluBziffer ~ = c,,/u, dem Saugrohrwirkungsgrad sm = 1—.m und dem 
Beiwert A (33) abhangt 


ee Ns 4/3 Nh [A p23 ae (r EA) 4/3 
fi 9 ]sm YP ]> (67) 


nimmt unter der Optimalbedingung do/dp = 0 bei einer fiir die Kavitationssicherheit optimalen 


DurchfluBziffer 
= Wn 1/2 
Po opt a 2 (Ysm =| (68) 


den folgenden Bestwert an: 


2333 ns \4/ 1/3 42/3 
Cot = = Me 77 [2 (nam + A)]YS 2208. (69) 


Andererseits macht bei einer Turbine von gegebenem n,, A, 7), 75m und gegebenem Gleitwinkel ¢ 
eine DurchfluBziffer von der Grobe 


= P A \4/313/7 
Phy opt 3 la —TNsm) h & | = 


den hydraulischen Gesamtverlust h, zu einem Minimum!. 
Setzt man die aus (68) und (70) sich ergebenden Durchflufziffern einander gleich, dann bekommt 
man die folgende, fiir die Kavitationssicherheit und Strémungsverluste optimale spezifische Dreh- 


zahl: 


my (L—v) V2 [ 3 € 3/4 [2 (ism + A)] 7/8 


Beispielsweise erhalt man fiir eine Turbine mit y = 0,45 7, = 0,9, sm = 0,9, ¢ = 0,01, A = 0,12 
einen Wert von n,, = 795. 
Aus (69) und (71) ergibt sich mit n, = Bs oot ein beziiglich Strémungsverlust und Kavitations- 


sicherheit optimaler Kavitationsbeiwert von 
Cage =F Tae Mom +A) A, (72) 
opt 4A il oie Nsm sim 


Hieraus erhalt man fiir die Turbine mit den oben angenommenen Daten 0%; = 0,949. 
Aus der bekannten Beziehung 


n (73) 


? 


30 y mn (1 —v?)]4/2 
i = (2 g)?/4 ae wig? pri? 


1 iene die Arbeit von Keller und Bleuler in FuBnote 1 auf Seite 8. 
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welche exakt fiir eine iiber den DurchfluBquerschnitt konstante Meridiangeschwindigkeit gilt, er- 


5 5 5 . . . . . * O und 
gibt sich, wenn man n, = n,,,, und @ = Ps ops Setzt, die hinsichtlich Kavitationsfestigkeit 


Wirkungsgrad optimale bezogene Umfangsgeschwindigkeit 
13 6. (Ueto Gems aes (74) 
lid, = — Io - ——|. 
opt j2 E ae A | 


Hieraus erhalt man beispielsweise fiir cine Turbine mit 7;m = 0,9, € = 0,008, A = 0,15 den Wert 
Wid, = 1,76. Aus (74) kann bei gegebenem H und n der mit Riicksicht auf Kavitationssicherheit 
Tbe ateeisad optimale Durchmesser eines Laufrades bestimmt werden. ; 
b) Der Laufradentwurf in zweidimensionaler Behandlung mit Beriicksichti- 
gung der Vertraglichkeit der Strémung vom Laufrad und Leitrad. Im Folgenden 
wird dem Laufradentwurf das vom Leitapparat bei fortgedachtem Laufrad vor diesem erzeugte 
Geschwindigkeitsfeld mit den Komponenten c,,,, und Cur. zugrunde gelegt. Die Leitschaufel- 


dffnung wird dabei so gewahlt, daB fiir die aus der Entwurfsdrehzahl n, dem LaufradauBendurch- 
messer D,, dem Entwurfsgefille H und einem zunachst geschatzten hydraulischen Wirkungsgrad 
ny, folgende bezogene Umfangsgeschwindigkeit u;g an einem mittlerem Stromfaden annahernd 
drallfreier Austritt vorliegt. Mit dem Gitterverstarkungsfaktor K aus einer durch (3) symbolisierten 
Beziehung zu {3,9 und t/L, dem Auftriebsbeiwert c, gemaB (2) dem Winkel der ungestérten An- 
strémrichtung gegen den Umfang (Abb. 1) 


1 
Uid —= Cuz (Eng, — Cua) 
id MM id i 15 
Boo = are tg |" = _“__*], (75) 


dem Gleitwinkel im Gitterverband 


m k sin (Boo — Bpo)’ 


Gon (Po Bp 0) - (76) 


é = are tg 


(wobei der Zahler der rechten Seite in symbolischer Schreibart die Abhangigkeit des Widerstands- 
beiwertes c,, vom Anstellwinkel 6,, — {po verdeutlichen soll), der Bauersfeldschen Gleichung 


2t (Cu, ae: Cuz.) sin [Bre 


shad § cm,q (1 — tg ectg Ba) ’ 


Cc 


(77) 


dem hydraulischen Wirkungsgrad als Produkt der Wirkungsgrade von Laufrad- und Saugrohrstufe 


1 2 Fe 2 
hh siden nema (78) 
T 


sin (Bo — &) Sin Boo wid 


und der Turbinenhauptgleichung 


Wid (Cur. a Cus.) = “Sy (79) 


kann man bei aus der Zustrémung bekanntem Cm, und Cur. bei erfahrungsgema8 bekanntem 
u 
Cur Ssm und m und wahlweise vorgegebenen Werten Ca? Cua? t/L usf. die verbleibenden sieben Un- 
u 


bekannten aus den sieben Gleichungen (2), (3), (75), (76), (77), (78), (79) bestimmen. Zur Auslegung 


eines Fliigels bei gegebener Zustré6mung vom Leitapparat also gegebenem Cm,, und c,,  werde die 
i id 
Aufgabe auf die drei folgenden, praktisch bedeutsamen Faille I, 1 und III beschrankt: 
Fall I. Berechnung von £, Bp0>67 pas Cus. und t/L, wenn c, gegeben ist. 
Fall If. Berechnung von Bx, Bpo» €:1),, K, ¢, und t/L, wenn Cus. gegeben ist. 
u 
Fall III. Berechnung von Boor Bpo» €& Ny» K, c, und Cup, 9 Wenn i/L gegeben ist. 
ne Selbstverstandlich umfat das symbolisch geschriebene Argument auBer dem Anstellwinkel Bo Bpe 
die Gitterparameter, die Profilparameter und die Reynolsche Zahl, womit der umfangreiche Pragenkomplex 


der Bestimmung des Gleitwinkels im Gitterverband nur angedeutet ist. Naheres dariiber findet man in dem 
Aufsatz von H. Schlichting und N. Scholz, Ing.-Arch. 19 (1951) S. 42. 
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Um den Werkstoff auszuniitzen, wird man bei der innersten Teilturbine zweckmabig den durch 
die Ablésung bestimmten, gréBten Auftriebsbeiwert ¢, vorschreiben, also t/L gema®B Fall I be- 
stimmen. Bei der mittleren Teilturbine wird man zweckmaBig die der Leitradstellung zundchst 
annahernd unterlegte drallfreie Laufradabstrémung (cup = 0) an Hand des Falles II prazisieren. 

Mit der so bestimmten Fliigellange vom inneren und mittleren Fliigelschnitt ist wegen der er- 
winschten stetigen Fliigelumrandung das Teilungsverhaltnis t/L fiir die tbrigen Fligelschnitte 
vorgeschrieben. Bei diesen ist also nach obigem Fall III vorzugehen. Zur rationellen Lésung kann 
man sich dabei des folgenden Naherungsverfahrens aus dem Abschnitt 5 bedienen. Hiernach ergibt 
sich bei vorliegenden Werten Bpo» t/L, Osmo Cou» X, und einem zunichst geschatzten Gleitwinkel ¢ die 
bezogene Meridiangeschwindigkeit Cm,, von einer Teilturbine aus (46) und (47). Im vorliegenden 


Fall ist x, und ¢m,, vom Leitapparat her vorgegeben, das Teilungsverhaltnis t/L ist unbekannt. Zu 


seiner Bestimmung wird man an Hand von (46) und (47) fiir einige Werte von t/L die bezogene 
Meridiangeschwindigkeit Cm;, berechnen und iiber t/L aufzeichnen, um dann graphisch jenen Wert 


t/L zu erhalten, bei welchen Cm,, dem vom Leitapparat vorgegebenen Wert entspricht. 


Fiir den Profiltropfen wahlt man zweckmaBig die Profildicke s wie folgt in Abhangigkeit von 
der von der Fliigelnase aus zaihlenden Skelettkoordinate é 


P 1 a ll é \p 2 
jarred) [Esa] " 


Dabei folgt die maximale Fliigelstarke eines Fliigelschnittes s,,,, aus der Beanspruchung des 
Fliigels. Da die Austrittskante in Praxis keine Schneide haben kann, ist a = 1,03 — 1,1 zu wahlen; 
p hat man bei den itiblichen Profilen im Durchschnitt mit etwa 0,5 anzusetzen. Die mit (80) be- 
schriebene Profilform enthalt nicht die Eintrittskantenform. Diese ist entsprechend der FuBnote 
von Abschnitt 4 mit Riicksicht auf die Eintrittskantenkavitation auszubilden. Aus po, t/L und 
dem Profil kann man mit der Naherung 6,9 = 6 an Hand der im Abschnitt 4 beschriebenen Me- 
thode denjenigen Verlauf des Fliigelskeletts berechnen, der einen konstanten Druck auf der Saug- 
seite sichert. Fiir das somit festliegende Profil kann nach Abschnitt 3 von der durch fo gekenn- 
zeichneten auftriebslosen Richtung auf die Richtung der Skelettsehne gegeniiber der Gitterachse 
geschlossen werden. 


c) Die Vertraglichkeitsbedingung fiir die Strémung von Laufrad und Saug- 
rohr. Es werde angenommen, da fiir die Flissigkeitsteilchen im Unterwasserbehilter die spezi- 


fische Strémungsenergie H = = +h+ ee konstant sei. Setzt man diese willkiirlich gleich Null 


é . . 
und iibertragt man den Energiesatz fiir die wirkliche Rohrstrémung auf einen Stromfaden im 
Saugrohr, dann entspricht die spezifische Strémungsenergie am Laufradaustrittsquerschnitt H, 


2 2 F 
der zu (sm st - bua angenommenen Verlusthéhe der Strémung innerhalb des Saugrohres. 
& 


Nimmt man weiter an, dafs H, annahernd der Eulerschen Gleichung in der Form 
ibe 
grad H, = F [Co rotsc,| 


geniige, und daf der Gradient von H, keine axiale! und tangentiale Komponente besitze, dann 
laBt sich zeigen, daB auch seine radiale Komponente verschwinden muB, daB also tiber den Lauf- 
radaustrittsquerschnitt der Zusammenhang 


Com get ap (Be ae -_= konst. (81) 


gilt. . . Se — . 
Die Vertraglichkeitsbedingung (81) schreibt bei gegebenen Werten von ¢,,, und Csm, ee ba 
obigen Annahmen aus Cmy., = mig bekannten Meridiangeschwindigkeit und bei der willkirli 


gewahlten Umfangskomponente c¢,, eines Fliigelschnittes diesen Wert fiir alle itbrigen Fliigel- 
schnitte vor, so da man hierbei nach dem in Abschnitt 7b) beschriebenen Fall II zu verfahren 


1 Ein Verschwinden der axialen Komponente von grad H kann natiirlich nur annahernd vorliegen, da Hy 
langs des Saugrobres um Csm (¢m3/2 g) + Csu (Cu3/2 g) abnimmt. 
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hatte. Dabei ist zu bedenken, dafB bei gegebenem c,, s und uja wegen der Hauptgleichung (79) mit 
Cus, auch der hydraulische Wirkungsgrad 7, bestimmt ist. Liegt nun 7, bei einzelnen Fliigel- 


schnitten auf Grund der Berechnung auferhalb verniinftiger Grenzen, dann sind Stellung oder 
Form des Leitapparates so zu dndern, dal} Cur. und ¢m,, in Verbindung mit (81) und (79) zu trag- 


baren Werten von 7, fiihren. 


8. Niherungsmethoden zur Beriicksichtigung der Einfliisse bei dreidimensionaler Laufrad- 
strémung. a) Modell fiir die gebundenen Wirbel der Restfliigel’ bei ruhendem 
Laufrad. Es werde angenommen, da8 sich die Zirkulation um einen Laufradfligel in Richtung 
des Halbmessers nicht andere und daB® die Induktionen der von den Fliigelenden ausgehenden 
freien Wirbel bei den iiblichen kleinen Laufradspaltweiten infolge der spiegelbildlich zu den Rand- 
stromflaichen zu denkenden, entgegengesetzt drehenden Wirbelfaden zumindesten fiir die maf- 
gebenden, vom Spalt entfernten Aufpunkte verschwinden, die freien Wirbel also vernachlassigt 
werden kénnen. Setzt man weiter eine Zustrémung von konstantem Drall langs des Halbmessers 
und drallfreie Abstrémung voraus, dann kann man einen beliebigen gebundenen Wirbelfaden, 
wegen dessen Quellenlosigkeit, vor dem Laufrad langs der Turbinenachse fortsetzen, wobei diese 
Fortsetzung aller gebundenen Wirbel den vorausgesetzten Eintrittsdrall erzeugt. Innerhalb des 
Laufrades schlieSt sich an diesen axialen Teil eines gebundenen Wirbelfadens ein die Naben- und 
AuBenmantelflache normal durchdringender bis ins Unendliche reichender Teil an, der von dort 
aus wegen seiner Quellenlosigkeit wieder im Unendlichen vor dem Laufrad auf die Turbinenachse 
zu verlaufen muB. Der ins Endliche reichende Teil dieses Modells von einem gebundenen Wirbel- 
faden besteht also aus zwei unendlich langen geradlinigen Schenkeln, von denen einer vom Unend- | 
lichen vor dem Laufrad kommend bis zum Laufrad mit der Turbinenachse zusammenfallt und der | 
andere hierzu um den Winkel (7/2) — v geneigt ist. 

Die Induktion des axialen Wirbelteils, welche dem konstant angenommenen LEintrittsdrall 
gleichkommt, wird beim Entwurf des fiir sich betrachteten Fliigelschnittes durch dessen mit Cy, t 
vorgegebene Zirkulation beriicksichtigt. Will man die raumliche Anordnung der gebundenen 
Wirbel in erster Naherung beriicksichtigen, dann hat man nach dem Vorgang von Betz? das Skelett 
dieses fiir sich betrachteten Fliigelschnittes infolge der Induktionen der Restfliigel zu verzerren, 
damit die aus den Geschwindigkeitsdreiecken folgende Zirkulation fir einen Fliigelschnitt erhalten 


y 


Abb. 9, Zur Induktion eines vom Punkt 0 ausgehenden geraden Wirbelfadens der Zirkulation Pr, 
welcher um den Winkel » gegen eine zur Turbinenachse senkrechte z-y-Ebene geneigt ist. 


bleibt. Dabei werde in erster Naherung ein heliebiger Restfliigel durch einen um 77/2 — y gegen die 
Turbinenachse geneigten Stabwirbel ersetzt, welcher in einer Axialebene verlaufe, die den Schwer- 
punkten der Zirkulation yon den Zylinderschnitten des Fliigels méglichst benachbart ist. Es 


1 Unter ,,Restfliigel,, seien alle Fliicel verstanden mit jeni 
; i é ig a it Ausnahme desjenigen, der den Aufpunkt a 
2 Siehe die Arbeit von Betz in FuBnote 5 auf Seite 1. aie 2 eres 
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Interessieren daher nur die Induktionen eines 


Wirbelfadens. 


derartigen zur Turbinenachse geneigten geraden 


Auf Grund des Biot-Savartschen Gesetzes induziert der in Abb. 9 dargestellte von 0 ins Unend- 


- gehende, gerade Wirbelfaden von der Zirkulation J” » dessen Richtung durch den Einheits- 
vektor 


Boat ith eae )° cos » cos Op + 4° cosy sin Oy, (82) 


bestimmt ist (wobei ¢°, )°, 3° die Einheitsvektoren in x, y und z-Richtung darstellen) im Punkt P 
(x, y =rcos 3, z=rsin ?) die Geschwindigkeit 


I 1l+cosy 


a ee ee (83) 


40a siny 
Dabei ist p° der Einheitsvektor der Strecke O—P 


x o reos? o rsin? 


p= jere ry jetr je@rr’ (84) 
a der Abstand des Punktes P von dem geraden Wirbelfaden 
a = psin y = |[p, 3°]| (85) 
und y der von den Einheitsvektoren p° und 8° eingeschlossene Winkel 
y = are cos (p°, 8°). (86) 


Stellt man das Vektorprodukt auf der rechten Seite von (83) in Komponenten dar und setzt 
aus v, und v, nach bekannten Transformationstegeln die radiale Komponente v, und die Um- 
fangskomponente v, zusammen, dann bekommt man mit den Abkiirzungen A) = 3, — 0 und 
_ €=4x/r die folgende axiale, tangentiale und radiale Komponente von ): 


Ce " cos v sin Ad Esinvy + cos AP cos v 87 
eta MOOS Gin vcs APF con) jite , (87) 

anges E cos » cos Ad — sin v Lot | Esinv + cosA’ cosy (88 
°u ~ far sin? 48 + (sin v cos AD — E cos v2 i jl+2 y ) 

ae HF Ecosyv sin AP " & siny + cos Ad cos v (89) 
ar ce avon 1 cosy irae.) 


Die raumliche Anordnung der gebundenen Wirbel kann mau beim F'liigelentwurf so beriicksichtigen, 
da man fir die vorgeschriebene Zirkulation J’ = Ac, t und eine vorgeschriebene Druckverteilung 
auf der Saugseite zunadchst ein Einzelprofil entwirft, dessen Eigenschaften im raumlichen Gitter- 
verband erhalten bleiben, wenn man nach dem Vorgang von Betz seine Skelettelemente in Ab- 
 hangigkeit von den Induktionen der ,,Restfliigel“* am Ort des betrachteten Skeletts verzerrt. Dabei 
geniigt es in erster Annaherung, wenn man die Skelettelemente um den Winkel Af. dreht, um 
-_-welchen sich das Element einer Stromlinie am Ort des betrachteten Fliigels bei Uberlagerung der 
Induktionen der Restfliigel dreht. 


Sind Yup und Vxp als Induktionen der Restfliigel in tangentialer und axialer Richtung klein 


gegen W.., dann ist 


ABo = — eR gin Ba oo C08 Boo (90) 


Bbc 


die Drehung eines Skelettelements im Zylinderschnitt. Dabei sind die Induktionen des Restgitters 
an einem Aufpunkt des betrachteten Fliigels gegeben durch 


if if Peel: poked 
Of 4ar Ev, B, é, 0) ° Vur a Amr FQ, od b, d) ’ i awe Anr EQ, e é, 3) ; oD) 


m R 


2 
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Hierin folgen die EinfluBfunktionen F' aus (87), (88) und (89), wenn man den Langenwinkel des 
n-ten Nachbarwirbels #p mit 2 7/z n ansetzt und fiir jede Komponente der Induktionen die Anteile 
der z — 1 Fliigel des Restgitters in nachstehender Weise summiert zut 


ener, 
a—1 cos sin (22 n— ® 
z 


Fy, 2,€, 0) = SS = = ——— —___— ———- = 
2 sat? n—6) + [inv cos “Fn — 8) —£ ones 


gsiny + cosy cos (72"— 0) 
j1+2 


: (92) ) 


2mn : 
z—l Ecosy cos|—,  —#]— sin v 


POehS ie 2 
2 sint (“F"— 8 -- sin v COs (-5*—0|—f os] 


ae Dae) 
— sin vy -+ cos v cos —¢ 


Zz i 
el Ee ee 3 93 
yite O23) 
z—1 s cosy sin (7 "— a) 
F(y, z, Ss 0) = - - = — x aE 
2 sit(“F"— 8) + |siny cos (“2 — 9) — cosy] 


Sse 2an 
Sin + 608» | —?o 
z 


1 aby ae 94) || 
; ! yi+2 Oe 


Der Verlauf der EinfluBfunktionen F’, = konst., F', = konst. und F’, = konst. ist fir z = 4 
und z = 6 Wirbelstrahlen in einer Achsnormalebene und fiir z = 4 Wirbelstrahlen, welche um 20° 
gegen die Achsnormalebene geneigt sind, in den Abb. 10 bis 18 in einem ?, -Koordinatensystem 
dargestellt. 

Auf diese Weise kann man fiir einen zunichst als Einzelskelett von bestimmter Zirkulation 
entworfenen Zylinderschnitt vom Radius r = 1, die gemaB (90) fiir die Drehung eines Skelettele- 
ments erforderlichen Induktionen der raumlich angeordneten Restfligel sofort abgreifen. 

b) Modell fiir die Wirkung der Rotation des Laufrades bei raumlicher Fligel- 
anordnung. Die absolute Strémung durch das Laufrad werde wirbelfrei und von konstanter 
Meridiangeschwindigkeit angenommen. Die betrachtete Relativstrémung in einer mit Meridian- 
geschwindigkeit durch das Laufrad bewegten Achsnormalebene sei stationaér. Die Strémung sei | 
aufen und innen durch koachsiale Zylinder begrenzt und das Fliigelskelett habe die Form einer 
geraden Schraubenflache von der Gleichung 


x=Origp, (95) | 


wobei § der Steigungswinkel eines Zylinderschnittes vom Halbmesser r, x die axiale, § die azimutale 
Koordinate sei und 


rtg 6 = konst. (96) 
fiir alle Zylinderschnitte gelte. 


Die radialen Schaufelspuren, sowie Nabe und Mantel begrenzen damit in einer Ebene normal 
zur ‘Turbinenachse einen ebenen Kreisringausschnitt. Bewegt sich dieser mit c,, = konst. in axialer | 
Richtung, dann besteht er stets aus denselben Fliissigkeitsteilchen, welche relativ zum Laufrad | 
mit einer in Drehrichtung positiven Winkelgeschwindigkeit | 


or =— “otgB =—(1—4) “ =— (1a (97) 


_ + Die Auswertung der Gleichungen (92) bis (94) wurde im Rahmen einer Diplomarbeit an der Fakultat 
fir Maschinenwesen und Elektrotechnik der T.H. Miinchen von Herrn Hansen besorgt. 
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rotieren, wobei angenommen wird, daf ¢,, ctg fp & u baw. € A 0 ist, die Schaufel also keinen senk- 
rechten Austritt erzwingt. Da die Flissigkeit voraussetzungsgemals relativ zum Laufrad mit der 
Winkelgeschwindigkeit @ rotiert, haben die Teilchen relativ zu der mit c,, bewegten Achsnormal- 
ebene in Drehrichtung eine Winkelgeschwindigkeit von 


We =—O—OoOr=—loO. (98) 

Damit lautet die Grundgleichung fiir die Geschwindigkeit tv der Strémung in einem mit Meridian- 
geschwindigkeit durch das Laufrad bewegten Kreisringausschnitt zwischen zwei Fliigeln 

rot ty =—2Co. (99) 

Driickt man den Rotor in bekannter Weise durch die radiale und tangentiale Geschwindigkeits- 

komponente w, und w, aus, und leitet man diese in der folgenden iiblichen Weise von einer Strom- 


funktion p abt 


fe (100) 

or ? 

ow 
cot fang ast 101 
oe roo” uy) 

dann bekommt man als Grundgleichung fiir die untersuchte Strémung: 

yy Lvov earl eC ea. 102 
or? r or | r 07 ac ee 


Lat man y auf dem Rand des Kreisringausschnittes verschwinden, welcher bei z Fligeln durch 
die Geraden # = (z/z) und die Kreise mit r =17;, r = 1a gebildet wird, dann lautet die Lésung 
von (102) mit den Abkiirzungen 9 = (r/r;), 0, = (T4/i) 


oe. : In nz td : 
p=—CZeri+Gme+ >a, sin [nx ino) © fee) Ge (103) 
n=l 
wobei 
oa 
C,=Cor ea ; (104) 
Baye 
(eee (105) 
und 
F a5 ; I a 
c COAT, [1 = ( 1) ] [1 r Gal fe 4 In? Qa + n? 7? 
OF =a sae oS Ss é (106) 
2 Co} poe 
z In Qa 


Abb. 19 zeigt den Stromlinienverlauf innerhalb des Kreisringausschnittes. Man erkennt dar- 
aus, daB bei € > 0 die Flissigkeit auf der Druckseite des Fliigels nach auBen strémt und auf der 
Saugseite nach innen. Vom Staupunkt der Eintrittskante aus entfernen sich druck- und saug- 
seitige Stromlinie in entgegengesetzter radialer Richtung. 

Liegt die beschriebene Sekundarstrémung innerhalb einer Achsnormalebene auch am Laufrad- 
austritt vor, dann ist damit ein Energieverlust verbunden, da die radiale und tangentiale Geschwin- 
digkeitskomponente in einem dem Laufrad im allgemeinen folgenden Saugrohr nicht in Druck 
umgesetzt werden kénnen. 

Da die untersuchte Ebene sich laut (97) mit einer Winkelgeschwindigkeit von wy relativ zum 
Fligel bewegt, setzen sich die absoluten Geschwindigkeitskomponenten c, und c, am Laufradaustritt 
wie folgt zusammen: 


¢,=u,+(or+o)r=w,+lor, (107) 
Cry (108) 
Der auf H bezogene Verlust infolge der oben beschriebenen Sekundarstrémung wird damit 
1 
h= aE [| [(wy + Cm 1)? + w?] dF (109) 
(F) 


! H. Lamb, Lehrbuch der Hydromechanik, Leipzig 1931. 
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oder, in der Stromfunktion ausgedriickt, 


ess SAL i Oeey ae. aw aw \2 
h 2ga(r2 —r?) H ‘hy les == Dae Or aS + Ga te C2 r2 os rdr dé. (110) 
(F) 


Saugseite 


Abb. 19. Relativstroémung in einer mit c, = constans durch das Laufrad bewegten, zur 


m 
Turbinenachse senkrechten Ebene bei z = 4 Fliigeln € > 1 und einem Nabenyerhaltnis von 0,45. 


Die partielle Integration des ersten und zweiten Gliedes im Integranden nach r und die des dritten 
Gliedes nach # ergibt im Hinblick auf den verschwindenden Randwert von y und mit Riicksicht 
auf (102) den folgenden Ausdruck fiir den Energieverlust h infolge Drehung des Laufrades: 


C2 wr? (02 —1 Bz, 0 | 7 =~ ; if 1 mn le ma \) 
= 4g H “In Ga oa (Qa t 1) In Qa »,, m? 7? + 4 In? 0g |m a 4 In? 0g + m? x? 2G 2 In ga}| © 
n=0 
(111) 


Dabei werden mit m = 2n + 1 die gemaB (106) allein existierenden ungeraden Glieder von (103) 
beriicksichtigt. Die Fliigelzahl z mu selbstverstandlich von Null verschieden sein. 

Im Gegensatz zu dem annahernd mit ¢ (Wia/em, ,) anzusetzenden Laufradverlust ist der durch die 
Drehung des Laufrades entstehende Verlust vom Quadrat der bezogenen Umfangsgeschwindigkeit 
abhangig. 

aw? r2 


2gH 


2 
= Uid - 
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Eine einfache geometrische Naherungsbetrachtung zeigt, daB der Beiwert € mit der Abweichung 
vom drallfreien Schaufelaustrittswinkel Af wie folgt zusammenhangt: 


Mey akon 0) (112) 


9. Versuche an einem in Anlehnung an Abschnitt 7 entworfenen Laufrad. Ein im Auftrag der 
Firma E. Ruch (Oberkirch) zu entwerfendes 4fliigeliges Laufrad fiir eine offene Kaplanturbine, deren 
Meridianschnitt durch Leitapparat und Laufrad in Abb. 20 zu sehen ist, bot Gelegenheit, die 
aufgestellten Berechnungsgrundsatze an einem 
Modellrad von D, = 350 mm LaufradauBen- 
durchmesser im Versuch zu tiberpriifen. 

Zunachst wurden fiir die im Grundrif dar- 
gestellte Leitschaufelstellung bei radialer Zu- 
stromung zum Leitrad, wie sie bei einer offenen 
Turbine in gro8em Oberwasserbehalter an- 
genommen werden darf, nach dem Vorgang 
von Strscheletzky! der in Abb. 20 dargestellte 
Verlauf der Meridianstromlinien bei einem 
Gesamtstrom von Q = 0,32 m?/sec berechnet. 
Der sich gleichzeitig hiermit ergebende Ver- 
Jauf der Komponenten der Eintrittsgeschwin- 
digkeit Ors und c,,., entlang der Normalen zu 


Zustromung 
radial — 


“lid 
den Meridianstromlinien n— n ist in Tabelle 3 
fir die drei betrachteten Stromlinien an- 
gegeben. Eine durch die in die Strémung 
hineinragende Leitschaufelkante verursachte 
steile Zacke in der Geschwindigkeitsverteilung 
wurde dabei unterdriickt, weil sie im Wider: 
spruch mit der Erfahrung gestanden hatte. 
Fiir das Laufrad wurde eine bezogene Umfangs. 
geschwindigkeit am Aufendurchmesser D, 
von uj, = 1,71 zugrunde gelegt, welche mit 


i dem Gesamtstrom im Entwurfspunkt eine 
nop gu, Moran durch Leteppert ond Laut eZetule speaifische Drehzahl von etwa n, = 600 U/imi 
Zustrémung zum Leitapparat und die abgebildete Leitschaufelstellung ergab. 
eae An dem der Nabe benachbarten Zylinder 
schnitt t vom Durchmesser 0,48-D, wurde 
das Teilungsverhaltnis t/L so bestimmt, da® bei einem geringen negativen Austrittsdrall der Auf 
triebsbeiwert nicht zu Ablésungen fihrte. 

Im mittleren Zylinderschnitt m vom Durchmesser 0,71 - D, sollte bei der entsprechend gewahlter 
Leitradéffnung und dem zu erwartenden hydraulischen Wirkungsgrad drallfreier Austritt vorliegen 

Am auBeren Zylinderschnitt a* vom Durchmesser 0,94 D, wurde ein geringer negativer Aus 
trittsdrall vorgesehen, wobei im Teilungsverhaltnis auf eine kontinuierliche Umrandung des Fliigel 
zu achten war in Verbindung mit einem festigkeitsmaBig erwiinschten Zusammenfall von Profil 
schwerpunkten und Fliigelachse. Der Austrittsdrall in den extremen Zylinderschnitten sollte di 
Strémung im Saugrohr stabilisieren und durch Herabsetzen von w,, in den auGeren Fligelpartie: 
deren Kavitationssicherheit erhéhen. AuSerdem konnte hierdurch bei Abweichung vom Kon 
struktionsgefalle nach oben oder unten fiir die innere oder auBere Teilturbine ein gréRerer Wirkungs 
grad als iiblich erwartet werden und damit ein flacher Verlauf des Gesamtwirkungsgrades iiber de 
Wallhohe. 

Die maximale Fliigeldicke an der Nabe ergab sich aus der Beanspruchung des Fligelwurze! 
querschnittes bei Hy,q,. Von hier aus zum AuBenrand wurde $mqx linear iiber dem Halbmesser au 
das technologisch vertretbare Ma vermindert. Die Profiltropfen wurden affin zum Profil Gét 
tingen 398 ausgebildet. Der Neigungswinkel f,, der auftriebslosen Anstrémrichtung gegen de 
Umfang wurde nach Abschnitt 7b) so festgelegt, daB die in Tabelle 3 angegebene Geschwindig 
keitsverteilung vor dem Laufrad mit der vom Leitrad her berechneten zusammenstimmte, Hierbs 
wurde der Saugrohrwirkungsgrad fiir die Umsetzung der Umfangskomponente 7,,, zu Null angenow 


1 Siehe die Arbeit von Surscheletzky in FuGBnote 6 auf Seite 1. 
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men. Der Saugrohrwirkungsgrad zur Umsetzung der Meridiangeschwindigkeit 7, wurde zu 0,8 
angesetzt, ein Wert, welcher annahernd durch nachtragliche Messung bestatigt wurde und der 
teilweise durch wasserbauliche Wiinsche bedingt war. Die Steigung m der c,(6))-Kennlinie wurde 


Abb, 21. Zylinderschnitte a, m, i, des Fligels. 


aus den auf unendliche Fliigelbreite umgerechneten Polaren annahernd ahnlicher Profile zu m—=4,8 
im Mittel angesetzt, wobei fiir die Profile i und m die der Gottinger Reihen 496 und 491 genommen 
wurden. Fiir die Werte B50, Smax, t/L wurde mittels der 

im Abschnitt 4 behandelten Methode der Verlauf der Tabelle 3. Daten der Zylinderschnitte. 


Skelettlinie fiir einen konstanten Druck auf der Saugseite Zylinder- : ¥ : 
festgelegt. Hieraus und aus dem Profiltropfen wurde schnitt 
nach der Methode von Pantell1, erganzt durch das im 
Abschnitt 3 angegebene Verfahren von der durch fo uid 0,8225 1,215 L0L0 
gekennzeichneten auftriebslosen Richtung auf den Staf- "h 0,88 0:87 Bee 
felungswinkel / geschlossen. Cutid 0,39 034 0,34 
Die Entwurfsdaten der so erhaltenen und in Abb. 21 Cmig 0,61 | 0,04 | Cont 
dargestellten Zylinderschnitte sind in Tabelle 3 zu- L/t Aol 1,07 0,935 | pe} 
sammengefaft. Dabei ist der hydraulische Wirkungsgrad Bpo@) re | age ne 
Ca | 3 | > ot 


im Mittel etwa 1°% gréBer als der gemessene Wert an- 
gesetzt, um das Verhalten der geplanten GroBausfiihrung 
des Rades zu kennzeichnen. Der fiir die Berechnung vom Auftriebsbeiwert verwendete Gleitwinkel 
wurde dem in der Tabelle 3 angegebenen hydraulischen Wirkungsgrad 7, und dem gemessenen 
Saugrohrwirkungsgrad 7,,, angepaBbt. 


Th 
—— — Entwurfspunkt 
a) Sol 
085 5: —17118° 
13° Hide 
Z ir r 
La 
we SS 
98 B° ike ae 
— 
i 
| 334 
G75 Pu =3° i ara 
i al 
at =] fs call Le Bi fs \y 4 
= y 


Abb. 22. Hydraulischer Wirkungsgrad 7, itiber dem Einheitsstrom Q,, bei konstanter Fliigelstellung 6, im Entwurfspunkt, 
Einheitsdrehzahl n,, = 145. 


Die Versuche des Modellrades von 350 mm Durchmesser wurden an einem Stand der Technischen 
Hochschule Miinchen durchgefiithrt, der durch 3,5 m Fallhéhe und etwa 0,90 m?/sec Vollwasserstrom 


3 Snes Arbeit von Pantell in FuBnote 1 auf Seite 2. 
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gekennzeichnet ist. Dabei wurde die Turbine in einem Laufradwinkelbereich von etwa 35°, in 
cinem Leitradwinkelbereich von etwa 60° und in einem Drehzablbereich von 500—1000 U/min 
abgebremst. Die Fallhéhe wurde durch pneumatische Ubertragung mit einem Differenzdruck- 
manometer, der Gesamtstrom durch die Volumendifferenz eines geeichten Behalters innerhalb 
einer gemessenen Zeit, die Drehzahl mittels Stichdrehzahlmesser von der Firma Hasler und das 
Drehmoment an der Kupplung mittels Pendeldynamo gemessen. Das zugehorige Laufradmoment 
wurde im Schleppversuch mit geeichtem Motor bestimmt. 

Abb. 22 zeigt den Verlauf des hydraulischen Wirkungsgrades tiber dem Einheitsstrom Q,, bei 
der Konstruktionseinheitsdrehzahl n,, fiir verschiedene Laufradstellungen, welche durch den 
Winkel gekennzeichnet sind, den die Skelettsehne des auBeren Zylinderschnittes mit der Umfangs- 
richtung einschlieBt. Man erkennt, daB der Auslegungspunkt mit einem gemessenen Wert? von 
Q,, = 1,37 etwas von der Einhiillenden der Wirkungsgradlinien und deren Maximum abweicht. 


4 on; 
40 —- 


0 10 60 


Abb. 23. Vergleich der gemessenen und berechneten Cm (u;4) Werte. 


Um das in Abschnitt 5 beschriebene Verfahren zur Berechnung der ¢n,, (uja)-Kennlinie einer 


Teilturbine beurteilen zu kénnen, wurde fiir den vom Leitapparat her berechneten Winkel «, der 
mittleren Stromflache die in Abb. 23 dargestellte ¢,,,(uia)-Linie berechnet. Dieser Linie sind in 


Abb. 23 die gemessenen Werte des mittleren Stromfadens gegeniibergestellt. 


10. Zusammenfassung. Es wird angegeben, wie eine aus Leitrad, Ubergangsraum, Laufrad 
und Saugrohr bestehende Kaplanturbine mittels 2-dimensionaler Methoden ausgelegt werden kann, 
wobei die Vertraglichkeit der Strémungsfelder von zwei aufeinanderfolgenden Bauteilen und die 
kavitationssicherste Druckverteilung auf der Saugseite des Laufradfliigels besonders beriicksich- 
tigt wird. 

Zur Berechnung des Verlaufes vom Gesamtstrom Q iiber der Drehzahl n wird auf der Grundlage 
2-dimensionalen Potentialstrémung durch ein Profilgitter fiir die Teilturbine eine tabellierte, ge- 
schlossene Naherungslésung der c,,(u)-Linie angegeben. 

Weiter werden geschlossene Naherungslésungen mitgeteilt fiir die Form eines Fliigelskeletts im 
Gitterverband von diinnen Profilen bei konstantem saugseitigen Druck, fiir die Geschwindigkeits- 
verteilung hinter dem Leitrad und hinter dem Laufrad bei wirbelbehafteter Strémung, fiir das 
3-dimensionale Geschwindigkeitsfeld der Restfligel bei raumlicher Fligelanordnung und fiir den 
Stromlinienverlauf sowie den Verlust infolge Drehung des Laufrades. 

In Anlehnung an die angefiihrte Berechnungsmethode wird eine Modellturbine von 350 mm 
Laufraddurchmesser und etwa 13 PS mittlerer Leistung gebaut und auf dem Priifstand des Hydrau- 
lischen Instituts der TH Miinchen durchgemessen. Das Wirkungsgradoptimum des Rades von 
n, © 0,87 weicht geringfiigig von den Entwurfsdaten ab. Der gemessene Verlauf des Gesamt- 


stromes Q iiber der Drehzahl n weicht etwa 5% von der fiir die mittlere Teilturbine berechneten 
Kennlinie ab. 


(Eingegangen am 24, Marz 1958.) 


Anschrift des Verfassers: Dr.-Ing. J. Raabe, Privatdozent der Technischen Hochschule Munchen, Torkenweiler 
bei Ravensburg, Tettnanger Str. 43. 


' Die Turbine war fiir einen Einheitsgesamtstrom von Q,, = 1,43 ausgelegt worden. 
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Die Kigenfrequenzen und Schwingungsformen von Turbomaschinen 
Von B. Jager 


1. Zweck der Untersuchung. Die vorliegende Untersuchung ist fiir die Berechnung von Turbo- 
Triebwerken fiir Fahrzeuge und Flugzeuge gedacht. Hier interessieren die kritischen Drehzahlen, 
welche durch die Unwucht des Rotors erregt werden. Die Schwingungsfrequenz entspricht also 
der ersten Ordnung der Drehzahl. 

Die schwingungsmaBige Analyse eines solchen Triebwerks ist in den meisten Fallen etwa die 
folgende: Im untersten Drehzahlbereich sind Aufhangeschwingungen des elastisch aufgehangten 
Triebwerks zu erwarten. Dabei kann die Maschine selbst als starr betrachtet werden. Im mittleren 
Drehzahlbereich ist das Triebwerk iiberkritisch beziiglich der Aufhangung. Hier treten Schwin- 
gungen des ganzen Triebwerks an nicht rotierenden Teilen auf, etwa Biegeschwingungen von auf 
dem Gehause befestigten Teilen. Dagegen diirfen Schwingungen, an denen rotierende Teile mit 
wesentlichen Verformungen teilnehmen, im Drehzahlbereich der Maschine nicht vorkommen. Bei 
der geringen Dampfung wiirden solche Schwingungen, welche der Biegekritischen einer Welle 
auf festen Lagern vergleichbar sind, erfahrungsgema®B unzulassig sein. 

Da méglichst schon bei der Konzeption einer Maschine eine solche Analyse wenigstens ange- 
nahert vorliegen sollte, wird hier gezeigt, wie ein Ersatzsystem aufgestellt und wie das dazu ge- 
hérende Gleichungssystem behandelt werden kann, um das Verhalten eines Triebwerkes voraus- 
zusagen. 

Im folgenden ist versucht worden, eine Grundlage fiir die praktische Rechnung zu bieten. 
Deshalb ist ein méglichst allgemeiner Fall gewahlt worden, an dem die Operationen gezeigt werden. 


2. Die Einflu®zahlen der elastisch gelagerten Welle. Im Schrifttum iiber Biegeschwingungen 
ist eine Anzahl von Rechenverfahren beschrieben, welche es erlauben, unter den verschiedensten 
Gegebenheiten die unterste oder auch noch héhere Eigenschwingungszahlen von Wellen zu be- 
rechnen. Diesen Verfahren ist gemeinsam, da die Welle sich auf steifen Lagern befindet, oder daB 
eine gewisse Lagerelastizitat noch mit beriicksichtigt werden kann. Die wesentliche Elastizitat 
des schwingungsfahigen Systems liegt aber in der Biegung der Welle. 

Im modernen Maschinenbau, vor allem im Fahrzeug- und Flugzeugbau mit den hohen Dreh- 
zahlen und dem immer kleiner werdenden Leistungsgewicht kommt man darauf, die Wellen durch 
Hohlwellen ja sogar durch Trommellaufer zu ersetzen, bei denen die Beriicksichtigung der Biegung 
auBer in den Zapfen und Endscheiben keine wesentlichen Beitrage mehr liefert. 

Dagegen ist man bestrebt, die Gehause, in denen die rotierenden Teile gelagert sind, immer 
leichter auszufiihren, einerseits aus Griinden der Gewichtsersparnis, andererseits aus strémungs- 
technischen Griinden. Zum Beispiel will man im Turbomaschinenbau die ringférmigen Durch- 
fluBquerschnitte nicht mehr als unbedingt erforderlich durch Arme einschrénken. Denkt man 
sich, etwa im Rahmen eines statischen Versuchs, die Gehause in geeigneter Weise in ihrer auBeren 
Kontur festgehalten und eine Kraft auf dem Laufer angebracht etwa an der Stelle, an der die 
groBte Masse sitzt, so kann man in einem solchen Fall feststellen, dal die Weichheit der Gehause 
einen Beitrag zur Verschiebung des Angriffspunktes der Kraft liefert, der durchaus in der Grofen- 
ordnung der Verformung des Laufers selbst liegt. 

Will man nun fiir diesen Fall, also unter der Annahme der festgehaltenen AuBenkontur der 
Gehause oder unter der Annahme einer unendlich groBen Masse in den AuBenringen der Gehause, 
das Schwingungsverhalten des Laufers betrachten, so wird man 
bei der Aufstellung der EinfluBzahlen des schwingenden Systems 
die Verschiebungen, die auf die Elastizitat der Gehause zuriick- 
gehen, den Verschiebungen hinzufiigen miissen, die durch die 
Biegeweichheit des Laufers zustande kommen. Es liegt der 
Rechnung die Annahme zugrunde, da® alle Elastizitaten CINE) (inet Biegeame Welle aul clastischen Dagern: 
lineare Charakteristik besitzen. Man kommt auf das Schema fiir 
den Laufer von Abb. 1. Dabei gilt folgendes: Punkt 1 und 2 sind fest; Punkt 3 und 4 stellen die 
Masse der Lager mit einem Teil der Masse der Welle dar; Punkt 5 und 6 stellen die Massen der 
rotierenden Teile mit Teilen der Wellenmasse dar. Die Federn agg und 44 setzen sich aus allen 


3 
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Weichheiten zwischen den Festpunkten J und 2 und der Welle zusammen. Man sieht, daB sich die 
Anzahl der Massen des Systems gegeniiber der klassischen Art des Rechnens vergroBert hat, da 
auch die an den Lagern befindlichen Massen an der Bewegung teilnehmen und in die Rechnung 
eingefiihrt werden miissen. 

Die Bestimmung der EinfluBzahlen zerfallt nun in zwei Rechengange: In die Ermittlung der 
BiegeeinfluBzahlen des Laufers unter der Annahme steifer Lager und die Ermittlung der Lager- 
einfluBzahlen unter der Annahme, dai der Laufer starr sei. 

a) Die Ermittlung der Biege-EinfluBzahlen des Laufers unter der Annahme 
starrer Lager. In den Lagerpunkten soll keine Verschiebung quer zur Rotationsachse des Lau- 
fers méglich sein; dagegen sei eine Neigungsanderung der Mittellinie des Laufers zugelassen, da 
in den Lagern keine Einspannung bestehen soll. Fiir die Berechnung der Einflu®zahlen sei an 
sich auf das Schrifttum verwiesen; es soll aber des Zusammenhangs wegen ein kurzer Uberblick 
gegeben werden. 

Auf einer Welle sollen sich einige scheibenférmige Massen mit der Masse m; und dem aqua- 
torialen Massentragheitsmoment 9; befinden. Dann greift in der biegekritischen Drehzahl, wenn 
die Welle also mit der ersten Ordnung der Drehzahl schwingt, in einem Massenpunkt die Flieh- 
kraft 


te Ree (1) 


an. Dabei ist y; die elastische Auslenkung und @ die Winkelgeschwindigkeit der Rotation. Da 
die Massen scheibenférmig sein sollen, greift in jedem Massenpunkt auch das riickdrehende Kreisel- 
moment 


M; = — 0; 0 9; (2) 


an, wobei g; die Neigung im Punkt 7 ist?. 

Die Einflu8zahlen sind dem Schrifttum? zu entnehmen; es soll hier nur noch folgendes zusammen- 
gefaBt werden. 

Die EinfluBzahlen «;, geben die Verschiebung im Punkti unter der Einheitskraft im Punkt k 
an. Dabei gilt der Satz der Reziprokitat o;, =a; Die EinfluBzahlen 6; bezeichnen die Nei- 
gung der Biegelinie an der Stelle i unter dem Einheitsmoment an der Stelle k, wobei ebenfalls 
gilt Bi, = Pri Die EinfluBzahlen y;, sind als Neigung im Punkti unter der Einheitskraft im 
Punkt k definiert. Gleichzeitig gilt auch, daB y;, die Durchbiegung im Punkt k unter dem Einheits- 
moment im Punktz ist. Bei der zweiten Definition scheinen die Indices vertauscht zu sein, da 
der zweite Index die Stelle der Verschiebung und der erste Index die Stelle des Angriffs des Mo- 
mentes angibt. AuBerdem ist nicht unbedingt y;; gleich y;,;. Dies fiihrt nun bei den EinfluB- 
zahlen y;;, welche zur Berechnung der Verschiebung unter dem EinfluB eines Momentes dienen, 
zur Vertauschung der Indices gegeniiber der sonst itiblichen Regel. Erst dadurch erhalt man 
iibrigens eine symmetrische Elastizitatsmatrix. 


Die Verschiebungen y; sind 


n n 
— 2 
5 6 vi = Dd) Kik M, Ve Des Op Or Ge (3) 

k=1 yal 

ES ES Die Neigungswinkel sind 

3 g 
n 

Abb. 2. Biegsame Welle auf starren Lagern. Oh = D5 Vik M, wo? Ds ere Day Bir —, Ww? Py + (4) 

k=1 


Damit kann man fiir das System, das eine Welle auf festen Lagern darstellt, das Gleichungs- 
system ansetzen. In Abb. 2 ist das System zu sehen; dabei bedeuten: Punkt 3 und 4 feste Lager, 
Punkt 5 und 6 rotierende scheibenformige Massen. Es gilt yz; = 0 und Yq — 9. Fiir die Punkte 5 
und 6 lauten die Gleichungen 


V5 = M55 Ms W" Ys + O5g Mg WO” Yg — V55 Os W? Ys — V65 Og © Y » 
Ye = Xgs Ms W Ys + Keg Mg W? Ve — og Os we? Ps — Vo Og © Y » 
Ps = V55 Ms WO ¥5 + V5 Mg W* Ve — Bos O5 w? Ps — Bog Og WM » ©) 
Pe = Vos Ms W ¥5 + eg Mg W Ve — Bes Os W Ys — Beg Og © Ge - 


' C, B. Biezeno und R. Grammel, Technische D ik, 2. Aufl., B P ie. F 
Heidelberg 1953. 1SC ynamik, ., Bd. 2, S. 184 und 185, Berlin/Gottingen/ 


2 Vel. in diesem Buche S. 193. 
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In der Matrizenschreibweise ergibt sich das einfache Kigenwertproblem 


Xs %55 %6 V5 Ves Ms 0 5 
a Vey 1 Mes) es Vag) « 7e8 Ms X6 
OV oe} \ ys Yee Bes Bee =O) Ps 
P6 Yes Yeo Pes Pes 0 —O, pa) 
oder 
231) = Urata Y (6) 


b) Die Ermittlung der EinfluBzahlen unter der Annahme, daB® die Welle un- 
endlich biegesteif und die Lager elastisch seien. In Abb. 3 ist die Welle als starres 
Gebilde auf zwei Federn gelagert dargestellt. Dabei bedeuten Punkt J und 2 Festpunkte, Punkt 3 
und 4 Massen der Lager und Teile der Wellenmasse an den 
Lagern, Punkt 5 und 6 rotierende scheibenférmige Massen und 
Teile der Wellenmasse an diesen Punkten, «3, und «4, alle Weich- 
heiten zwischen den Festpunkten und der Mittellinie der Welle 
an den Punkten 3 und 4. 

Ks gibt auch hier drei verschiedene Arten von EinfluBzahlen. 
Zuerst seien die Verschiebungen unter der Einheitskraft behandelt. Abb. 3. Starre Welle auf elastischen Lagern. 
Ks ist charakteristisch, daf nur die Verschiebungen in den 
Punkten 3 und 4 bestimmt werden miissen. Da die Welle gerade bleibt, ergeben sich die Ver- 
schiebungen in den tibrigen Punkten allein aus den Langenverhaltnissen. AuSerdem sind die 
Verschiebungen in 3 und 4 nicht gekoppelt, da eine Kraft in 3 keine Verschiebung in 4 hervorrufen 
kann und umgekehrt. Es ist unter der Annahme von Fliehkraften in den Punkten 3 bis 6 


I I 
Ses OE (" ae i ae Fe). 
34 o4 


(7) 


hk I 
Yau (Fat 7 Fo t 7 Fo) 


I34 
AuBerdem ist 


Lys 1 
Ys 5 Ps 1, Ya? 
I I (8) 
46 
LS ag mee CA Gap 
In der Matrizenschreibweise wird dies 
l I Ei. 3 
/ ls __ tas 
¥3 os O\f1 O Iss 34 F, (9) 
NO eet WOR SAV I 
I34 4 lang 
eee 
Me 0 l 
NES) Ope (") (10) 
Ys I34 I54 Ya 
Ye is bse 
Isq I34 
Setzt man (9) in (10) ein, so erhalt man 
1 0 
i 0 1 145 pel lag Be 
Xa — Lis Is5 X33 0 i 0 Be ig Ey (11) 
- Tog Toa |\O O%gg/ 1 0 tes Ise F; 
Y6 Tis ae lop ae F, 
I34 Ioq 
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Das zwischen dem Gleichheitszeichen und dem Kraftvektor stehende Produkt erlaubt eine véllig 
schematische Berechnung der «-Werte fiir die starre Welle auf elastischen Lagern. 

Nun folgt die Behandlung der Verschiebungen unter der Wirkung der Momente. Es ist in 
Abb. 3 leicht zu sehen, daB ein in irgend einem Punkt der Welle angreifendes positives Kinheits- 
moment im Punkt 3 die Kraft — 1/I,, und im Punkt 4 die Kraft 1/l,, hervorruft. Die Verschiebung 


in Punkt 3 ist dann — cgq/J,, und im Punkt 4 o44/I,,. Dafiir kann man schreiben 


/‘M, 
Heareotaae 
yi} eo let O Px! \y Leta) Sat aes 

6 
Da (10) auch hier gilt, kann man setzen 
furs aural 
Ys Oia M, 
OT ey ee! ax Thos fe )( 1 1 1 ) M, (13) 
¥5 I5q Ing Us4 [XO X44 Il tl it I M,; 
X6 as bee M, 
\ Iq In J 


Das zwischen dem Gleichheitszeichen und dem Momentenvektor stehende Produkt entspricht 
den Werten y;; fiir Verschiebungen unter der Wirkung von Momenten. 


Es mu8 nun untersucht werden, wie die y-Werte sich darstellen lassen, wenn sie die Neigung 
unter der Wirkung von Kraften darstellen. Am starren Stab sind die Neigungswinkel in allen 
Punkten einander gleich: 


G3 Vase es == Ver (14) 
Sie lassen sich aus den Verschiebungen in den Punkten 3 und 4 darstellen: 


Cg = (15) 


54 I34 


In der Matrizenrechnung wird dafiir geschrieben 


Pa} _ a sib UL 3 

ise re) Fees (16) 
Y5 rl Na 
P6/ ral Ray 


Die Verschiebungen in 3 und 4 ergeben sich aus dem Angriffspunkt der Kraft nach dem Hebel- 
gesetz. Es gelten die Gleichungen (7). Damit werden die Neigungen g; aus (9) und (16) 


ie ele 0/1 oo a 

Pa ee cl ALN as Iq I34 Fy (17 
Ys Se VES LF \ Oo evr aol lan 56 FF, ) 
OR —l I I3q Iq E 


Vergleicht man die Matrizengleichungen (13) und (17), so sieht man, daB die Faktoren vertauscht 
sind; die Zahlenwerte sind jedoch die gleichen, die aus beiden Ausdriicken errechnet werden. 


Es sind noch die Formeln fiir die Neigungen unter der Wirkung von Momenten herzuleiten. 
Auch hier sind alle Koordinaten g; einander gleich, da die Welle ja starr ist. AuBerdem ist die 
Verschiebung in 3 und 4 unabhangig von der Stelle, an der das Moment M; angreift, da man Mo- 
mente am starren Kérper verschieben darf. Die Verschiebungen in 3 und 4 sind 


<7 1 
Ye — 7 tas (Ms + M,+ M;+ M), Ya = 7, Maa (Ms + My + Ms + M,). (18) 
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Mit dem Ausdruck (15) werden die Winkel 


9 


4 


1 
Os a Pe Ps = Pa — (i. (33 + O44) (My 4 M,+M;-+ M,). (19) 


In der Matrizenschreibweise kann man setzen: 


P3 ealgerh M, 

a ae (ay ar ont h Hs ‘ 1 1 --l 1\ | M, 
Ys 84 —l 1jJ\0 a, 1 1 | | M; | 20) 
M6 / 1 1/ M,/ 


Hier ergibt das Produkt auf der rechten Seite der Gleichung zwischen Gleichheitszeichen und 
Momentvektor die EinfluBzahlen f fiir die starre Welle. Der Ausdruck lieBe sich einfacher schreiben, 
es soll jedoch ein schematischer Ausdruck gefunden werden, in dem alle EinfluBzahlen zusammen- 
gefaBt sind. Dazu kann man nun (11), (13), (17) und (20) zusammenschreiben und erhalt fiir die 
Elastizitatsmatrix einer starren Welle auf clastischer Lagerung 


ea \ 
0 1 
ls Iss. 
Ing Igy 
lye L3 / las Lie | 1 1 i 1 
9 el ioe 0 Nome 4 I34 | Is4 I34 U4 Iq 
mg 1 1 : | 0, 1 «Tes 156 | 1 1 if gd BS (21) 
3 Is4 Tyg Iya I34 | Iq I 34 ls4 I34 
il 1 
dag dag 
1 1 
dsp Usa 
1 1 
dg Ugg 


Die Ausrechnung ergibt eine quadratische symmetrische Matrix. Die Elastizitatsmatrix fir die 
biegsame Welle hat zwei Zeilen und zwei Spalten weniger. Fiir die Punkte 3 und 4 sind dort also 
Nullen zu setzen, da diese Punkte ja als starr angenommen waren. So kommt man auf zwei qua- 
dratische Matrizen der gleichen Zeilenzahl. Die gesamte Elastizitatsmatrix ist dann die Summe 
dieser beiden Matrizen: 


2 aes Qieeer = vette 2 


und die Eigenwertaufgabe lautet entsprechend (6) 
1 © 
rae 9) = (QUireeer a Welte) My . (23) 


Diese EinfluBzahlen gelten nun allgemein und sind nicht daran gebunden, daf$ man nur 
die biegekritische Drehzahl (lavalkritische) berechnet. Es wiirde ja auch fiir die Berechnung der 
Kigenfrequenzen der stillstehenden Welle gelten ebenso fiir jede andere Ordnung der Drehzahl. 
Es soll angenommen werden, da die Welle in allen Querschnitten eine kreisférmige Tragheits- 
ellipse habe, also in allen Richtungen quer zur Wellenachse die gleiche Steifigkeit besitze. Wenn 
man die Gehauseelastizitat beriicksichtigt, fiigt man die Weichheit der umlaufenden und still- 
stehenden Teile zusammen. Es ist zunachst einmal angenommen worden, da auch die Gehause 
in allen Richtungen die gleiche Steifigkeit besitzen sollen. Sind sie nun etwa in senkrechter und 
waagrechter Richtung verschieden steif, so kann man den Fall in beiden Richtungen getrennt 
behandeln, da dann eine Aufspaltung der kritischen Drehzahlen eintritt.1_ Bei verschiedener 


1 Vel. Seite 168 bis 171 des in der FuBnote | auf Seite 34 zitierten Buches, 
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Steifigkeit der nicht rotierenden Teile tritt keine Koppelung in den beiden Richtungen ein. Der 
Fall unterschiedlicher Steifigkeit in den rotierenden Teilen (der Fall der ,,ovalen Welle), soll 
hier nicht behandelt werden. 
Wahrend bei der biegekritischen Drehzahl einer biegsamen Welle auf festen Lagern die Welle 
als deformiertes und erstarrtes Gebilde umlauft und dabei die Dampfung sehr gering ist, nehmen 
bei elastischen Gehausen auch nicht rotierende Teile mit ihren Verformungen an der Schwingung 
teil. Sie werden bei der Bewegung laufend gespannt und entspannt und verbrauchen dabei Damp- 
fungsenergie. Man kann aus Erfahrung sagen, dal} eine kritische Drehzahl in dem MaBe weniger 
gefahrlich ist, als nichtrotierende Teile mit ihren Deformationen daran teilnehmen. Durch die 
Berechnung der Schwingungsform und den Vergleich zwischen den Verschiebungen an rotierenden 
und nicht rotierenden Teilen kann die Gefahrlichkeit einer kritischen Drehzahl beurteilt werden. 


3. Das schwingende System in elastischer Aufhangung. Es sollen nun die Punkte I und 2, 
welche die auBeren Teile der Gehause einer Maschine darstellen, nicht mehr unbeweglich sein. 
Damit werden die Massen an diesen Stellen nicht mehr als unendlich groB betrachtet (Abb. 1). 
So tritt aber nun die Frage auf, ob nicht auch noch weitere Massen, die an den Punkten I und 2 
elastisch befestigt sind, an der Bewegung teilnehmen. Damit werden entsprechend dem in Zift. 2 
Gesagten rotierende und nicht rotierende Teile zu einem schwingenden System vereinigt. Dies 
ist bei punktférmigen Massen ohne weiteres méglich. Hier wird eine Drehung im Punkt selbst 
gar nicht registriert; der mathematische Ausdruck fiir die Schwingung im Stillstand ist identisch 
dem fiir die kritische Drehzahl. Scheibenférmige Massen dagegen sind im Stillstand mit Mo- 
menten belastet, welche die Ausschlage zu vergréBern versuchen und mit dem gleichen Vorzeichen 
wie die Fliehkrafte in das Gleichungssystem einzufiihren sind. In der biegekritischen Drehzahl 
ereifen an scheibenférmigen Massen Momente an, welche die Ausschlage zu verkleinern bestrebt 
sind. Diese Momente miissen mit entgegengesetztem Vorzeichen eingefiihrt werden. Dies hat zur 
Folge, daB eine so eingesetzte scheibenformige Masse trotz der beiden hinzukommenden Koordi- 
naten nur zu einem positiven Eigenwert AnlaB gibt; es erscheint bei der Eigenwertbestimmung 
noch ein negativer Eigenwert, der eine imaginare EKigenschwingungszahl bedeutet. 

Es soll nun ein System, wie in Abb. 4 gezeichnet, zugrunde gelegt werden. Dies ist nur eines 
der vielen méglichen Systeme. Indessen soll unter der Verwendung von nicht zuviel Massen- 
punkten das Prinzip erklart werden; die Aus- 
dehnung auf andere Faille ist dann leicht méglich. 
In Abb. 3 bedeutet die Gerade 3—5—4—6 eine 
mit der Winkelgeschwindigkeit ~ rotierende Welle, 
welche elastisch auf den Federn «3, und «4, in dem 
Gehiuse I1—2 gelagert ist. Dieses sei zwischen 
den Punkten J und 2 als biegesteif angenommen, 
was in den meisten Fallen gelten kann. An 
diesem Gehauseblock sei biegeelastisch noch eine 
Masse m, angehingt. Der Gehauseblock J—2 
sei elastisch auf den Federn x ,, und a, aufgehangt. Die EinfluBzahlen der Aufhangung a,, und 
Og Sollen nicht mit den EinfluBzahlen der Welle zusammengerechnet werden. Die EinfluBzahlen 
der Welle seien, wie in Ziff. 2 gezeigt, unter der Annahme fester Punkte 1 und 2 berechnet. 

Man kann fiir die Verschiebungen in den Punkten 1 bis 7 und fiir die Neigungswinkel in den 
Punkten 5 und 6 das Gleichungssystem von Tabelle 1 anschreiben, wenn man in den einzelnen 
Punkten die Fliehkrafte (1) und auf der Welle in den Punkten 5 und 6 die riickdrehenden Kreisel- 
momente (2) angreifen laBt. Die Verschiebungen der Punkte 3 bis 7 setzen sich aus den Verschie- 
bungen des Systems als starrer Kérper auf den Federna,, und a, und aus der elastischen De- 
formation des Systems selbst zusammen. Das erste ergibt sich aus der Lage der Punkte 1 und 2; 
das zweite aus den Kraften, Momenten und Elastizitaten des Systems. Ebenso setzen sich die 
Neigungswinkel y, und @, zusammen. Es sind das die schon in Ziff. 2 beschriebenen Grundsatze. 

Definiert man als Amplitudenvektor ty 

1) 
V2 


Veale (27) 


Abb. 4. Biegsame Welle elastisch in Gehause gelagert, 
das seinerseits elastisch auf dem Fundament aufgehiangt ist. 
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und als Langenmatrix 


( 0» 0-420 20 0 Omen : 
OO.) O20) eset 0 OO 
= [EY eee eee 
hs hs | 
lie lie | 
Ig dy | 
lie lie | 
Ios Ls | 

ae ee tee ara 0 (28) 
ee le 

de he | 
eats 
le lie | 
1 1 | 

Set Os hie 

| 
| 
| 


und als Elastizitaétsmatrix, die in Tabelle 1 angeschriebene Matrix 2{ und als Massenmatrix 
2 
Ms 


( m, 2 
m 
Hh 


m = ms . (29) 


so kann man fiir das Gleichungssystem anschreiben 
y =(L + OAM) y. 


Dies stellt die allgemeine Eigenwertaufgabe dar. Bringt man & auf die linke Seite der Gleichung 
und dividiert man durch © — &, so erhalt man die einfache Eigenwertaufgabe, welche durch 
Iteration lésbar ist.) ? 


+) =(C—QHAAMy. (30) 


@ 


Die Kehrmatrix (© — Q)~ ist leicht zu berechnen, da gilt 


(E—Yt*=C +8, (31) 
wovon man sich leicht tiberzeugen kann. Fihrt man die Multiplikation 
(EF + QA (32) 


aus, so kommt man wieder auf eine symmetrische Elastizitatsmatrix. Auf die gleiche Matrix 
ware man gekommen, wenn man alle Elastizitaten des Systems gleich beriicksichtigt hatte. 

Bisher war angenommen worden, daB die Elastizitat der Gehause x3, und «,, schon in den 
KinfluBzahlen der Welle enthalten seien. Man kann aber das Gleichungssystem fiir den Fall von 
Abb. 4 auch anschreiben, indem man das oben gezeigte Prinzip zweimal anwendet. Festpunkte 
fiir die Berechnung der EinfluBzahlen der Punkte 3, 4 und 7 sind dann die Punkte 1 und 2. Die 
Punkte 3 und 4 sind Festpunkte fiir die Berechnung der EinfluBzahlen der Welle. In diesen sind 
die Weichheiten «3, und a4, nicht mehr enthalten. 


* R. Zurmiihl, Praktische Mathematik, 2. Aufl., S. 133, Berlin/Gottingen/Heidelberg 1957. 
2 R. Zurmihl, Matrizen, S. 152, Berlin/Géttingen/Heidelberg 1950. 
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70 70 1 
ag ot | 9 srt Co Sau To mia 7620 Tan 1p a 2 oy Sup tho ah 0 Thott = 
| ye z @ “Ul rd Fe 6 @ 
°b 0 °G a °h 20 °@ iho 2 2c “wu Mo Ty x 09 Su 9) z iy s vy ey | 
DE 911999,.L, 


42 B. Jager: Eigenfrequenzen und Schwingungsformen von Turbomaschinen Ingenieur-Archiv 


Diese Uberlegungen fiihren auf das Gleichungssystem der Tabelle 2a. Die Gleichung fiir y, 
und y, sind die gleichen, wie in (24) der Tabelle 1. Die Verschiebungen y3, y, und y, setzen sich 
aus den Verschiebungen als starre Kérper, gegeben durch die Lage der Punkte 1 und 2, und aus 
den durch die Weichheiten x33 und a4, gegebenen Deformationen zusammen. Die Verschiebungen 
der Punkte 5 und 6 sind beschrieben durch die Lage der Punkte 3 und 4 und durch die Verbiegung 
der Welle selbst. Die Neigungswinkel p, und q sind durch die Stellung der Punkte 3 und 4 und 
durch die Deformation der Welle definiert. 


Amplitudenvektor und Massenmatrix bleiben gegeniiber (27) und (29) unverandert. Die Langen- 


matrix ist 


(ae 02700 bed) get Um 
0 0 eae kt ey 
Iog hs 
fog Ga a 
Ly, Ls 
Ned iy ea 

Iys 135 
or Pe Le Pin 0 (36) 
1 Isg | 
a ee 
Is4 Is4 
ewes We at 
12 lis 
1 iL 
(all Ove oe 
sq sq 
1 1 
Nae ae ee 
I34 Is, 


Die Elastizitatsmatrix umfaBt wieder alle mit jm? behafteten Glieder (Tabelle 2b). Ihre Auf- 
stellung ist einfacher als im Fall der Tabelle 1, da die Federn a3; und o4, nicht in die EinfluB- 
zahlen der Welle hineingerechnet zu werden brauchen (34). Dafiir hat die Matrix (€ — Q)7 einen 
etwas komplizierteren Aufbau. Sie lat sich nicht mehr so einfach 
anschreiben wie im Fall (25), kann aber noch in allgemeiner Form 
angegeben werden: Ausdruck (35) in Tabelle 2b. Durch Ausmulti- 
plizieren von (© — %)-+9{ kommt man auch hier auf die eigentliche 
Elastizitatsmatrix des Systems. Die Ermittlung der EKigenwerte ge- 
tn, schieht wieder wie bei (30). Die Eigenwertaufgabe hat auch hier die 
Abb. 5. Laufer mit elastisch Form der Gleichung (30). 
ante sheng tees CcBeure: Es soll noch auf eine Méglichkeit hingewiesen werden, eine ver- 
einfachte Rechnung im Falle eines Laufers durchzufithren, bei dem 
mit wesentlichen Verbiegungen im Drehzahlbereich nicht zu rechnen ist. Dies gilt vor allem, 
wenn man die Aufhaéngeschwingungen einer Maschine berechnen will. Das System Abb. 5 soll 
gegeben sein. Dabei ist unter den Punkten I und 2 ein Gehause auf den Federna,, und o» zu 
verstehen; der Punkt 3 stelle einen relativ steifen Laufer dar. Das Gehause soll noch ein aqua- 
toriales Massentragheitsmoment @,, besitzen. Dazu gehért dann ein Moment 


Pane 2 
Myy = Ojy © Pp - 
Das riickdrehende Kreiselmoment des Liufers ist 
2 
M, = — 0, w*g, . 


Driickt man nun die Neigungswinkel durch die Verschiebungen in I und 2 aus: 


1 
C05, UD Te 


12 


Wey) 
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Tabelle 2b 
l, l l 
og 0 28 ——10 ab 26 ly at Oy 
iat Ls 11 1 : 1 f Ay i Oy ik oael i = ie 
0 oD) ao op na Ooo bs O59 at Os ae Ooo G22 Sea 
he hy be ~ Le Liga tee he he 
0 0 One 0 fas ep ae Lig O35 0 __ Og X33 
Ia Iq lsq Iq 
eter 0 0 ee lessee > Sn Oy 4 
Iq Iq Iq Isq 
0 0 0 0 O55 O56 V55 V65 
f : ¢ 0 65 %6 0 56 Y66 
0 0 0 0 0 0 dey 0 0 
0 0 0 0 55 V56 0 Ps5 Bre 
0 0 0 0 ie Ves 0 Bes pe 
( 1 0 0 0 Ui 0 
0 if 0 0 Ge W 0 
I, l 
rie ri 1 i--"01 oO 
12 12 
E; l 
= — 0 La 00, O's LO 
lis hye 
[yeh (Bye! Teelne ey hf I toe 
7 7. is = 7 i 4 7 i ‘ie fe 1 0 0 0 
(G ais gy = 12 "34 12 "34 12 “34 12 “34 34 34 
I, Ing Tn, Ls6 his Iyg | lig Is6 Ing Ty6 0 il 0 0 
lio lay lie bs li Iy4 ly» Is, 34 154 
ise l 
23 = 0 0 OVROM ear O. 
ie hie e 
oyiigy Bel Eglo digest 1 1 
= = Wy x0) I 
li» I34 lio 34 li Isy ly» Io, lo, ley 
ee ee ae ee a 1 1 
: ata gt We =e LOS 20. 80 
lis I34 lig 54 li lo4 he Is4 Is4 toy 
so erhalt man fiir die Momente 
1 1 ; 
My, = im O1, 0 (¥2 —41) » M, = — me O, w (Yo — 91) 
und fiir die Krafte in den Punkten J und 2 
0d. Or, 5 i @ Ory 
P, =a? ma Te (V2 ¥1) o “ery P, —= we oe rea i. (Ye =) . 
12 12 12 12 
Mit 3,, = O,,/l?, und J; = O,/I?, kann man folgendes Gleichungssystem anschreiben: 
5 Ne 
Yr = My (mm + Dyg — Vy) © ¥y — O41 (yg — Vg) © Yq + O41 ie M3 0" Ys » 
l 
Ya = — gg (Ppp — Oy) W? Vy + Og (Mg + yg — Dy) w Yo + Og9 is Ms WO” ¥5 5 
is l 
¥3 = = V4 ae fa Vo + X35 Ms WY » 
12 12 


Die auch auf erhalb der Hauptdiagonalen besetzte Massenmatrix ist 
m, + Oy, — Us — (Fy, — Ds) Os 
Ne= | —(,.— 9) m, + 0,,— Vs 0 
0 0 Ms 


(34) 


(35) 


(37) 
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Die Elastizitatsmatrix ist 


fae 0 fas all 
12 | 
ee I : (38) 
0 Xp ne Oo9 
0 0 0 
und die Matrix (© — %)* ist 
1 0 0 
©=94+=64y2=( "9 (39) 
lbs is i) 
lig he 


An den Stator 1—2 kann man noch beliebig Massen anhangen. 
Die einfache Eigenwertaufgabe lautet nun wieder wie bei Gleichung (30) 


+) =(—QHAAMy. 


4. Das schwingende System frei im Raum. Man stelle sich jetzt vor, daB die Federn a,, und Oy. 
sehr weich seien im Vergleich zu den iibrigen im System vorkommenden Elastizitaten. Eine zu 
diesem System gehérende Maschine wiirde im unteren Drehzahlbereich (bei oder unterhalb von 
30% der Nenndrehzahl) Aufhangeschwingungen durchfahren. Diese sind im allgemeinen unge- 
fahrlich. Die eigentlichen stérenden kritischen Drehzahlen im héheren Drehzahlbereich erscheinen 
aber in der Rechnung erst als héhere Eigenwerte. Man miiBte also bei einer mehrfachen Durch- 
fiihrung der Rechnung jedesmal erst die beiden unteren Eigenschwingungszahlen berechnen, um 
die wirklich interessierenden Werte zu bekommen. AuBerdem wiirden bei der Rechnung zu den 
groBen Werten fiir die Weichheit der Aufhangung der Maschine die kleinen Werte fiir die Elastizitat 
des Systems selbst hinzukommen. Die Rechnung miBte deshalb mit groBer Genauigkeit durch- 
gefuhrt werden, um die héheren Eigenwerte noch geniigend genau zu ergeben. 

Es ist einfacher, in einem gewissen Abstand von den Aufhange-Higenfrequenzen das System 
als frei im Raum zu betrachten. Es ist wohl eine Transformation des Gleichungssystems erforder- 
lich. Es scheint zuerst die Schwierigkeit zu bestehen, die Einflu®zahlen fiir diesen Fall zu 
berechnen. Diese werden aber wie in Ziff. 2 und 3 unter der Voraussetzung berechnet, dafB die 
Federn o,, und a. unendlich steif seien. 

Es soll zur Erklarung des Verfahrens das Gleichungssystem Tabelle 1 zugrunde gelegt werden. 
Dividiert man die erste Gleichung fiir y, durch «,, m, w? und die zweite Gleichung fiir y, durch 
Xoo My, w*, so erhalt man 


se Noa Ye alts ak TS) ess las We oe Paty TO ats LOY eas eae 
4, mM, w? : (Re Te Peer lhe, ite V1 he TP Ps Lig) Ty v6? 
__¥2 (ore Ihe Ms a La m4 le m7 1 0; i 0, 
Oly Ms WW” 0 hy ms Yar ly ms v4 25 ly mg Ys lz mg ae Ig my oe 


LaBt man o,,—> co und o> oo gehen, so werden die beiden linken Seiten der Gleichungen zu 
Null, und man kann fiir y, und y, schreiben 


ee om Tiley Un Tiap lez mm, I Oy. 1 0, 
a le my oe Iya my a ' ye m se hy my, we Ig my ae 
re (bs Ms Va ™4 ly, m7 1 O; i OP 
vA lig Ms es la my a Ty me Ks ly ms a ly, Ms AS 
In der Matrizenschreibweise ergibt dies 
ae tite In, mM, ae ly, mz 1 0; 1 @O, (73 
e) aed ly my he my [ig my hy my hh, m | (40) 
\Ye they THe Iya ™, ae liz my l @O; 1; V7 | 
Ii ms he my Io My hy ms Ig ms eA 
Pe 


oder kurz 
Yie Sa nue Yg~7 . 
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Tabelle 3a 
In, lig 
(92) ( i Ly (2s Asq O35 Age 0 V53 vss (m, 0) ( v0) 
Tog La 
Ya i iss Hyg Oyq %y5 Age 0 Ysa Yea 4 Ya 
Ins ys 
Xs im "Ex ; X53 Aq X55 Ase 0 V55 65 Ms ¥5 
=x l I yar: 2 
SiGel [eae a ema re i, fe + @ Xs Xe 5 Age 0 Vs6e 66 M6 Y6 (41) 
ee ba hz 0 0 640 <0 
V7 % E; Okan 0 0 my Yq 
il 1 
Ps ae G ie Y53 Ysa V55 56 0 Bss  Bse —0; M5 
1 1 
%6 a a cm Yes Ves Ves Vee 0 Bes Bes 0 —O; Po 


Fir die ibrigen Koordinaten ist das Gleichungssystem in der Matrizenschreibweise auf Tabelle 3a 
angegeben. Fihrt man nun in dieses System die Gleichung (40) ein, so kommen die Koordinaten y, 
und y, nicht mehr vor, und man erhalt ein Gleichungssystem von der Form: 


| Yung = — & Vm Ygzz + ©? Ug. 7 Mege7 Yge7- 
Bringt man & &, auf die linke Seite der Gleichung, so erhalt man 
{f 
we (E + @ Vm) Yge7 = Use 7 Megz7 Ya=7 
oder 


1 
2 KY 327 = Usz7 Mesz7Ys=7- 


Die Matrix § ist auf der Tabelle 3b angeschrieben. 


Tabelle 3b 
l. l 
l 0 28 13 
lis li 
I l 
if 24 it 
lie dye 
1 Lys Ls Z : 
Fey Ths Lg Mz Ipg mg Las Ms log Mg Ion M7 1. © 1 O¢ 
] Le my, ly my by, my Ig my I. my ly. my Lig my 
@ = Sey) et eee ae (42) 
le bie Lg mg yy my lis ms lig me yz my ay @; 1 9 
1 lay La Ig Mz tg mg yy my ho Mz Iya My lig My Iyy Mz / 
lye Lys 
1 il 
1 i, 
le is 
1 1 
0 1 ——— 
lie ie 


Dies stellt das allgemeine Eigenwertproblem dar. + 
Durch Bilden der Kehrmatrix fiir den Ausdruck (42) kommt man auf die einfache Aufgabe 


= Yge7 = R72 As Megz7 Ysa - (43) 


Das Gleichungssystem hat nun zwei Koordinaten verloren. Praktisch lauft die Ua eus auf 
eine Transformation der Elastizitatsmatrix durch Vormultiplizieren mit dem Kehrwert §% einer 
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aus Langen- und Massenverhiltnissen entstandenen Matrix hinaus. Um die Schwingungsform 
vollstandig zu bestimmen, mu8 man die Koordinaten y, und y, aus den Gleichungen (40) be- 
rechnen, wenn die Werte fiir die iibrigen Koordinaten bestimmt sind. 2 

Liegt ein Gleichungssystem von der Form der Tabelle 2a vor, und soll dieses als frei im Raum 
betrachtet werden, so mu die Elastizitatsmatrix ermittelt werden, fiir welche die Punkte J und 2 
Festpunkte sind. Um von der Annahme der Festpunkte in 3 und 4 auf die Festpunkte in J und 2 
zu kommen, mu die in Ziff.3 beschriebene Transformation, namlich die Gleichung (32) an- 
gewandt werden. Man hat mit 2{ aus (34) und { aus (36) ohne die ersten beiden Zeilen und Spalten 


(Tabelle 2b): 
| \( ive la _% _%s) 
| | 
le 0. a 
0, 19) 0 0° 0 0 9 0) Ooh tual ee ee eee 
eee LS ie raed Wey Lie eae Pagers re 
I | | 
Lys 35 (Jd 0!) 0 ORs Xsg 0 W565 Y65 
Iy4 Is4 | | 
I l 
a 2 7 | it 0 0 1 Ges X66 0 156 Y66 
34 34 | 
| 
Ona l yew w 0 Sj FORMS 
1 | | 
Pays aes 1 0 0 V 55 56 0 Bss Bse 
34 34 | | 
1 ef | 
—- ;!0 1c} Oe Orecee Yee 0 Bes Bes 
34 34 | | 


In den Ausdriicken fiir y, und y, treten Bedingungen auf, welche nur noch Langen, Neigungen, 
Massen und Drehmassen enthalten. Dies stellt die Bedingungen dar, die fiir einen Punkthaufen 
im Raum gelten, namlich daB der Schwerpunkt des Haufens keine Verschiebung erfahren kann, 
was auch immer die einzelnen Punkte fiir Bewegungen ausfiihren. Damit kann die Rechnung 
kontrolliert werden. Bei der errechneten Schwingungsform miissen die an den einzelnen Massen- 
punkten angreifenden Fliehkrafte und Momente die beiden Bedingungen erfiillen: 


a) Die Summe aller Krafte, die am System angreifen, mu Null sein; 
b) die Summe aller Momente mu Null sein. 


Zu a): Als Krafte kommen nur Fliehkrafte in Frage, die Massentragheitsmomente ergeben 
keine Beitrage dazu. Es gilt also 


= my; =0 
i=l 


oder als skalares Vektorprodukt geschrieben 


(m,mz...m,!0...0)]y, | =0. (45) 


Pn 


Zu b): Die Momente setzen sich aus folgenden Anteilen zusammen: Aus den Momenten der 
Fliehkrafte um einen Punkt; dazu wird eine Langskoordinate x eingefiihrt; es ist dabei gleich- 
giiltig, wo man den Nullpunkt fiir « wahlt. Um keinen Vorzeichenwechsel zu haben, ist es zweck- 
maBig den Nullpunkt in die ganz links liegende Masse zu legen. Das Moment der an den einzelnen 
Massenpunkten angreifenden Fliehkrafte ist dann 


n 
mi Po 1 ie 
i= 
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a eee die Momente der Drehmassen hinzu, wobei die Werte QO; fiir rotierende Massen wegen 
es Kreiseleffektes negativ, die Werte fiir nicht rotierende Massen positiv eingesetzt werden: 


M = 5 6; w? Qi; . 


Dies ergibt 


2 mi UV+ > O:9; = 0. 
i= =i 


Als skalares Vektorprodukt kann man dafiir schreiben 


(m,%,...m,x, | @,...0,)|--|=0. (46) 


Pn 


5. Zur Rechenpraxis. Es ist in der bisherigen Ausfiihrung versucht worden, zu zeigen, daB 
sich die Gleichungssysteme, auch wenn sie noch so umfangreich sein mégen, bei allen behandelten 
Fallen tibersichtlich schreiben lassen. Die Rechnung kann also schon unmittelbar nach dem Auf- 
stellen der Gleichungssysteme auf véllig schematische Rechenoperationen iiberfiihrt werden. Diese 
sind ohne grundsatzliche Uberlegungen auszufiihren und kénnen von einer Rechenhilfskraft auf 
einer elektrischen Tischrechenmaschine oder bei gréBerem Umfang auf einem elektronischen 
Rechenautomaten durchgefiihrt werden. 

Bei der Eigenwertbestimmung gibt es verschiedene Verfahren. Bei den vorliegenden Aufgaben 
sind Iterationsverfahren besonders geeignet. Diese fiihren zuerst auf den gréBten Eigenwert, 
also die unterste Kigenschwingungszahl. Dann kann man auch die héheren Eigenschwingungs- 
zahlen bestimmen. 

Die Schwierigkeit der Rechnung liegt darin, zu einer gegebenen Maschine ein System zu finden, 
das die gewiinschten Eigenwerte geniigend genau ergibt. Dafiir ist eine gute Kenntnis der Ver- 
teilung von Masse und Elastizitat notwendig. Man wird mit méglichst wenig Koordinaten auszu- 
kommen suchen. Dabei ist zu bedenken, da im allgemeinen nur ein bestimmter Drehzahlbereich 
interessiert. Man kann also sehr weich am Gehause befestigte Massen in der Annahme weglassen, 
daB sie im interessierenden Drehzahlbereich schon iiberkritisch sind. Ebenso ist zu untersuchen, 
ob in diesem Bereich schon Verschiebungen zwischen relativ steif verbundenen Gliedern auftreten. 
Indessen ist der Genauigkeit der Rechnung durch eine geniigende Anzahl von Koordinaten keine 
Grenze gesetzt. 

Eine gréBere Schwierigkeit liegt in der genauen rechnerischen Erfassung der Elastizitat der 
einzelnen Glieder. Dies setzt nicht nur Kenntnisse der neueren Festigkeitslehre voraus; vielfach 
gibt es noch keine rechnerische Lésung fiir ein Problem. Dann hilft nur noch ein statischer Versuch. 

Es soll nun noch gezeigt werden, wie man mit der Aufstellung des Systems beginnen kann, 
wenn bei der zu untersuchenden Maschine ein relativ steifes Gehduse als Mittelstiick vorliegt, 
an dem die iibrigen Teile rotierend oder nicht rotierend befestigt sind. 

Diesen Block nimmt man als unendlich steif an, wie dies bisher mitdem Se 2 
Strang 1—2 geschehen ist. Man mu nun ein Ersatzschema suchen, das 


‘5 . “ Q, Qe 
sowohl fiir eine Auslenkung in y-Richtung als auch fiir eine Drehung um 
G o x O * Oa P. Abb. 6. Reduktion 
den Winkel py die Masse reprasentiert. Man kénnte nun die Gesamtinasse Fie Cacia ie 


im Schwerpunkt des Gehauseblocks vereinigen und noch das aquatoriale auf zwei Punktmassen, 
Massentragheitsmoment dazu angeben. Da aber in den Festpunkten 1 
und 2 zwei y-Koordinaten gegeben sind, wird man Masse und Massentragheitsmoment durch zwei 
Punktmassen darstellen. Es sind die Gesetze der Bewegung des starren Kérpers namlich drei 
Bedingungen zu erfiillen: 

1. Die Gesamtmasse muB erhalten bleiben; 

2. der Schwerpunkt mu an der gleichen Stelle bleiben; 

3. das Massentragheitsmoment muB erhalten bleiben. 


Dies ergibt nach Abb. 6 folgende Bedingungen: 
m, +m,=)>m, Dida = Ns Ay, Tit atm, da ==8 >, Our. (47) 
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Dabei ist r der Abstand der Schwerpunkte der einzelnen Teile vom Gesamtschwerpunkt S,, der 
Massen I—2. Aus diesen Bedingungen ergeben sich durch einfache Umrechnung folgende For- 


meln, wobei man a, beliebig annehmen kann: 


mM, = + uae eo? m, = > m—m,, ti es Ay « (48) 
ym T Ons 
Die Durchfiihrung der Rechnung geschieht mit einer Tabelle, wie sie im folgenden Beispiel an- 
gewandt ist. 

6. Rechenbeispiel. Es sei der in Abb. 7 gezeichnete Axialverdichter als Beispiel gegeben. Die 
Gehiuse a, b und ¢ seien so steif, daB sie fiir die zu berechnenden Schwingungen als starrer Block 
1—2 betrachtet werden kénnen. In den Gehausen a und c sei ein Laufer gelagert, der aus einer 
steifen Trommel mit relativ weichen End- 
scheiben und Zapfen bestehen soll. Die 
elastischen Verschiebungen der Gehause- 
TT Lagerstem naben mit den Lagern gegeniiber den 
mit 1WArmen AuBenringen der Gehaiuse werden in der 
Rechnung beriicksichtigt. Damit hat man 
zur Berechnung der EinfluBzahlen des 
Laufers unter der Annahme festgehaltener 
GehauseauBenringe das Schema von Abb. 8. 
Trommeltduter Die Federn «,; und xg, der Punkte 5 und 6, 
fr $5 45,5 = welche masselos sein sollen, bestehen aus 
[ esl a allen zwischen der Laufertrommel und den 
AuBenringen der Gehause befindlichen 
es eS Oe see iT ae Weichheiten. So kénnen fiir den Schwer- 


0,=19,76 ,-3000— Ug =28 em 
Ugg =24 em 


U5 =52 CM 


= Sip= 4448 


Abb. 7. Beispiel: Reduktion einer Maschine auf ein Ersatzsystem Abb. 8. Beispiel: Schema zur Berechnung 
(Langen in cm). der EinfluBzahlen des Trommellaufers, 


punkt der Laufertrommel mit den beiden Koordinaten y, und g, die EinfluBzahlen berechnet 
werden. Am Gehause c sei noch ein Tragarm befestigt, in dem eine sehr steife Antriebswelle ge- 
lagert sein soll. Die Welle mit dem Tragarm ergibt den Massepunkt 4. 

Der ganze Verdichter sei auf dem Priifstand so weich aufgehangt, daB er im normalen Drehzahl- 
bereich als frei im Raum befindlich betrachtet werden kann. Von den vier Massepunkten mit 
den fiinf Koordinaten wird man also nach der Transformation (Ziff. 4) auf drei Koordinaten kom- 
men. Da eine Kreiselmasse vorliegt, wird man zwei positive und einen negativen Kigenwert er- 
halten. Dies entspricht zwei kritischen Drehzahlen. 

Die Aufgabe enthalt also alle wesentlichen Gedanken; sie ist aber so einfach abgefaBt, daB sie 
leicht nachgerechnet werden kann. 


Tabelle 4. Daten der Teile des Gehduseblocks 1—2 


| | 
Hee te s ms r = S,;— Sy, mr® Oge | mr*+ Og, 
| = Ss? | cm kg s? cm em kg s? em kg s? 
a 0,011 1251 7.307133 Lin le—=3 2036 I Asy333) isa 13,23 
b | 0,032 41,5 1,3280 | — 2,98 0,284 Dal 5,98 
© | 0,016 _ 72,7 | 1,1632 28,22 12,742 2,9 15,64 
Summe | 0,059 2,6242 3 34,85 
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Zur Reduktion der Gehause a, b und ¢ auf die beiden Massen I und 2 wird die Tabelle 4 aufge- 
stellt. Darin bedeuten s (cm) die einzelnen Schwerpunktsabstiande von einer Bezugsebene, und 
Sj. ist der Abstand des Gesamtschwerpunktes von der gleichen Bezugsebene. Ferner ist r der 
Abstand der einzelnen Schwerpunkte vom gesamten Schwerpunkt S,,. Wie in Ziff. 5 beschrieben, 
werden nun die Massen m, und m, und die Lange a, berechnet. Die Lange a, wurde zu 30 cm 


- angenommen. 


_ Sms __ 2,624°3 


Sm = 09 = 4448 cm. 


Sie 


Die Anwendung der Formeln (48) mit a, = 30 cm ergibt 


a A kg 5 
My T = se 0,0234 eae 
34,85 


0,059 
m, = 0,059 —0,0234 = 0,0356 52 3, 
cm 


__ 0,023°4 


a = 99956 30 = 19,72 em. 


Tabelle 5. Daten der Punkte 3 und 4 


m™m 
: ke 
Bezeichnung 2 2 em | cm kg s? 
3 0,025 439° | oa 
4 0,020 104,77 | 0 


Die Lage der Punkte 1 bis 4 ist in Abb. 7 aufgezeichnet. Daraus folgt 


bg == 49572, ei Ls 
Le = 18.54, Ly = 30522, 
Ei —— 79,94, 
hs _ 0,373, bos _ 0,627, Wes _iyaiiyy, a1 nse 
lie he my, IMs 
fa = 1.608. Ion _ 0.608, Tae (862 ; PET 55% 
lig hye aT His 
BES) 090 - Os _ 58 088, 4 _ 90,744. 
he ae WG 
Die Matrix § errechnet sich nach dem Ausdruck auf Tabelle 3b: 
( Ing Lys 
1 0 O | he lis “ Ms _ TT, 1 0, 
i I 1g My 1g My ly, my 
G=T0 1 0)4 a ; 
| he he | bes, Tis ha m4 A QO, 
i gi eee BY Iya Ms Iho Mz ly Mg 
| hi» i) 
‘1,424 0,298 0,070; 
ee UAE ee 0,440 —0,342 1,180 
R=10 1 0} +1 —0,608 1,608 oor eee ae =| 0,372 3,419 —3,604 ]. 
@ Wa —0,020 0,020 : $ ; 0 0,034 0,940 


Die Kehrmatrix §~ ist 
0,7179 —0,0597 —0,282'6 


St = | —0,0752 0,288 0 1,1098 
0,0027 —0,0105 1,023°7 
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Die Massenmatrix Jt ist 


0,025 0 0 
Maa =| 0 0,020 0 
0 0 —2,100 
Es sind folgende EinfluBzahlen gegeben: 
cm ech 
oss = 14,5 - 10°F Hep = 7,6 10-8 X44 = 51-10% F. 


Nach dem Schema Abb. 8 werden die EinfluBzahlen fiir den Punkt 3 berechnet: 


& i =~ Is6 I56 G ‘ Is6 I56 
Ysa P33 gerks 0 Xe Ty5, i 


I56 Is Is6 I56 


— 


24 28 24 
i. 52 52 14,5 O\ [| 52 52 Wee 5,2916 —0,0500 roca 
ilwey "A 0 7,6 —0,0500 0,0082 
oe Ss 52 52 


Setzt man in der Elastizitatsmatrix den Punkt 4 anschlieBend an den Punkt 3, und bringt man die 
Koordinate gy, zuletzt, so ergibt sich die Matrix 


‘ 5,2916 0 —0,0500 
Us4= | 0 Sil 0 | Omer 


—0,0500 0 0,0082, 


Die Ausgangsmatrix fiir die Kigenwertbestimmung ist 
0,09532 —0,06089 0,08024 
K As Mts. = | —0,01128 0,29376 —0,02698 | 10-°. 
—0,00087 —0,01071 —0,01732 


Es ergeben sich folgende Kigenwerte und Eigenvektoren: 


Ay = 0,29819 An = 90,090°97 ‘V5 —0,313'1 ‘V5 / 18,555 
wi = 3,353 - 108 wir = 10,992°6-108 | | ¥,) = 1,000°0 |; | = 1,000 
Or — 1831s oy = 3316 1/s P3/ 1 —0,033'1 3/11  \—0,248 
Am = —0,0161 
Die Errechnung der Koordinaten y, und y, ergibt nach Gleichung (40) 
ele ee MLE SS Be Ss 
°") Sees Ig my Ip my hj. my a 
2/1, 1 hs Ms lym 1 : 
hg Mg Ia my he m, Ps 
a oa ce —0,342 ae Bey 0,518'8 
0,398 1,375 —1,795 7. ia ) 
Ye —0,033'1 1,309°8 
18,555 


yy 0,440 —0,342 1,180 —T,529 
cana 1,000 | = 
0,398 1,375 —1,795 
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Die Schwingungsformen sind in Abb. 9 aufgezeichnet. Die Probe nach Gleichung (45) und (46) ist 
fiir die Schwingungsform I: 


(m, m, m, m, 0)p; =0 


0,518'8 
—1,3098 
(0,035.6 0,023.4 0,0250 0,0200 0) —0,3131 | = 9,000008 ~ 0, 
1,000°0 
—(),033 1 
(0 mghy mghgs ml, —O@ ) yr = 0 
De | 
—1,309°8 
(0 1,163'4 0,4635 1,5988 —2,100°0) —=(.213'1 = 0,00063 ~ 0, 
1,000°0 
—0,033'1 
Neigung im 
Punkt 3 
=i 
s 
ie ; 
7 
Neigung im 
| Punkt3 
-10 
eS 
| OF b 
70 


Abb. 9. a) Beispiel: Schwingungsform1I; b) Beispiel: Schwingungsform IT. 


fiir die Schwingungsform II 
=—(-529 


—9,205 
(0,035°6 0,0234 0,0250 0,020°0 0)]| 18,555] =0,00044 ~ 0, 
1,000 
—0,248 
—7,929 
—9,205 
(0 1,1634 0,4635 1,5988 —2,1000)] 18555] =0,0107 x0. 
1,000 
—0,248 
Die auf Seite 490 des Buches von Zurmiihl! erwahnte Bedingung der Orthogonalitat kann auch hier 
zur Kontrolle der Schwingungsformen I und II angewandt werden: 
pM yn =0, 


1 Siehe FuBnote 1 auf Seite 40. 
4* 
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also 
0,035°6 0\ (/—7,529 
0,023°4 —9,205 
(0.5188 —1,3098 —0,3131 1,000°0 —0,033'1) 0,025°0 | 18,555 | = 0,000505 ~ 0. 
0,020°0 | 1,000 


lo —2,1000 ) |_o.248 ) 


Es soll noch eine Eigenart des Systems, das als frei im Raum betrachtet wird, untersucht werden. 
La@t man die Masse 3 mit ihren beiden Koordinaten zu null werden, so hat man ein Dreimassen- 
system frei im Raum mit einem Freiheitsgrad. Das Quadrat der Eigenschwingungszahl ist 


1,4 \2 m 1, \2 m 1 
= |) ee a | 
ie i a (4) | a (7) = O44 My” ee 
Waren die Punkte I und 2 fest, so hatte man gehabt 
1 
2 
Em reais (50) 


Da alle Ausdriicke in der eckigen Klammer der Gleichung (49) nur positiv sein kénnen, ist der Wert 
fiir w? in Gleichung (49) auf alle Falle gréBer als in Gleichung (50). In dem vorliegenden Beispiel 
ist der Wert des Ausdrucks in der eckigen Klammer gleich 3,419. Die Eigenfrequenz erhéht sich 
also auf das 1,85fache. 

Streicht man die Koordinaten ¢~, und y,, so erhalt man fiir die Masse 3 den Wert der eckigen 
Klammer in (49): 1,424. Die Eigenfrequenz erhéht sich also auf das 1,20fache. Es ist damit be- 
wiesen, da} man schon fiir eine einigermaBen genaue Rechnung die Maschine als Ganzes unter 
Beriicksichtigung der Aufhangung betrachten muB8. 

Die Schwingungsform I kommt durch die elastische Verformung von nicht rotierenden Teilen 
zustande. Bei einigermaBen gutem Wuchtzustand der umlaufenden Teile ist die zu ihr gehérende 
kritische Drehzahl ungefahrlich. 

Die Schwingungsform II dagegen besteht in einer starken Auslenkung des Trommellaufers. 
Fiir diese Schwingung sind die EinfluBzahlen a;; und ag, maBgebend, welche, wie schon erwahnt, 
alle Weichheiten zwischen der Laufertrommel und den AuBenringen der Gehause darstellen. Je 
gréBer nun der Anteil der umlaufenden Teile an der Auslenkung ist, um so weniger gedampft ist 
die Schwingung. Wenn also die Gehause relativ steif sind und die Auslenkung des Laufers im 
wesentlichen durch die Elastizitat der Zapfen und Endscheiben der Trommel zustande kommt, ist 
die Schwingung im Drehzahlbereich unbedingt zu vermeiden. Stellen dagegen die Zahlen «,;, und 
Xgg die Weichheit der Gehause dar, sind also die umlaufenden Teile sehr steif, so ist die Schwingung IT 
ebenso ungefahrlich wie die Schwingung I. Es ist dabei vorausgesetzt, da® die Abstande zwischen 
den Schaufelenden und dem Stetor groB genug sind, um ein Anstreifen in der Resonanz zu ver- 
meiden. 


(Eingegangen am 27. Marz 1958.) 


Anschrift des Verfassers: Dipl.-Ing. Berthold Jéger, Stuttgart-Wangen, Kemptenerstr. 22. 
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Die Beschleunigungsanderung 
I. Mitteilung 
Von W. Meyer zur Capellen 


J. Einleitung. Will man fiir das geradgefiihrte Abtriebsglied eines Getriebes die Extrema der 
Beschleunigungen b ermitteln, so findet man zunachst die zugehérige Winkelstellung « des Antriebs- 
gliedes, wenn man die Beschleunigungsanderung db/dt gleich Null setzt. In gleicher Weise wird 
man auch bei drehendem Abtrieb verfahren. Hierzu bedarf man aber nicht der Beschleunigungs- 
anderung eines einzelnen Punktes sondern der des ganzen Getriebes, d. h. es wird notwendig sein, 
die GréBe allgemein bei der ebenen Bewegung zu definieren. Da die dritte Ableitung des Weges nach 
der Zeit auch, wohl zuerst von Melchior, als Ruck bezeichnet worden ist!, wollen wir uns, auch wegen 
der Kiirze des Ausdruckes, dieses Wortes bedienen, und unsere Aufgabe sei, Satze iiber den 
Ruck i = db/dt aufzustellen und diese 
auf die Ermittlung der Beschleunigungs- 
extrema in Getrieben anzuwenden. 

In der Getriebetechnik wird beilaufig 
auch der Ruck fiir einen anderen Fragen- 
komplex verwendet: Ein Getriebe ist 
stoBfrei, wenn keine Spriinge im Ge- 
schwindigkeitsverlauf auftreten, also die 
Beschleunigung nicht unendlich wird, 
und es ist ruckfrei, wenn keine Spriinge 
im Beschleunigungsverlauf auftreten, 
also die dritte Ableitung (hier der Ruck) 


nicht unendlich wird. 


I. Schiebung und Drehung 


2. Erweiterung des Satzes von Euler 
auf den Ruck. Ist die Bewegung eines 
Getriebegliedes durch die zweier Punkte 
A und B gegeben, und haben A und B 
die Geschwindigkeiten 4, Dg bzw. die 
Beschleunigungen 6, und bg, so gilt 
(Abb. la) nach dem Satz von Euler 


mit AB =r einerseits 
bp =v=v4t+0g4=04+o Xt, (la) 


andrerseits, durch Differentiation ge- 
wonnen, 


bg =b=b, +5544 ¥,pa a Re 
=b,+éexr—tro, (tb) 


. Abb. 1. Der Satz von Euler arn Beispiel einer Zweipunktfibrung 
wobei @ md ¢— oO — do/dt (Punkte a) fiir die Geschwindigkeiten, b) fiir die Beschleunigungen, c) fiir den Ruck. 


bedeuten Ableitungen nach der Zeit) die 


momentane Winkelgeschwindigkeit bzw. — Beschleunigung darstellen (Abb. 1b). Differentiiert 
man nun diese Gleichung nach der Zeit t, so folgt 
b=byte Xt +e xt—to*—21eo (2a) 


1 P, Melchior, Z.VDI. 72 (1928) S. 1842; ferner H. Finkelnburg, Maschinenbau/Betrieb 1935, Beilage Reu- 
leaux-Mitteilung, S.520; fee Peneee von H. Blaise, F. Maul und Entgegnung H F Ei inkelnburg, Maschinen- 
bau/Betrieb 1936, Reuleaux-Mitteilungen S.220.—Dem Vorschlag von H. Blaise, die GréBe db/dt ,,Schwellung 
zu nennen, kann ich mich aus sprachlichen Griinden nicht anschlieBen, da zusammengesetzte Ausdriicke wie 
,,Polschwellung“ und , Schwellungspol** unschén wirken gegeniiber ,,Polruck“ und ,,Ruckpol“ (vgl. u.). 
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oder, dat =@ XY und € Xt =—1WE ist, auch mit h = R 
f,p= R= Ry Nea = Mat Mops + Ripa: (2b) 


Hierin ist St,—=b, und sind ,p 4, Ntypa die Komponenten des Ruckes Jig, von B gegen A in nor- 
maler Richtung bzw. tangentialer Richtung. Fiir diese Komponenten gilt 


Nie, =e Xt—tro®*=é Xt—o* dgy, (3a) 
worin t den Einheitsvektor in Richtung », 4 darstellt, und wonach der Betrag von }t,p4 den Wert 


Rizgg HT é—o') = re — vg 4@" (3b) 


hat, ferner 
Nepa=—Iwet, 


Ripa = 3WEr=—d3vpyé- 


(3) 


(Abb. Ic). Fiir den Fall eines Extremums von ¢ reduziert sich, da dann é verschwindet, die Tangen- 
tialkomponente R,p4 auf Rizp4 = — ro. 


3. Bahnkomponenten des Ruckes. a) Reine Drehung. Bei reiner Drehung ist der Punkt A 
fest und fallt mit dem Drehpunkt 0 (Abb. 2) zusammen, d. h. es wird Si, gleich Null und die Kom- 
ponenten St, Jt, des Ruckes {tt eines bewegten Punktes stimmen dann unmittelbar mit den in 
(3a) bis (3c) wiedergegebenen Komponenten iiberein. Setzt man die Komponenten geometrisch 
zusammen, so erhellt sofort, daB 


1) der Ruck R proportional dem Abstand r des bewegten Punktes vom Drehpunkt ist, und 
2) daB der Ruck mit dem Strahl r = BO einen konstanten Winkel x bildet, welcher aus 


a R: ea eé— @3 
Alaa Biko as (4) 
folgt. Hat ¢ ein Extremum, so wird & = 0, also tg x = —q@?/3 «, d. h. es gilt tg (x + ~ =3tgy 


= 3¢/w*, wenn y der Winkel ist, den die Beschleunigung 6 mit dem Strahl BO bildet (vgl. 4a). 
Wenn zufallig ¢ = @° wird, so ist der Ruck nur normal gerichtet. 

Ist schlieBlich w konstant, d.h. ¢ = 0, so ist der Ruck nur noch tangential gerichtet, d. h. 
es gilt 


R= MR, =—tro® =—po?, (5) 


wobei der Kinheitsvektor { die Richtung von p hat, d.h. }t ist entgegengesetzt zu v gerichtet. 
Zeichnet man in diesem Fall die Geschwindigkeit v gleich dem Kurbelhalbmesser, d. h. setzt man 


Abb, 2. Der Ruck bei beschleunigter Kreisbewegung (E + 0). Abb. 3. Der Ruck bei nicht beschleunigter Kreisbewegung ( ¢ = 0). 


w gleich | (was bei kinematischen Konstruktionen immer gern gemacht wird), so ist auch » gleich 
dem Vektor tis OB’ (Abb. 3), der senkrecht zum Vektor t = OB steht. Infolgedessen stellt der 
Vektor t’’ = OB’, welcher wiederum senkrecht zu ry’ steht, die Beschleunigung § des Punktes B 
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(oder die Geschwindigkeit des Punktes B’) dar! und schli Blich 
th d é oe 
Vektor 1” = OB’” den Ruck 9 dar (vel. zit Jo OH ai ah ipa 


Ee on ) nicht konstant, aber € gerade gleich Null, d. h. hat m ein Extremum, so existieren auch nur 
1e l'angentialkomponenten mit den angegebenen Werten R, = v (é —o®). 


Yo, vy 


Abb. 4. Der Ruck am elliptisch bewegten Punkt (Kardankreispaar) 


Auch wenn der Ausschlagwinkel (z. B. in der Totlage einer Kurbelschwinge ein Extremum hat, also 
- @ momentan verschwindet, existieren nur Tangentialkomponenten des Ruckes. Fiir diesen gilt dann 
We". 

b) Sonderfall einer elliptischen Bahn. Bei einer Kreisbabn bewegt sich ein Punkt 
nach dem Gesetz x =acosa,y =asina mit a =wt(w =konst.) und a als Halbmesser des 
Kreises. Bei einer elliptischen Bahn (Abb. 4) mége diese entsprechend durch die Gleichungen 

x= bcoso, y =asina 

1 Der Hodograph ist wieder ein Kreis, und die Geschwindigkeit, mit welcher die Geschwindigkeitsspitzen 
wandern, ist die Beschleunigung. Ebenso wie die Beschleunigung bei konstanter Winkelgeschwindigkeit einen 


konstanten Betrag hat, ist auch der Ruckbetrag konstant, wie auch dort im Gegensatz zum Sprachgebrauch des 
taglichen Lebens. 
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gegeben sein, wobei a die groBe und b die kleine Halbachse darstellen (und wobei « als Kurbelwinkel 
in dem zugehérigen Kardangetriebe gedeutet werden kann). 
Da dx/dt = — bw sina und dy/dt = aq cos a gilt, so ist die Geschwindigkeit v durch 
v? = (b? sin? « + a® cos? x) w? = r7w? 

gegeben. Hierin ist r,; = OB’ der zur konjugierte Halbmesser; denn man erhalt diesen, wenn man 
in der Ellipsengleichung « durch « + 2/2 ersetzt. Fiir die Beschleunigungskomponenten folgt ahnlich 

d?x d’y ee 

: Fz ew sing , 

so daB die Beschleunigung durch 6 = — tq? oder durch b = tr, 0? dargestellt wird, wobei r, den zu 
OB’ =r, konjugierten Halbmesser bedeutet. SchlieBlich stellt Jt = — 1, @* = 1,03 den Ruck dar, 
wobei wieder r; zu r, konjugiert ist. Zeichnet man nun die Geschwindigkeiten so, dabia — lists 
so stellen die Vektoren t,, t, = — 1, t3 = — 1, die Geschwindigkeit, die Beschleunigung, den Ruck 
dar (Abb. 4). Bei dieser Zentralbewegung ist auch t, die Geschwindigkeit des Punktes B’ und 1, 
die Geschwindigkeit des Punktes B” oder die Beschleunigung des Punktes B’. 


Abb. 5. Die Einheitsvektoren an einer belicbigen Bahn. 


c) Beliebige Bahn. 1) Komponenten: Um bei einer beliebigen Bahn die Ruckkompo- 
nenten zu ermitteln, sei von der Entwicklung der Beschleunigungskomponenten ausgegangen. 
Nach Abb. 5a gilt fiir die Geschwindigkeit des Punktes 4 doch ») = t v = t ds/dt, wobei t den Ein- 
heitsvektor in Richtung der Bahntangente und t die Zeit bedeutet. Dann folgt fiir die Beschleunigung 


dv dt de 
din ah! 


oder da nach Abb. 5b doch dt =n dy ist mit n als Einheitsvektor in Richtung der Normale 
(zum Kriimmungsmittelpunkt hingerichtet) und mit dy als Winkel zwischen zwei unendlich be- 
nachbarten Normalen, so folgt unter Einfiihrung des Kriimmungshalbmessers 9 gemaB o dy = ds 
auch 


6= 


ON Bride v 
5 eon re o- 


also 


2 


b=th+n—=b, +h, (6) 


mit 
(b,| = b, = dvfde , 
|b,.| = v?/o 


in bekannter Weise. Die Differentiation von 6 aus (6) liefert dann den Ruck. Hierbei tritt der 
Vektor dn auf, und fiir diesen ergibt sich gemaB Abb. 5¢ doch 


dt dpe 
to) 
und 


dn v 
ge aad 
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Damit erhalt man 


db an v2 2vb 2 dt 
H ee | Uv UF v & db, 
a ab io | 0 alpen ia 
oder schlieBlich auch 
RSet tif RK, -ew- Ro (7a) 
wobei 
pee dbt es d*v v3 ES 1 /ds\3 
adie da ge: de ot (a) 2 oe) 
do v? v f/f d 
Fe ee gee, OY ee 80 
n 3 t 0 dt 2 0 i b, 0 dt (7c) 


ist. Wenn Q = konst., so ergeben sich die Werte wie beim Kreis bzw. die Komponenten it, g 4, St 


beim Satz von Euler. ig 


2) Die Flachengeschwindigkeit zweiter Art. Nun kann man den Wert fiir R, noch anders schrei- 


ben; denn es gilt auch! 
1 d /v® 
Be a | ee) 


wie eine Nachrechnung leicht erkennen laBt, und dem Klammerausdruck kommt eine besondere 
Bedeutung zu: Im Hodographen (Abb. 6) hat der Flachensektor d® zwischen den (unendlich) 
benachbarten Geschwindigkeiten ) und », =» + dv den Wert 
iL eeecas I) Ge 
— — 72 SS tS 
d® 5 Ui dy tay og 


Also gilt fir die Flachengeschwindigkeit zweiter Art? 


1@ Liters ds alae 
5. garde = Deore (8b) 


und somit auch? 


T° (8c) 


Ist also ® konstant, so ist R, immer gleich Null oder, wenn @ ein Extremum hat, momentan 


gleich Null (vgl. die folgende Ziff. 3). 
yh 


Tangente 


Abb. 6. Der Hodograph und die Flachen- Abb. 7. Der Ruck in Komponentendarstellung 
geschwindigkeit 2. Art. im festen und bewegten Koordinatensystem. 


3) Parameterdarstellung. Ist die Bahn eines Punktes in Parameterform, d.h. durch x = x(t) 
und y = y(t) gegeben, so sind die Komponenten des Ruckes in den Achsenrichtungen durch R, = 
x, R,=¥ gegeben und die Normal- bzw. Tangentialkomponenten folgen dann gemaB Abb. 7 zu 


bf 
Boat COS Oa R, sin a , 


R, =— R, sina + R, cosa 


1 Die Dimension von R ist cm/sec? bzw. Lange/(Zeit)’. 
2 Zur Unterscheidung von der tiblichen Flachengeschwindigkeit erster Art. 
3 Es bestehen keine Schwierigkeiten, den Ruck auch durch polare Komponenten darzustellen. 
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oder mit tga = y/x, also 


(os =e ee sino = > 
x? +9 v v 
auch zu ak 
Rot ee a ee (9a) 
: 2 +5? ue Vx? +y%? 
Da bekanntlich 9 = —. a se ist, hat die Flachengeschwindigkeit auch den Wert 
Ce aoe a 
$= F=s 5-7-9), (9b) 


und Einsetzen in (8c) liefert den gleichen Wert wie in (9a). 


In dem obigen Beispiel ittber die elliptische Bahn wird, wie leicht nachzurechnen, D = + ab w® 
= konst., also R, = 0 wie dort bereits gefunden. Ferner wird noch ganz einfach 


R, = —o° 0? sin? x + a? cos?a = — v0, 
wie auch dort angegeben. Bei einer hyperbolischen Bahn gemaB x = bWojwt, y =a Sinwt 
wird ebenfalls @ = konst., d. h. R, = 0 und beilaufig R, = vo”. 


Als Beispiel fiir das Verschwinden von R, sei das Gesetz 
Pian NAA Stein) 


genannt (Parabel); es verschwinden x und y, also R,und R,. Bedingung fiir das Verschwinden von 
R, ist y-y +x-x =0. 

4, Geometrische Orter. a) Bekannte Ortskurven. Wahrend zuletzt der Ruck auf der 
Bahn eines Punktes betrachtet wurde, soll nun der momentane Bewegungszustand der gesamten, 
zum Momentanpol P hin orientierten bewegten Ebene untersucht werden. Dieser ist zunachst 
durch die Bresseschen Kreise, d. h. durch Wendekreis K,, und Gleichenkreis K,, welche sich im 
Momentanpol P als singularem Punkt und im Beschleunigungspol J rechtwinklig schneiden (Abb. 8a), 
gegeben. Da der Momentanpol als bewegter Punkt die zum Wendepol hingerichtete Beschleunigung 
b, = Da? hat, so folgen fiir Normal- und Tangentialbeschleunigung (Abb. 8a) nach dem Satz von 
Euler die Werte 


b= —b, sn 9 -- To = ©" (D sm 7) = —@ (7) (10a) 
b, =— b, cosp +re=—Do*cosp +re; (10b) 
hierbei ist PA =r, PA, =w = Dsing mit D als Wendekreisdurchmesser und stellen w, ¢ die 
momentane Winkelgeschwindigkeit bzw. — Beschleunigung der Ebene dar. Die Normalbeschleu- 


nigung verschwindet fiir r = w, d. h. wenn A auf dem Wendekreis liegt; und verfolgt man einen 
Punkt auf dem Wendekreis, so hat dort 6, einen periodischen Verlauf gemaB 


b, =— Dw cosp + Desing = Dw + & sin (p — y;) , 


wobei tg py = w*/e ist, sodaB fiir py = yy demnach b, verschwindet (Beschleunigungspol J) und auf 
dem dazu senkrechten Strahl sein Extremum hat (Punkt Z,,) von der GréBe 


b= Dyotte=e WE. 


Die Tangentialbeschleunigung verschwindet, wenn r = g cos p mit g = Dw?2/e ist, also der Bahn- 
punkt auf dem Gleichenkreis liegt. Auf diesem hat wiederum b,, einen periodischen Verlauf gemaB 
b, = — got + e sin (p — yj); 

b,, verschwindet fiir p = py, (Beschleunigungspol) und hat auf dem dazu senkrechten Strahl (Punkt 
Z,) sein Extremum von der Grobe —b, = WG w?. Im Punkt W ist beilaufig b, = D ¢ und im PunktG 
entsprechend b, = gw?. 

In den oben angegebenen Formeln, z. B. in (7c), kommt auch der Ausdruck do/dt vor. Verschwin- 
det dieser momentan, d.h. hat @ einen stationaren Wert, so liegen alle solche Punkte auf der 


* Es wird also von der Bewegung des Massenpunktes in die Ebene iibergegangen zur ebenen Beweguag 
eines Korpers. 
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Kreisungspunktkurve k,, welche im allgemeinen dadurch definiert ist, da®B der betreffende 
Kriimmungskreis mit der Bahnkurve vier unendlich benachbarte Punkte gemeinsam hat. lis wer- 
den also fiir die weiteren Betrachtungen die Daten dieser Kurve sowie auch die Angelpunkt- 
kurve k, (dem Ort der zugehérigen Kriimmungsmittelpunkte) eine Rolle spielen. Die Kurven 
haben die Gleichungen? ? 


it 1 it 


Saas | i 
- Earn ene (Kurve k,) , (11a) 
1 1 1 4 
ieee Pare (Kurve k,) . (11b) 
a 


y$ (Polbahnnormale) 


-Wendekreis 


L 
Polbahntangente 


Gleichenkreis 


Abb. 8. Zur Euler-Savaryschen Formel: 


a) Beschleunigungen eines Punktes A der bewegten Ebene auferhalb des Wendekreises ; 
b) Punkt A der bewegten Ebene liegt innerhalb des Wendckreises 


1 Vel, 7B. R. Miiller, Einfuhrung in die theoretische Kinematik, Berlin 1932 oder R. Beyer, Kinematische 


Getrieb these, Berlin/Gottingen/Heidelberg 1953. 
n° 2 W. Meyer zur an Ing.-Arch. 17 (1949) S. 308 sowie 25 (1957) S. 140. 
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Hier ist m der Durchmesser der Kriimmungskreise von k, und k,, welche diese langs der Polbahn- 
normale in P beriihren, wahrend I bzw. |, der Durchmesser des Kriimmungskreises ist, welcher 
k, baw. k, lings der Polbahntangente in P berihrt, und ist ferner 


DACRE Cy Een pee hee 
gi? Eee: Oe er a 


i (12) 
Die Striche bedeuten Ableitungen nach dem Bogenelement der Rastpolkurve k,, und d? ist die 
momentane Winkeldrehung der Ebene, so daB 3! = 1/D und tr’ = 1/R, = Kriimmung der Rast- 
polkurve wird. Gemaf der Euler-Savaryschen Gleichung 


1 1 1 


r To D sin p 


gilt mit p = 90° und 1/R,, als Kriimmung der Gangpolkurve k, 


Pe arpa 1 te cee ; 
i aD wh ee 

; Cae (13) 
Be lt Antal 1,1 ue eal 2 
VR ONE ae 1 RR ee te 


Auf all diese Beziehungen wird gleich zuriickgegriffen. 


b) Komponenten des Ruckes in der bewegten Ebene. 1) Die Normalkomponente und 
der Tangentialkreis. Um die Normal- und Tangentialkomponenten des Ruckes zu finden, kénnte 
man 4bhnlich wie oben bei den Beschleunigungen den Satz von Euler anwenden, sofern der Ruck 
des Poles bekannt ware. Man kann jedoch auch unmittelbar von den Gleichungen (7b) und (7c) 
ausgehen, da ja bis auf do/dt zunachst alle EKinzelwerte, auch der Kriimmungshalbmesser o gemaB 
der Euler-Savaryschen Gleichung, bekannt sind. Es gilt nach Abb. 8b, wenn A innerhalb des Wende- 
kreises angenommen wird, 


r2 r2 


QE Se '—Dsnp=r sugore - ee 
Ferner ist v = r@, und gemaB (10b) betragt der erste Summand in R,, (7c) somit 
30 30 : 
Fy mr EXE) (15a) 
Fiir den Faktor dg/dt im zweiten Summanden von (7c) gilt zunachst 
do do ds , / 
Eo Laie OY =P ee: (15b) 


wenn ds hier das Bogenelement der Rastpolkurve, Striche Ableitungen nach s bedeuten und somit 
ds/dt die Polwechselgeschwindigkeit u = Dw darstellt. Differentiiert man (14) nach s, so erhalt man 
, r 


= Ging orgy 2" © sing +10" sing —rr' 8% —r I! g 008 9). 


? 


Nun ist nach Miiller} » 


r=—cosp, gp ==?_( 47/), 


und beachtet man ferner die Abkiirzungen (12), so folgt schlieBlich 


} 3rD i : 
r fo e -[- race) sin @ Cos °| 3 (15 c) 


e = (Dsing —r? 


Setzt man beilaufig die eckige Klammer gleich Null, so erhalt man (lla), wie zu erwarten. Fiihrt 
man nunmehr (15a) bis (15c) in (7c) fiir R, ein, so liefert dies geordnet 
| 
KRo==3 © [Dro® sin (cep ete w* cos (7 — D) — rel 


m 


' Vel. FuBnote | von Seite 59. 
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oder unter Beachtung der Definition fiir g, den Durchmesser des Gleichenkreises, und der ersten 


Gleichung (13) auch 


Re 2,250P pe .cosp _ 
3 30(D on re — Dat SP re (16a) 
oder 
= at OS Sh D ae: 
to. oe lin sin @ — jie As a (16b) 
mit der Abkirzung 
1 il 1 
a en (16c) 
L 
Polbahntangente 
Abb. 9. Verhiltnisse am Tangentialkreis und Normalkreis. 
Fiir einen weiter unten zu betrachtenden Sonderfall sei auch noch die Form 
1 A i 
Re 2 Dos Ee (D sin g —r) — D & Pe cal (16d) 
0 


angegeben. 

Jetzt kann die Frage nach dem Ort derjenigen Punkte gestellt werden, fiir welche die Normal- 
komponente des Ruckes verschwindet, also nur die Tangentialkomponente R, (s. u.) vorhanden ist. 
Nullsetzen von R,, in (16b) liefert dann 


r=slie sin /~ 7, 008 | = Bsinp— A cos, (17a) 
0 
mit 
Aes and | pee. (17b) 
I, m 
Hierfiir 1aBt sich auch schreiben 
r = T sin (y — q) (17c) 
mit 
enol \2 4d m* 
T=De\/(7) +(=,) und 8% = 37 (17d) 


Gleichung (17c) stellt also die Polargleichung eines Kreises, des Tangentialkreises Kpp vom 
Durchmesser T dar; sein Mittelpunkt liegt auf dem Strahl, welcher mit der Polbahntangente den 
Winkel y, + 2/2 bildet (Abb. 9). Da nach (17a) doch r(0) = — A und r(x/2) = B ist, schneidet der 
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Tangentialkreis auf der Polbahntangente die Strecke —.A und auf der Polbahnnormale die Strecke B | | 


ab. Der auf dem Strahl gy + 2/2 gelegene Punkt des Tangentialkreises sei mit W bezeichnet. 


2) Die Tangentialkomponente und der Normalkreis. Zur Ermittlung der Tangentialkomponente 
benutzen wir die Gleichungen (7b), (10b) und (14), d. h. es wird 


db, v3 d v3 


Peete, a 2 She 
fae oP ae (re — Dw* cos @) re 
oder mit 
d d ds d 
= ow ee 
auch 


3 
R,=Doer +ré—Do' cos p- D'’ —2 Dwe cos + D? wsin p- 9 ——(Dsin p—r)*. 


Beachtet man die oben angegebenen Werte fiir r’ und gy’, ferner, daB fiir den im Glied mit sin y | 
vorkommenden Faktor 2D — (# + 1’) D? auch D®(9’ — 1’) = — 3 D®/I, gemaB (12) geschrieben | 


werden kann, so erhalt man schlieBlich 


2 3 2 a3 
R, =r (¢é —o*) : = “cos Pp Se in Q- (18) | 


ii ly 


Diese Tangentialkomponente verschwindet auf dem Normalkreis, so daf auf diesem nur die | 


Normalkomponente R, auftritt. Er hat die Gleichung 


3 D? @® /cos sin x aa i 
Pas ( — a ol =Acosg + Bsingy = N cos (py — q) » (19a) 
also den Durchmesser 
3 D? w? 1 1 
ya = (19b) 


Sein Mittelpunkt liegt auf dem Strahl py — qo mit pp aus (17d), und er trifft die x- baw. y-Achse in 
den Strecken (Abb. 9) 
A= ale _— und B= se De 


6 o* Ke é—o> (Poe (1940 
c) Der Ruckpol. 1) Definition. Der Tangential- und der Normalkreis schneiden sich auBer 
im Momentanpol, dem singuliren Punkt, noch im Punkt Q, und zwar rechtwinklig (Abb. 9). Dort 
verschwinden R, und R,,, es ist R gleich Null, d. h. Q stellt den Ruckpol dar. Setzt man die Werter 
aus (19a) und (17c) einander gleich, d. h. setzt 


T sin (pg — Yo) = N cos (pg — Po) » 
so folgt mit pg — Py = Op die Richtung des Strahles PQ aus 


2 

tg 0g = tg (Pg —%) = a 5 se = oe . (20) 
Wahlt man nun Q als Bezugspunkt und wendet den Satz von Euler fiir den Ruck an, so haben 
wir (momentan), da Ry = 0, die gleichen Verhiltnisse wie bei der Drehung um einen festen Punkt: 
Die Ruckbetrage verhalten sich wie die Abstande der bewegten Punkte vom 
Ruckpol, und die Ruckvektoren bilden mit dem vom Ruckpol zu den bewegten 
Punkten gezogenen Strahlen den gleichen Winkel Yo =1/2— og. (Vgl. hierzu die 

spdter angegebenen Werte fiir R in den Punkten W, , Gy und P.) 


2) Verlauf von R, auf dem Tangentialkreis. Man kann dem Wert von R, aus (18) auch die Form 
R, = (é —o) [r — (A cos p + B sin ¢)| 
= (é—o%) [r — Neos (p — 9) 


geben, und setzt man hierin die Gleichung r = T' sin (py — gp) des Tangentialkreises ein, so folgt 
mit dem vom Strahl PG,, aus gerechneten Winkel o — Y — Po einfach 


(21) 


R, = (é —o) (T sino — N cos 0) = (é —o*) Wy Gy sin (o— og); (22) 
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da /T? + N? gleich der Strecke W7Gy ist (Abb. 9). Wie erwartet, verschwindet R, fir o = 09 


(Punkt Q); ferner hat R, sein Extremum R, = ( —) WrGy auf dem Strahl ¢ =o + a2 
also im Punkt Z,,, wobei PZ rallel QW, a i : at 
T r parallel QW > verlauft. Im Punkt W, gilt 


R, = (¢ —@*) W7Gy cos 09 = T (e—). 
3) Verlauf von R,, auf dem Normalkreis. Man kann dem Wert R,, aus (16a) auch die Form 
R, = 3 ¢ [(B sin y — A cos y) —r] = 36 [T sin py —q) —r] (23) 


geben, und setzt man hierin die Gleichung r = N cos (p — gy) des Normalkreises ein, so wird wie 
oben mit o = py — gp 


n = 3 &(T sino — N coso) = 3m WpGy sin (o — a9) . (24) 


Abb. 10. Der Ruck eines Punktes A der bewegten Ebene auferhalb des Normalkreises. 


Im Punkt Q(o = o9) verschwindet R,, und R, hat sein Extremum R, = 3w¢ W,Gy im Punkte Zy 
auf der Verlangerung des Strahles PZ, (Abb. 9). 


Im Punkt Gy gilt mito = 0 der Wert 


nm 


R =—3we-W,Gysinog =—3 Noe. 


d) Der Ruck des Momentalpoles. Wenn der Ruck }i,, des Momentanpoles P die Kompo- 
nenten R,,,und R,, in Richtung der Achsen x und y (Polbahntangente und Polbahnnormale) hat, 
so liefert die Anwendung des Satzes von Euler gemaB Abb. 10 die Komponenten! 


R,=—3ro@e+ R,,cosp + Ry sing, (25a) 
R, =r (¢ —o*) — R,, sng —K,, cos, (25b) 
und der Vergleich mit den Gleichungen (16b) und (18) zeigt, daB 
Day 3 D? w 
toe = a ha 5 ny = i ie (26a) 
sein muB. Der Ruck ft, hat also, da R,,/R,, = —|,/m* ist, die Richtung von PW, (Abb. 10) 
und den Betrag 
ea ei : 7 
R, = 3 D?o® ee Oe (26b) 


unter Beachtung von (17d). 


1 Man gehe in Abb. 10b von einem positiven, d. h, nach rechts gerichteten Wert Rpx aus. 
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e) Das Extremum der. Tangentialbeschleunigung. Nach (7b) ist db,/dt = 0, wenn 

R, = — v®/e, also auch nach (18) 
4 D2 3 3 D? 3 
r(é—%) 3 ae ae a) 

0 
ist. Bildet man unmittelbar db,/dt, wie bei der Entwicklung von R, angegeben, so tritt dabei der 
Ausdruck 1/D + 1/R, auf, welcher nach (13,2) aber 1/R,, bedeutet, und dann kann auch geschrie- 
ben werden 


sing + = (D sin y — r)? = 0 (27a) 


2 2 


3 D2 : 1 é ae 

— = COs Y — itis sin p + 7 [7 58 + D? sin? ®) =0 (27b) 
als geometrischer Ort r = r(p) derjenigen Punkte, welche momentan ein Extremum ihrer Tangen- 
tialbeschleunigung b, aufweisen.! Es handelt sich hierbei um eine Kurve vierter Ordnung, wie 
Einsetzen von cosy = «/r, sing =y/r und r? = a? + y? zeigt. Auf die Diskussion im einzelnen 
sei jetzt verzichtet, nur wird die Anwendung in Ziff. 4 der II. Mitteilung bei der zyklischen Bewe- 


gung erortert. 


Il. Relativbewegung 


5. Die Relativbeschleunigung. Hier erweist es sich als zweckmaBig, noch kurz auf die Entwick- 
lung der Sitze iiber Relativbewegung einzugehen: Die Lage des Punktes A einer bewegten Ebene 
sei durch den Vektor OA = a und die des Punktes B durch den Vektor OB = t angegeben, wobei 
O ein fester Punkt sei (Abb. 11). Dann ist die Lage des Punktes B gegeniiber A durch den der 
GréBe und Richtung nach veranderlichen Vektor 4B = 8 —t—a gekennzeichnet, und die Relativ- 
geschwindigkeit des Punktes B ist gegeben durch 

d8 (28) 
wobei wie itblich? das Zeichen ) auf die relative Anderung hinweisen 
soll. Ferner gilt fir die Fuhrungsgeschwindigkeit 


Dp =Dyg +X 8, (29) 
so dah 
Abb, 11. Die relative Lage des wird mit 
Punktes ae Saeed dem dg D3 Pm 
Breen ew (30b) 
Aligemein gilt also fiir die absolute und die relative Anderung eines Vektors & 
d3 dy le 
Fir die Absolutbeschleunigung des Punktes B gilt ferner 
ela ete Tae (32) 


worin 
Sane , : 
b= <a die Relativbeschleunigung, 
bp =b4—w?8 + € X 8 die Fiihrungsheschleunigung (Satz von Euler), 
6, =2@ X v, die Coriolis- oder Zusatzbeschleunigung 
darstellen mit € =@ = do/dt. 


6. Der Relativruck. Um nun den Ruck {t = db/dt zu finden, mu (32) unter Beachtung von 
(31) nochmals differentiiert werden. Dann folgt 


TR iene age wei Ce a ae 
a= (ae +2 X 84 = oF FOX 8)—2wes+Exs +Ex (FE +Hxs) 
= = D 
+2E Xd, +20 x (+O xv,), 


ae Punkte, welche ein Extremum der Geschwindigkeit v = ufwei legen j i 
Pires ahr g r@ aufweisen, liegen ja auf dem Gleichen- 


2 Vel. z. B. G. Hamel, Elementare Mechanik, Leipzig 1923 oder I. S infil ine di ; 
Mechanik, Berlin/Gottingen/Heidelberg 1954, ee aetna ie Kee 
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Beachtet man, da8 db,/dt = Hi, als Relativdruck und db,/dt = Kt, als Ruck des Punktes 4, ebenso 
Db, /Dt = b, geschrieben werden kann, ferner daB 
€xX @ X 3)=—en$ und © X @ Xb) =—o?bd, 
ist, so bleibt 
R= MR, + (8, +e x 8 —w* Dpy— 33) + (30xb,+3exv,—3?b,), (33a) 
oder, in ahnlichem Aufbau wie bei der Beschleunigung, 


R= R,+R- + RK, > (33b) 
worin §t der Absolutruck, St, der Relativruck = db,/dz, Mt, der Fiihrungsruck gemaB (26a) bis 
(26c) oder (33a) und 

Haus = 3(@ X b, + € Xv, —@* b,) (33¢) 
gesetzt sind. Der erste Summand von i,,, steht senkrecht 6, der zweite senkrecht », und der 
dritte ist b, entgegengerichtet. 

Die Bedingung ¢ = 0 fiir eine Extremum der Winkelbeschleunigung! der Fiihrung ergibt dann 
zeichnerische Verfahren, wie sie in der IT. Mitteilung an den Beispielen 7 und 8 gezeigt werden sollen. 
Betrachtet man beilaufig das bekannte Beispiel des sich auf einer drehenden Scheibe vom Mittel- 
punkt radial nach auBen bewegenden Punktes, so ergibt die Differentiation der Komponenten 


x =r(y) cosy, y =1r(~) sing nach der Zeit und die vektorielle Zusammensetzung der dritten 
Ableitungen den gleichen Aufbau fiir Jt wie in (33b). 


(Eingegangen am 21. April 1958.) 
Anschrift des Verfassers: Professor Dr. Ing. Walther Meyer zur Capellen, Aachen, Technische Hochschule. 


1 € tritt nur in Jt auf. 
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Uber die Grenzamplitude einer Klasse selbsterregter Schwingungen 
Von Hans Schubart 


1. Einleitung. K. Klotter hat in einer Arbeit, die vor kurzem in dieser Zeitschrift erschienen ist’, 
zusammen mit E. Kreyszig eine Klasse selbsterregter Schwingungen behandelt, die der Differential- 
gleichung 


i — (isn) |B —orgy ge] +f =9 (8%) (1.1) 


geniigen. Bekannte Transformationen (V = q’) fiihren mit V’ = dV/dq zu 
Vere BAL 08 Gy aay FG) a Oe (1.2) 


Als Lésung ergibt sich nach den ittblichen Regeln: 
Vig) == = (1(q) + F) (k eine Integrationskonstante) . (1.3) 


Dabei ist 
q 


sq) = exp |q(1—3 @)), 1a) = fo fo at. (1.4) 


Es handelt sich also im allgemeinen Fall um nicht elementar integrierbare Funktionen. 
Ein wesentliches Problem der Schwingungslehre stellt die Frage nach der GréBe der Grenzampli- 
tude q* in Abhangigkeit von der Riickstellkraft f(q) und von dem Koeffizienten « dar. Dabei wird 


q* bekanntlich durch 
(I(q*) =—I(C-9*).  (q* > 9)) (1.5) 


definiert. In der oben zitierten Arbeit ist gezeigt, daB q* = 3/0 fiir polynomiales f(q) gilt, wobei 


noch f(q) = 0 fiir g > 0 und f(q) S 0 fiir q <0 gefordert wird; auBerdem wird dort die Be- 
ziehung lim q*(I) = )/3/|a| fiir f(q) = C/q?'+! bewiesen. Weiter werden mit Hilfe transzendenter 
Gleichungen und nicht elementarer Funktionen obere Schranken fiir qg* gewonnen. 

In dieser Note sollen einige Aussagen tiber obere und untere Schranken von q* gemacht werden, 
wobei die Rickstellkraft 


f@) = S Gat} (CG, = 0) 


gesetzt ist. Es handelt sich dabei also um Riickstellkrafte, wie sie fiir die Schwingungslehre charak- 
teristisch sind. Bei der Berechnung der Schranken werden nur elementare Funktionen benutzt. 


2. Formulierung der mathematischen Aufgabe. Nach (1.5) muB fiir die Grenzamplitude g* > 0 


die Gleichung 


Fr Sy pe ( = e)| dt = — ft HOS [ (1 a e) dt (2.1) 


bestehen. Ersetzen wir rechts t durch (—?’) und beachten, daB 


ES eG G429,fgm=—fC®, f@+9 fir q+0 (2.2) 


ist, so erhalten wir 


also 
Ts f0 ae =f s(—e) pry ae. 


0 
1 K. Klotter u. E. Kreyszig, Ing.-Arch. 25 (1957) S. 389. 
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Schreiben wir nun wieder t statt ¢’ und setzen definitionsgemab 
s(t) — s(—1) = 2 Gin (ids a #)) 
so gilt, | 
q* 
2 f f(t) Sin 
0 


Wir untersuchen daher die Funktion 


Bt f a #) dt =0. (2.3) 


F(x) =) f(o) Sin ee(1 a) de ine x (2.4) 


Nach (2.2) hat mit r(x) = 6 x (1 —jo? «*) die Ableitung dF/dx = f(x) Gin (r(x)) fir «> 0 nur 
eine einzige Nullstelle & = V3/|a\. Da diese auBerdem einen Vorzeichenwechsel bedingt, so nimmt 
F(x) an der Stelle « = eine positives Maximum an, weil der Integrand > 0 im Intervall 0 <t <é 


ist. Aus (2.2) und (2.4) entnehmen wir schlieBlich, da bekanntlich Exponentialfunktion bzw. die 
Hyperbelfunktionen schneller fiir « —> co anwachsen als jede Potenz von x, 


lim F(x) = — oo. (2.5) 
Damit ist wegen F(€) > 0 und dF/dx < 0 fiir x > & die gesuchte Grenzamplitude q* die einzige 


positive Nullstelle unserer Funktion F(x) und es gilt qg* > & = //3’|a|, also das in Ziff. 1 zitierte 
Ergebnis. 


Im folgenden werden wir daher q* in der Form 6 & schreiben (6 > 1) und F(6&) = F(é) + 


[F'(6 €) — F(§)] abzuschatzen versuchen. Definitionsgema8 gilt dann fiir sign F'(6,€) = — 1 baw. 
sign F(6,é) = 1 die Ungleichung g* < 6, é baw. q* > 6, €. 
Um die Rechnung zu vereinfachen, beachten wir noch, das wegen Gin (— x) = — Ginx die 
Beziehung 
| 6é 
“Sf f@ Sin rH qi = Jf t) Gin r,(t) de (2.6) 
& 


mit r,(t) = 6 t (- «2? — 1) besteht. Damit kénnen wir uns bei der Bestimmung des Vorzeichens 
von F(6&) auf die Berechnung der Differenz zweier positiver GréBen F,(6 é) und F(é), die durch 


ji 6& ; & : 

F,— F = f f@ Gin r,t) de — f f(t) Stn r() dt (2.7) 
Gi 0 

definiert sind, beschranken. 


3. Obere und untere Schranken fiir F(é). Um eine obere Schranke von F(&) zu erhalten, er- 
setzen wir in (2.4) die Gin-Funktion durch ihr Maximum im Intervall 0 <t <&. Da aber Gin x 
fiir x > 0 monoton wachst, so muB das Maximum von r(é) fiir 0 <t < & bestimmt werden. 

Dieses wird nach bekannten Regeln fiir t = angenommen. Wir erhalten mit 


und (2.2) 


(3.1) 
i=0 

Soll eine untere Schranke fiir F'(€) angegeben werden, so miissen wir zunachst das Integrations- 

intervall so verkleinern, da®B in ihm keine Nullstelle der Gin-Funktion enthalten ist. Wir werden 


(Lan) & 
daher mit 0 <0, <1—o, <1 die GriBe F(€) durch ze f(t) Sin (r (Q) dt nach unten abschatzen. 


Im Falle der allgemein nach (2.2) gegebenen Rickaairat fg) wird sich spater (Ziff. 7) der Ansatz 
0: =0,=1/4 empfehlen, wahrend fiir f(q) = C q. f(q) = C@ sees = Cyq+ C,@ sich andere Ansatze 
als vorteilhafter erweisen. 

Eine wesentliche Schwierigkeit besteht darin, daB nichts iiber die GréBe von € bzw. von a vor- 
ausgesetzt werden soll. Es sei hier nur vermerkt, daB® sich u. U. unsere Formeln noch wesentlich 
verbessern lassen, wenn « weiteren Beschrankungen unterworfen wird. 


5* 
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In diesem Abschnitt behandeln wir nur die untere Schranke von I (é) fiir eine nach (2.2) allge- | 


mein gegebene Riickstellkraft und schatzen daher F'(é) durch 


E)> T(t) Gin r(t) dt (3.2) | 


&/4 
nach unten ab. Dabei werden wir f(q) in seine einzelnen Summanden zerlegen, was auch fiir (2.7) 


méglich ist, sodaB wir fiir unsere Betrachtung die GréBen C, nicht zu beriicksichtigen brauchen, | 


wegen ©in i ) < ©Ginr ie 7 :) erhalten wir mit r(E \= B &-15/64 


3 \2!+2 1 \2!+2 
Gy —() 


pé- = E2142 
64 21+2 


4 pit 2 Sin r(t) dt > Gin ( 


b/4 
g \Il41 red 
Gs) — (Gs) 


und damit nach (2.2) 


214 2 64 


F(é) > ae CG; 


4. Bestimmung einer oberen Schranke der Grenzamplitude q*. Um die Methoden, die uns weiter | 
unten eine bessere Bestimmung der Schranken von q* gestatten werden, zu erlautern und vorzu- | 


bereiten, soll hier zunachst der Satz S,: 


k 
Gao2e ee cur. 1} (9) =o CG) 
=0 


bewiesen werden. 

Hierzu geniigt es offenbar zu zeigen, daB in (2.7) die Ungleichung F(2 &) > F(E) richtig ist. 
Dazu benutzen wir (3.1) und zerlegen f(q) in seine einzelnen Summanden. Wir miissen die Un- 
gleichung 


28 é 
ie Gin r,(t) dt > if ?'+1 Gin r(t) dt (4.1) 


2€ é 
beweisen. Wir ersetzen dabei f durch f{ (0 <o <1) und wahleno = a Es wird dann sehr leicht 
é (l+o)& 
2 2g é 
fooa> food>fypQd mit g=F'1 Sins), y=?'*! Sinrl) (4.2) 
E 3/2€ 0 


zu zeigen sein. 
Da r,(t) fiir t > € monoton wachst, so ist 


Gin r,(t) > Gin Ae :)= Sin ie ) ee 2 E<¢1<2é 
und in (4.2) ergibt sich 


B 3 | E2142 ie abe 
= OR ies a jaca (Vial 
ie in(4 Bay oraea ie ;) 
Zum Beweis von Satz (S,) mu8 also nur noch zufolge (3.1) 
lee ae ATG E2142 B V3 15 ; 3 \2!+2 
srg Gin FF SE gin (215. eR GS 


erfillt sein. Da aber offenbar fiir alle 1] = 0,1,2,...,k die Beziehung 


3\2142 
yay ff 
l<2 (=) 


richtig ist, so haben wir (S,) vollstandig bewiesen. 


5. Verbesserung der oberen Schranke von q*. Wir miissen, um die in Satz (S,) gefundene obere 


one zu eos: einmal auf die GréBe Gin Ge a) aus (3.1) achten und zum anderen versuchen, 


fi aus (4.2) durch Of . “0 <y <o <1 méglichst giinstig nach unten abzuschatzen!. Als neue — und 
verbesserte—obere Sante von q* ergibt sich dann (1 + 9) &. Dabei wird in F, aus (2.7) Gin r,(t) 


‘1 Dje hier und in folgendem figurierende GréBe y ist nicht mit y aus (1. 2) bzw. (1.3) identisch, 


(3.3) 


eee (3.4) 
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durch Gin r,(1 + y) € nach unten abgeschatzt. Unsere Betrachtungen fiihren damit gleichzeitig, 


da fiir allgemeines € nicht Potenzen und Gin-Funktion verglichen werden kénnen, zu den beiden 
Ungleichungen 


fw eh 2 
Sin a S Gin r,[(1 + y) é], (5.1a) 
§ (L+e)é 
oaths? dese af lt 2 de. (5.1b) 
0 (l+y)é 


Wegen 


| ae 
néd+y) =P yy 4yo +2) 
folgt aus (5.la) y > 0,1547... Dabei kann y = 13/84 gesetzt werden. Mit diesem Wert ist die Un- 
gileichung erfiillt, wahrend y = 0,1545 dies nicht leistet. Der Wert y = 13/84 1aBt sich héchstens 
noch in der vierten Dezimalen verbessern. 
Nach (5.1) mu8 auBerdem ¢ so bestimmt werden, daf 


1) aie Sete Ce (5.2) 
erfillt ist. Da wir g (0 > y) méglichst klein machen wollen, ersetzen wir (5.2) durch 
1 Lt y\2t+2 
oo (+ o)2+1 _ (i + 20. os) 
Gleichung (5.3) muB fiir ] = 0,1,..., k gelten. Aus dieser Ungleichung ist unschwer zu erkennen, 


daf} fiir wachsendes | sich g immer mehr y néhert. 
Damit folgt zunachst eine etwas unscharfere Begrenzung der Granzamplitude nach oben, als 
sie in der oben zitierten Arbeit angegeben ist. 
Ks gilt der Satz S,: 
Ist f(q) = C q?'*?, so geniigt q*(l) der Ungleichung 


(E pyle < lim q*(I) < 1,1547€. 
l=o 


[| 


Soll dagegen (5.3) auch fiir ] = 0, 1,...,k gelten, so ergibt sich wegen 
l L+yzt2 1 1+y\2 
a ties) Saye tlre) 


1 Dee ey 
(1+)? ' \l+e 


eine fiir alle J giiltige Schranke von 9. Als Wurzel dieser quadratischen Gleichung erhalten wir 
Lte=siFO+ 
1+o 21,5276. 


aus 


oder (mit y = 13/84) 


Damit gilt der Satz S,: 

Hat die Riickstellkraft f(q) die Form (2.2), also insbesondere auch fiir f(q) = Cq, so ist q* > 
1, 5276. 

Ist Ss f(q) = C@, so muB die Gleichung (1 + 0) =1-+ (1 + y)* gelést werden. Die 
hieraus zu errechnende obere Schranke ergibt sich zu 1,2987 &. 

Die in Satz (S,) angegebene Schranke gilt naturgemaf auch fiir die spezielle Riickstellkraft 
f(q) = Coq + C, @. Wollte man diese noch verbessern, so geht in die Formel die GréBe von € ein. 
Ersetzen wir in (5.1b) die t-Potenz durch f(t) = Cyt + C,#, so geniigt x = (1 + @)? der quadra- 
tischen Gleichung 


Cy 2 


Die positive Wurzel x, von (5.4) liefert dann 1 + @ = |x, und mit diesem Wert 0 ist (5.1 a) sicher 
erfiillt; denn mit y = 13/84 bedingt (5.4) die Ungleichung «, > (1 + y)?, also auch (1 + @) > (1 + 7). 
Wir bemerken noch abschlieBend, da8 fiir f(q) = C q?'*? aus (5.3) 
fog (le 7) Cd) 
1 +o = exp( Seiaae 
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folgt,! und fassen diese Ergebnisse in dem folgenden Satz S, zusammen: 
a) Fir f(q) = Cq gilt q* <1,5276€, 
b) fir fq) = Gq? silt q= =< P2987 

logs (ae a euae 


; \)} 
ce) fir f@Qe G>T) gilt ye payee 


d) fir f(q)=Cyq+C,q@ — kann die in Satz (S;) angegebene obere Schranke von q* durch 
|x, € ersetzt werden, wenn x, die positive Wurzel der quadratischen Gleichung (5.4) ist. 


6. Berechnung einer unteren Schranke von q* in drei Sonderfiillen spezieller Riickstellkrifte. 
Wenn eine untere Schranke der Grenzamplitude bestimmt werden soll, so muf in (2.7) 


Fe) Sin de S ff Sine ra (6.1) 


gelten. Wir werden hier bestrebt sein, das links stehende Integral méglichst giimstig nach oben 
und die rechte Seite der Ungleichung méglichst giinstig nach unten abzuschitzen, ohne dabei das 
Zeichen > in (6.1) zu andern. Wir wollen zunachst den Fall 


S(Q) = God (6.2) 
behandeln. Die Abschatzung (3.2) laBt sich fiir diese Riickstellkraft durch eine andere ersetzen, die 
wesentlich zweckdienlicher ist. Da fiir <t = 0& sicher Gin r,(t) < Gin r,(6 &) ist, so gilt mit (6.2) 
zunachst 


Orem Ges = (8 —1) Gin [BE 6 (62 — D)]. (6.3) 


& 


(L—o2)é 
AnStelle von (3.2) werden fiir f(q) = Coq die GréBe F(é) durch =f tGinr(t)dtmit0 <e,<1—e, <1 


Q1 
nach unten abschatzen, dabei aber 0, und 9, so wahlen, daB dieses Integral immer noch gréBer bleibt 
als die in (6.3) angegebene obere Schranke fiir F,(0). Bei der eben geschilderten Abschatzung von 
F(é) nach unten, ist dann ©in r(t) durch 
A = Min (Gin r(o, €)/Gin r (1 — 0s)) 


zu ersetzen. Wegen 


£2 
EG) Aq 0s) Oa 
miissen wir 0, und Q, so bestimmen, daf nach (6.1) beide Ungleichungen 
(®—1) <—o,)'—@, (6.4a) 
s(@—1) <B mit B=Min(g, (1—e?)/e,(1—0) @—2)) (6.4b) 
gleichzeitig erfiillt sind. Damit ist (6.1) unter der Annahme (6.2) Geniige getan. 
Kine kurze Rechnung zeigt, daB (6.4) mit 9, = 0,5, 0, = 0,25, 6 = 1,125 erfiillt ist. Somit gilt 
gq. > 11256, 
Diese untere Schranke von q* laBt sich aber noch vergréBern (und damit verbessern), wenn wir 
0 <0, =e <@ +e =1—e@, <1 setzen und beachten, das damit aus (6.4a) und (6.4b) 
&@—l<2o0e+ 2, (6.5a) 
d(®—1)<B, mit B,=Min[(e—oyle +) —(@ + 2)I] (6.5 b) 
folgt und die Funktion f,(«) = « — «3 ihr Maximum im konkurrierenden Bereich fiir « — 1/3 an- 
nimmt. 
Wegen 6 > 1 mu8 nun go und ¢ so gewahlt werden, da 2 oe +e < B, gilt und endlich ist 


noch zu beriicksichtigen, unter welchen Bedingungen fie) SfAe+2e oder f, (0 +8) < ee 
W. ob B, = f,(0) oder B, =f, (0 + 8) ist. u a\@ ) <file) 


Hierbei erhalten wir oie Beziehungen: 


ee — 3a 302 . 
Fir 0 <@ Soe und Ves to RTESIS on 5 ge Huleent on 


l eer 
_ ca tO + V4=3 re . 
fir ye <Qo<1 md eS “SEES ist fie) 2 fileo +e). 


' Hier und in folgenden bedeutet log stets den (natiirlichen) Logarithmus zur Basis e. 
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Wenn wir dies bedenken, ergibt sich 6 = 1,1421 fir 0) = 0,43, € = 0,27, wobei es noch darauf an- 
kommt, die GréBe 20 ¢ + &2 moglichst geschickt zu variieren. Die angegebenen Werte fiir 0, ¢, 6 
lassen sich aber auch ohne diese Uberlegungen durch Rechnung bestatigen. 

Der hier angegebene Wert von 6 ist iibrigens, wie sich aus Betrachtungen tiber mégliche Extrem- 
werte der in (6.5) auftretenden Funktionen ergibt, héchstens noch in der dritten Stelle zu verbessern, 
wobei auch noch etwas iiber das Maximum aller nach (6.2) und (6.1) méglichen 6-Werte gesagt wer- 
den kénnte.t In Verbindung mit Satz (S,) a) haben wir somit den Satz (S,) a): 

Ist f(q) = Cq, so erfiillt die Grenzamplitude q die Ungleichung 1,1421 & < q* < 1,5276é. 

Wir betrachten nun einen weiteren Spezialfall, namlich 


fQM= Ce. (6.6) 
Analog zu (6.4) und (6.5) miissen hier jeweils gleichzeitig die beiden Ungleichungen 
a) #1 <(1—o)'—of, —d) (#1) SB, (6.7) 
bzw. 
a) #H—1<(@+e)'—o, —ib) 9 (° 1) SB, (6.8) 


erfillt sein. 

Die Schranken B bzw. B, sind von uns schon genauer analysiert worden; weiter ist zu beachten, 
daB fiir 6 > 1 stets dt— 1 > 6 (6? — 1) gilt. Nach (6.8) a) ist zu erwarten, da8 fiir 9 © « sich eine 
besonders giinstige Schranke fiir 64— 1 ergibt. Damit lassen sich allerdings keine besonders 
giinstigen Werte fiir B, erzielen. In der Nahe des Maximum der 6- Werte liegt 6, mit 6 > 1,1021 = 6,, 


= a + 3) errechnen 1a8t. Damit erhalten wir in Verbindung mit Satz 
S, a) und b) die Satze S; b): 

Ist f(q) = C q’, so gilt fiir die Grenzamplitude q* 1,1021& < q* < 1,2987€. 
und S; c): 

Ist f(q) = Cog? + C, gq’, so gilt fiir die Grenzamplitude q* 1,1021¢ < q* < 1,5276€. 

Kine dem Satz S, d) entsprechende Aussage aft sich bei unserem Verfahren mit Hilfe algebrai- 
scher Gleichungen vermutlich nicht gewinnen. 


was sich fiir@ =¢ = 


7. Generelle Bestimmung einer unteren Schranke der Grenzamplitude. Setzen wir fiir 1> 1 
f(Q) = Cg@'*, (7.1) 
so erkennen wir sofort, daf in der (6.4 a) bzw. (6.7) a) entsprechen Ungleichung auf der rechten Seite 
ein Faktor auftritt, der mit wachsendem / immer kleiner wird. Es liegt daher nahe fiir (7.1) in (2.7) 
6 = 1-4 G(l) mit lim Dl) = 0 zu setzen. Aus Griinden, die sogleich? evident werden, empfiehlt 
loco 
es sich, 

gaa SAS 1.2 
Ft ety (0 < o(l) <1), (7.2) 
zu wahlen, wobei dieses o von dem o in Ziff. 4 verschieden ist, und eine méglichst groBe obere 
Schranke fiir o zu bestimmen. Wir erhalten nunmehr analog den bisherigen Uberlegungen fiir F’, 

aus (2.7) entsprechend (6.3) 


8é - 
fete Gin n@) ar E21 = (1 4 o neta y Sin n(l at aes 7} g 


27+ 2(1+ 1) 
péo o \2 
Pl ( ergo (I +) | ce 


Ent see 
1k e- 

Bei der Abschatzung von F’(&) nach unten benutzen wir (3.4) und zerlegen f(q) in seine Summanden. 

Damit F(6 €) => 0 gilt, miissen wieder die zwei Ungleichungen 


a ee Fl = fc ae (7.4a) 


(e° — 1) Gin 


= 16 16 
(oy i | Oo 2 << 15 (7.4b) 
[oe r2d-+-1) = 64 


gleichzeitig erfiillt sein. 


1 Eine sehr ungenaue Abschatzung liefert (6.5) mit 6 (0? — 1)<fi (1/3) zu 6 < 1,16, wahrend wohl 6 S 
1,147 vermutet werden kann. 
2 Siehe die Abschatzung in (7.3). 
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Wir betrachten zunachst (7.4a). Wegen 


g\I+1 wees 9\l+1 il 5 /9\! Pd pa 
ays) Saf ors 


ist (7.4a) sicher fiir 


o(l) = los (1 sesh (nzeihae) 


erfiillt. s ee 
Da o(l) < o(2) fiir 1 => 2 und o(2) < 0,16, ist aber auch (7.4b) fiir dieses o(l) sicher richtig und es 
ergibt sich in Verbindung mit Satz S, c) der S, d): 
Ist f(q) = Ger +" 20 elt fur le2 
o(l) * 1 
( te aes oe qe<rt(é, 
wobei o(l) = log (1 “2 (5 )) und 7(1) in Satz S, c) definiert ist. 


SchlieBlich erhalten wir mit Satz S, und Satz S, d) den Satz S,: Hat die Riickstellkraft die 


allgemeine Form 
k 
fQ = = Cy git? (C, 9), 
=0 


so gilt fiir die Grenzamplitude q* 


* k 
( ns ier) Eecgt<1527%& mit o* =log i zm site| | 

8. Zusammenfassung. Wird die nichtlineare Differentialgleichung (1.1) einer Klasse selbst- 
erregter Schwingungen untersucht, so lassen sich untere und obere Schranken fiir die Grenzampli- 
tude q* angeben. Die Schranken kénnen fiir eine groBe Zahl solcher Riickstellkrafte die in der 
Praxis der Schwingungslehre auftreten, allein mit Hilfe bekannter elementarer Funktionen, gewonnen 
werden. 

Es gelingt, fiir die Riickstellkrafte der allgemeinen Form 


k 
F(q) ae i. (C, 
obere und untere Schranken fiir q* zu bestimmen (siehe Satz S, und Satz S,), wahrend fir 


fM9=Cq, f@=CHP,  fQ=—GQiraE 
sich diese Schranken noch verbessern lassen (siehe Satz S, und Satz S;). Die bereits explizit be- 


kannte untere Schranke € = )/3/|x| kann dabei fiir diese speziellen Riickstellkrafte nicht ganz un- 
wesentlich verbessert werden. Im Falle allgemeiner Rickstellkrafte ist dagegen die Verbesserung 
der unteren Schranke weniger erheblich. Sie kann jedoch fiir kleine |x| durchaus ins Gewicht fallen. 

Die oberen Schranken errechnen sich fiir die betrachteten Riickstellkrafte in allen Fallen ohne 
die Benutzung transzendenter Gleichungen und lassen sich explizit angeben. 


IV 


0) 


(Eingegangen am 10. Mai 1958.) 


Anschrift des Verfassers: Privatdozent Dr. H. Schubart, Karlsruhe/Baden, 
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Inh at stubersicht: Der Matrizenkalkiil: Grundbegriffe und einfache Rechenregeln. Das Matrizenprodukt. 
Die Kehrmatrix. Komplexe Matrizen. Lineare Abbildungen und Koordinatentransformationen. — Lineare 
Gleichungen: Der GauSsche Algorithmus. Lineare Abhangigkeit und Rang. Allgemeine lineare Gleichungs- 
systeme. Orthogonalsysteme. Polynommatrizen und ganzzahlige Matrizen. — Quadratische Formen 
nebst Anwendungen: Quadratische Formen. Einige Anwendungen des Matrizenkalkiils. — Das Eigenwert- 
problem: Eigenwerte und Eigenvektoren. Diagonalahnliche Matrizen. Symmetrische Matrizen. Normale 
Matrizen. Die Matrix 2(’9[. Abschatzungen. Eigenwerte spezieller Matrizen. — Struktur der Matrix: Mini- 
mumgleichung, Charakteristik und Klassifikation. Die Normalform. Hauptvektoren und Hauptvektor- 
ketten. Matrizenfunktionen und Matrizengleichungen. — Numerische Verfahren: Eigenwertaufgabe: 
Iterative Verfahren. Eigenwertaufgabe: Direkte Verfahren. Iterative Behandlung linearer Gleichungs- 
systeme. — Anwendungen: Matrizen in der Elektrotechnik. Anwendungen in der Statik. Ubertragungs- 
matrizen zur Behandlung elastomechanischer Aufgaben. Matrizen in der Schwingungstechnik. Systeme 
Iinearer Differentialgleichungen. — Sachverzeichnis. 


Die Matrizenrechnung als jiingster Zweig der Ingenieurmathematik wird hier dem Ingenieur in einer 
ihm zuganglichen Darstellung und mit ausdriicklichem Bezug auf technische Anwendung nahegebracht. 
Matrizen sind die ebenso kurze wie sinnfallige Formelsprache linearer Beziehungen, sie sind Formelzeichen 
und Rechenvorschrift zugleich. Sie beherrschen Theorie und Praxis linearer Gleichungssysteme, und sie 
fiihren zu dem auch fiir zahlreiche Anwendungen so bedeutsamen Eigenwertproblem. Diese beiden Auf- 
gaben, ihre Theorie und ihre numerische Behandlung, bilden den Hauptinhalt des Buches, wobei gréBter 
Wert auf eine enge Verbindung zwischen theoretischer Entwicklung und zahlenmaBiger Realisierung 
gelegt wird. Die Darstellung wird immer wieder durch Zahlenbeispiele, Abbildungen und Rechenschemata 
erganzt und belebt. Insbesondere die breit behandelten numerischen Verfahren, bei denen auch die neuere 
Entwicklung gebiihrend beriicksichtigt wird, sind durch ausfiihrliche Beispiele erlautert. Anwendungen 
in Elektrotechnik, Baustatik, Elastomechanik, Schwingungslehre und auf Systeme linearer Differential- 
gleichungen beschlieBen das Buch und regen zu eigener praktischer Arbeit mit dem so niitzlichen Hilfs- 
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Aus einer soeben in ,,Engineering“ erschienenen Besprechung : 


», His ist fraglich, ob es in England heute mehr als einen Mann gibt, der fahig ist, Professor Beyer’s Buch 
iiber kinematische Getriebeanalyse hinlinglich zu rezensieren (zu beurteilen). Das zeugt von einer nahezu 
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verzeichnis, in dem deutsche und amerikanische Abhandlungen aufgefiihrt sind, schlieBt Arbeiten ein, die 
erst im Friihjahr 1958 veréffentlicht wurden. AuBer den durchgearbeiteten Beispielen ist dem Buch eine 
groBe Anzahl von Aufgaben — zum Teil mit Fingerzeigen — zur Lésung beigegeben.“ 
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